%pendulum
function dxdt=pendulum(t, x)
L=1; %ingahossz
g=9.8; %gravitáció
theta=x(1); %kitérés
gamma=x(2); %sebesség

dtheta=gamma; %egyenletek
dgamma=-(g/L)*sin(theta);

dxdt=zeros(2,1);
dxdt(1)=dtheta; %output
dxdt(2)=dgamma;
end
