y0=[1; 0];  tspan=[0 10]; N=2; f=inline('[-y(2); y(1)]', 't', 'y')
y='[cos(t); sin(t)]';
for k=1:10
  [tt,u]=cranknic(f,tspan,y0,N);
  t=tt(end); e(k)=norm(u(:,end)-eval(y)); N=2*N;
end
semilogy(e); title('Erreur');