function [y,t]=mysolve(A,y0,T,N)
%
%
%

h=T/N;
n=length(y0);
y=zeros(n,N+1);
if nargout>1
    t=linspace(0,T,N+1);
end
y(:,1)=y0;
for n=1:N
    y(:,n+1)=cnstep(A,y(:,n),h);
end