function dzdt = A1_P1_ode_2025(t,z,M,k,c,A,w,L0,g)
% computes the derivative of the states z

if w~=0 % Check if the input force is harmonic
    F = A*sin(w*t);
else
    F = A;
end

dzdt(1,1) = z(2);
dzdt(2,1) = (F-M*g-c*z(2)-k*(z(1)-L0))/M;
end
