function y=cnstep(A,y,h)

I=eye(2);
z=(I+0.5*h*A)*y;
y=(I-0.5*h*A)\z;

