No BSD License
function yh = euler(f, x, y, h) k1 = h*feval(f, x, y); yh = y + k1;
Contact us at files@mathworks.com