# convert psolution to ODE;


toode(psolution(a_,b_,c_,d_,y_,x_,q_)) := y(q,x)-a*y(b,x)-c*y-d;
toode(psolution(c_,d_,y_,x_,q_)) := y(q,x)-c*y-d;