This preview shows pages 1–2. Sign up to view the full content.
This preview has intentionally blurred sections. Sign up to view the full version.View Full Document
Unformatted text preview: ... rinit = sqrt(512/8)=8 % r/x = rinit/xinit = 8/8=1 (conical shape) % r(t) = x(t) % A(t) = pi * r^2(t) % A(t) = pi* (x(t))^2 % gfun = inline('-.6*pi*(.1)^2*sqrt(-2*-32.17)*sqrt(x)/(pi*x^2)','t','x'); g % changed final time to 100000000 mins so it will loop a while and then % added the break for the loop within the rkmid function to create % rkmidbreak % fourteenn = rkmidb(gfun,t0,100000000,xinit,h2); [m,n]=size(fourteenn); for i = 1:m if(mod(i-1,3)==0) fourteen((i-1)/3+1,1) = (i-1)/3; fourteen((i-1)/3+1,2) = fourteenn(i,1); end; end; fourteen(27,1)=25.20; fourteen(27,2)=fourteenn(25*3+2,1); f fourteen...
View Full Document
This note was uploaded on 09/19/2009 for the course MATH numerical taught by Professor Ford during the Spring '09 term at FAU.
- Spring '09