genPendulum.m
ffd58092ca58master

genPendulum.m
View Options

			beep off
			%% para
			writeInFile = 0;
			nb_steps = 5000;%number of steps
			t0 = 0;%initial time
			t1 = 50;%final time
			dt = (t1-t0)/nb_steps;%time-step
			q0 = 0;%initial position
			% qDot0 = 10^(-6);%initial speed
			qDot0 = pi/1.6;%initial speed
			%% integrator
			a = @(t,q,qDot)pendulum(t,q,qDot);
			[t,q,qDot] = RK2(nb_steps,t0,t1,q0,qDot0,a);
			%% write in file
			if writeInFile == 1
			fileName = '../input/pendulum.csv';
			generatedData = [t, q];
			csvwrite(fileName,generatedData)
			end
			%% plot
			figure(1)
			plot(t,q)
			figure(2)
			plot(q,qDot)