MAT 275 LAB 3.pdf

LAB 3 - Your Name - MAT 275

Exercise 1

Part (a)

Define ODE function for your version of the lab.f

f =@( t , y ) 0.7* y;

Define vector of time-values over the interval in your version of the lab, to compute analytical solution vector.t

t = linspace (0 ,3 ,100); y = -4* exp (0.7* t );

Create vector of analytical solution values at corresponding t values.

[ t60 , y60 ]= euler (f ,[0 ,3] , -4 ,60);

Solve IVP numerically using forward Euler's method with Nhuge timesteps (use variable names as instructed).

[ t6000 , y6000 ]= euler (f ,[0 ,3] , -4 ,6000);

In the following steps, we define error as exact - numerical solution value at the last time step

%eN=y(end)-yN(end)

Compute numerical solution error at the last time step for forward Euler with Nsmall timesteps (use variable

names as instructed).

eNsmall= y(end) - y6(end);

Compute numerical solution error at the last time step for forward Euler with Nmed timesteps (use variable

names as instructed).

eNmed= y(end) - y60(end);