> restart;> pendphaseplot := proc (c, k, ui, uf, n1, vi, vf, n2, ti, tf, n, vu, vv) local s, i, a, j, b, sys, fcns, q, p; end proc;> s := {};>> for i from 0 to n1 do;> a := ui+i*(uf-ui)/n1;> for j from 0 to n2 do;> b := vi+j*(vf-vi)/n2;> sys := diff(u(t), t) = v(t), diff(v(.
