R2096/eulerProg.m
b14423f6276bmaster

eulerProg.m
View Options

			function [t,q]=eulerProg(n,t0,t1,q0,qDot0)
			h=(t1-t0)/n;
			t = t0:h:t1;
			q = ones(size(t));
			qDot = ones(size(t));
			q(1)=q0;
			qDot(1) = qDot0;
			for i=1:n-1
			qDot(i+1)=qDot(i)+h*pendulum(t(i),q(i));
			q(i+1) = q(i) + h*qDot(i);
			end;