Problem set 15

Problem set 15 - >>

Info iconThis preview shows page 1. Sign up to view the full content.

View Full Document Right Arrow Icon
Problem 15-1 part 2 function dp = AmpConc(t,p) dp = zeros(3,1); % a column vector dp(1) = (-.5)*p(1)+(.477)*p(2); dp(2) = (1.25)*p(1)-(2.4125)*p(2)+(.569)*p(3); dp(3) = (.741)*p(2)-(7.34)*p(3); end >> options = odeset('AbsTol',1e-8,'RelTol',1e-5); >> [T,P] = ode45(@AmpConc,[0 10],[.04167 0 0],options);
Background image of page 1
This is the end of the preview. Sign up to access the rest of the document.

Unformatted text preview: >> plot(T,P(:,1),'-',T,P(:,2),'-.',T,P(:,3),'.') 1 2 3 4 5 6 7 8 9 10 0.005 0.01 0.015 0.02 0.025 0.03 0.035 0.04 0.045 Time (hours) Concentration (g/L) Ampicillin Concentration Over Time central ISF target...
View Full Document

Ask a homework question - tutors are online