function [u]=beonestep(t,u,y,h,f,fn,varargin) % BEONESTEP un passo del metodo di Eulero all'indietro u = u + h*feval(f,t,y,varargin{:}); return