This preview shows page 1. Sign up to view the full content.
Unformatted text preview: h = Z(end, 2); if(h < -htol) dt = 0.5 * dt; if (dt < 1e-9) break end elseif(h > htol) z0(:) = Z(end, :); %change maybe t = t + dt; i = i + 1; time(i) = t; x(i) = Z(end, 1); H(i) = Z(end, 2); V(i) = Z(end3); gamma(i) = Z(4); end end e figure(1) plot(time, x) figure(2) plot(time, H) figure(3) plot(time, V) figure(4) plot(time, gamma)...
View Full Document
This note was uploaded on 09/18/2011 for the course ASE 167M taught by Professor Staff during the Spring '10 term at University of Texas at Austin.
- Spring '10