function f = kepler ( x )
f = x - 3.0 * sin ( x ) - 1.0;
