% dati epsilon ed eta calcola lo zero della funzione di keplero
a=eta-abs(epsilon)
b=eta+abs(epsilon)
fa=keplero(a,eta,epsilon)
fb=keplero(b,eta,epsilon)
numpassi=0;
while (b-a>eps) 
    passo2;
    numpassi=numpassi+1;
end
numpassi
c
