System of PDE's Satisfied by
We compute a system of PDE's for
>
by first computing an ODE with respect to the variable by using gfun[holexprtodiffeq] , and next encoding that does not depend on .
> holexprtodiffeq(f,y(x));