This preview has intentionally blurred sections. Sign up to view the full version.
View Full DocumentThis preview has intentionally blurred sections. Sign up to view the full version.
View Full DocumentThis preview has intentionally blurred sections. Sign up to view the full version.
View Full Document
Unformatted text preview: tS enErrS 2.) dT/dx = q dq/dx = .15 * T function df = heat(x,y) df=[ y(2); .15*(y(1))]; function b = bal(qp) [x,y]=ode45(@heat,[0 10], [240 qp]); b = y(length(x),1)  150; >> fzero(@bal , 90) ans =86.5858 >> [x,y] = ode45(@heat, [0 10], [240 fzero(@bal,90)]); >> plot(x, y(:,1)) >> fzero(@bal,90) ans =90.6147 >> [x,y] = ode45(@heat,[0 10],[240 fzero(@bal, 90)]); 1 2 3 4 5 6 7 8 9 10 40 60 80 100 120 140 160 180 200 220 240 >> plot(x,y(:,1)) >> v = [2.15 1 0 0 0 0 0 0 0; 1 2.15 1 0 0 0 0 0 0; 0 1 2.15 1 0 0 0 0 0; 0 0 1 2.15 1 0 0 0 0; 0 0 0 1 2.15 1 0 0 0; 0 0 0 0 1 2.15 1 0 0; 0 0 0 0 0 1 2.15 1 0; 0 0 0 0 0 0 1 2.15 1; 0 0 0 0 0 0 0 1 2.15]; >> w = [240 0 0 0 0 0 0 0 150]'; >> T = v \ w T = 165.7573 116.3782 84.4558 65.2018 55.7281 54.6136 61.6911 78.0223 106.0569 T(0) = 240 T(10) = 150 1 2 3 4 5 6 7 8 9 10 40 60 80 100 120 140 160 180 200 220 240...
View
Full Document
 Spring '08
 KRACZEK
 Harshad number, tspan, energy error, Trajectory Error, larger ending time, Error StormerVerlet

Click to edit the document details