% Script to solve assignment one
diary sol_one_data;
f = inline('1/(1-0.2*cos(y))','x','y');
yzero = 0;
xspan = [0 2*pi];
[x1 y1] = ode23(f,xspan,yzero);
[x2 y2] = ode45(f,xspan,yzero);
res_23 = norm(y1 - 0.2.*sin(y1) - x1)
res_45 = norm(y2 - 0.2.*sin(y2) - x2)
steps_23 = length(x1)
steps_45 = length(x2)
diary off;
