X0 = [x0; v0]; X = ode(X0,0,t,DUFF); scf;plot(X(1,:), X(2,:)) scf;param3d(t, X(1,:), X(2,:)) I, the copyright holder of this work, hereby publish it under...
