Il suffit d'utiliser l'opérateur logique ~ sur le vecteur en entrée :


>> a = [1 0 0 0 5 4 1 1 2 2 0 8 0 1 3 ]
>> [snz, nbnz] = nonzeros2(a)

snz =

1 6 1 2