% Igor Yanovsky
% Teaching Assistant
% Math 151B

%#####################################################################
%
% This program solves the ordinary differential equation
% 
% dy/dt = f(y,t) y(t_0) = y_0
% 
% using Euler's method. Euler's method is implemented as a separate 
% routine where the function defining the ODE is passed in as a 
% function pointer. 
%
%
% Version 03/04/07
%#####################################################################

clear all

format long
% ODE Specification
%

t0 = 0.0;
tF = 2.0;
y0 = -1.0;

N = 100;     % number of timesteps

dt = (tF - t0)/N

tn = t0;
yn = y0;

for(i=1:N)
    ynp1 = yn + dt*f(yn,tn);
    yn = ynp1;
    tn = tn + dt;
end

yn
tn

clear all