function [u]=cnonestep(t,u,y,h,f,fn,varargin)
% CNONESTEP un passo del metodo di Crank-Nicolson
u = u + 0.5*h*(feval(f,t,y,varargin{:})+fn);
return