function dzdt = A1_P2_ode_2025(t,z,M,k,c,L0,l0,alpha,beta)
% computes the derivative of the states z
g = 10;

dzdt(1,1) = z(2);
dzdt(2,1) = (M+alpha.^2.*beta.^2.*M.*sin(beta.*z(1)).^2).^(-1).*(alpha.*beta.* ...
  g.*M.*sin(beta.*z(1))+alpha.*beta.*k.*((-1).* ...
  alpha+(-1).*l0+L0+alpha.*cos(beta.*z(1))).*sin(beta.*z(1))+(-1).* ...
  alpha.^2.*beta.^2.*c.*z(2).*sin(beta.*z(1)).^2+(-1/2).*alpha.^2.* ...
  beta.^3.*M.*z(2).^2.*sin(2.*beta.*z(1)));
end

