{[ promptMessage ]}

Bookmark it

{[ promptMessage ]}

BENG449 Problem 8 6

# BENG449 Problem 8 6 - 10:41 PM C\Documents and...

This preview shows pages 1–3. Sign up to view the full content.

4/10/08 10:41 PM C:\Documents and Settings\Alex Lemon\My...\BENG449 Problem 8 6.m 1 of 1 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % Alex Lemon % % BENG449 -- Biomedical Data Analysis % % Problem Set 8, Question 6 % % Due April 11, 2008 % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % The file derivative.m provides a function for calculating the following % differential equation: % dy/dt = a * t * (b * exp(- c * t) - y); y(0) = 0 % A proposed "analytical" solution to the problem is: % y(t) = a*b*exp(-1/2*a*t^2)*int(u*exp(1/2*a* u^2-c*u),y,0,t) % The integral is included because it cannot be evaluated numerically. % Using ode45, verify the solution. % First, choose some sample parameter values a = 1.5; b = 2.3; c = 6.1; % Next, calculate the values of the "analytical" solution: t = 0:0.01:10; % Compute values of the integrand in the expression for y i = t .* exp(1/2 * a * t .^ 2 - c * t); % Compute the cumulate integral values using cumtrapz I = cumtrapz(t,i); % Compute the values of the analytical solution y = a * b * exp(-1/2 * a * t .^ 2) .* I; % Compute a numerical solution using ode45:

This preview has intentionally blurred sections. Sign up to view the full version.

View Full Document
This is the end of the preview. Sign up to access the rest of the document.

{[ snackBarMessage ]}

### Page1 / 3

BENG449 Problem 8 6 - 10:41 PM C\Documents and...

This preview shows document pages 1 - 3. Sign up to view the full document.

View Full Document
Ask a homework question - tutors are online