function harmonic;

x0=[0.1;0.];
[t,x]=ode45(@rpp,[0 100], x0);
plot (x(:,1))function f=rpp(t,x)
f=[x(2);-2.*x(1)]
