1;
y0=[1; 0];  
tspan=[0 10]; 
f=inline('[-y(2); y(1)]', 't', 'y');
y='[cos(t); sin(t)]';

N=2;
for k=1:5
 
% call good
disp('get w and info')
[tt,u]=cranknic_good(f,tspan,y0,N);
t=tt(end); 
e_good(k)=norm(u(:,end)-eval(y)); 
disp('get w and info OK')

% call bad
disp('get only w')
[tt,u]=cranknic_bad(f,tspan,y0,N);
t=tt(end); 
e_bad(k)=norm(u(:,end)-eval(y));
disp('get only w OK')


N=2*N
end