function theta=reversefit(P,y,theta0);

TOL=1e-8;
F='Func1(x,P,y)';
FP='FFunc1(x,P)';
theta=Newton(F,FP,theta0,y,P,TOL);