%2D1240 Numeriska metoder gk II
% Sergej Gratchev, e99_sgr@e.kth.se
% Mikael Simonsson, micke@monter.com
% December 2000
%
% Delfunktion för labb 4, uppgift 1: Futten

function f=fnewton(u, c)
r=[u(1) u(2)];
fi=[u(3) u(4)];

G = c(1);
R=c(2);
alfa=c(3);
g=c(4);

f=[r(2), G*cos(alfa)-g*(R/r(1))^2+r(1)*fi(2)^2, fi(2), (G*sin(alfa)-2*r(2)*fi(2))/r(1)];