function dzdt = EX3_NLode_2025(t,z,M,m,k,c,cd1,cd2,cd3,A,T)
 % computes the derivative of the states z
    F = 0;
    P = ; % COMPLETE
    if t>T/2
        P=; % COMPLETE
    end

    dzdt(1,1) = ; % COMPLETE
    dzdt(2,1) = ; % COMPLETE
    dzdt(3,1) = ; % COMPLETE
    dzdt(4,1) = ; % COMPLETE

end