			% It calculates ODE using Runge-Kutta 2th order method
			function [t,q,qDot]=RK2(n,t0,t1,q0,qDot0,a)
			dt = (t1-t0)/n;
			t = (t0:dt:t1)';
			q = ones(size(t));
			qDot = ones(size(t));
			q(1)=q0;
			qDot(1) = qDot0;

			for i=1:n
			% intermediate steps
			q_1 = dt*qDot(i);
			qDot_1 = dt*a(t(i),q(i),qDot(i));
			q_2 = q_1 + (dt/2)*qDot_1;
			qDot_2 = dt*a(t(i) + dt/2, q(i) + q_1/2,qDot(i) + qDot_1/2);

			q(i+1) = q(i) + q_2;
			qDot(i+1) = qDot(i) + qDot_2;
			end

R2096/generate/RK2.m
ffd58092ca58master

RK2.m
View Options