%  function y=fib10(a,b,p)  finds first zero mod p  (<1000000);
function y=fib10(a,b,p)
if nargin==2, p=b; b=a; a=1; end
if nargin==1, p=a; a=1; b=1; end
disp('  Data:')
disp([a b p])
disp('________')
u=mod(a,p); v=mod(b,p);
if u==0, y=1; return, end
if v==0, y=2; return, end
limit=min(2222222, 2*p);
u=mod(a,p); v=mod(b,p);
for n=3:limit
w=u+v; u=v; v=mod(w,p);
% if mod(n,50000)==0, disp(n), end
if v==0, y=n; return, end
end
y='    No  zeros !';


