Equations of Motion for Linear and Quadratic Models

> IC := { x(0)=x0, v(0)=0 };

> EOM[lin] := { diff(x(t),t) = v(t),

> m*diff(v(t),t) = -m*g - k*v(t) };

> EOM[quad] := { diff(x(t),t) = v(t),

> m*diff(v(t),t) = -m*g + k*v(t)^2 };

>