BENG449 Problem 8 6

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

Info iconThis preview shows pages 1–2. Sign up to view the full content.

View Full Document Right Arrow Icon

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

View Full DocumentRight Arrow Icon
This is the end of the preview. Sign up to access the rest of the document.

Unformatted text preview: 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...
View Full Document

Page1 / 3

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

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

View Full Document Right Arrow Icon
Ask a homework question - tutors are online