function f=Kepler_Eq(x,M,e)
f=x-e*sin(x)-M;