P21_11 - f c1 = cnet(1,1); c2 = cnet(2,1); c3 = cnet(3,1);...

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

View Full Document Right Arrow Icon
function P21_11 f tspan = [0:1:50]; t c0 = [1; 1; 0]; c [t, cnet] = ode45(@dcnetdt, tspan, c0); [ [t, cnet] = ode45(@dcnetdt, tspan, c0); [ x1 = cnet(:,1); x x2 = cnet(:,2); x x3 = cnet(:,3); x plot(t,x1, t,x2, t,x3) p function new1 = dcnetdt(t,cnet)
Background image of page 1
This is the end of the preview. Sign up to access the rest of the document.

Unformatted text preview: f c1 = cnet(1,1); c2 = cnet(2,1); c3 = cnet(3,1); c dc1dt = -0.013*c1 - 1000*c1*c3; d dc2dt = -2500*c2*c3; d dc3dt = -0.013*c1 - 1000*c1*c3 - 2500*c2*c3; d new1 = [dc1dt; dc2dt; dc3dt];...
View Full Document

This note was uploaded on 09/21/2011 for the course EGM 3344 taught by Professor Raphaelhaftka during the Spring '09 term at University of Florida.

Ask a homework question - tutors are online