% Usage: [y t] = euler1(f,a,b,ya,n) or y = euler(f,a,b,ya,n)
% Euler's method for initial value problems
%
% Input:
% f - Matlab inline function f(t,y)
% a,b - interval
% ya - initial condition
% n - number of subintervals (panels)
%
% Output:
% y - computed solution
% t - time steps
%
% Examples:
% [y t]=euler1(@myfunc,0,1,1,10);          here 'myfunc' is a user-defined function in M-file
% y=euler1(inline('sin(y*t)','t','y'),0,1,1,10);
% f=inline('sin(y)-cos(t)','t','y');
% y=euler1(f,0,1,10);

function [y t] = euler1(f,a,b,ya,n)
h = (b - a) / n;
y(1) = ya;
t(1) = a;
for i = 1 : n
    y(i+1) = y(i) + h * f(t(i),y(i));
    t(i+1) = t(i) + h;
end;