finitediff.lib
LIB "finitediff.lib"; list D="Ux","Uy","Ut","U"; list P="a","b"; list V="t","x","y"; setinitials(V,D,P); u(Ux); ==> [1] u(Ux,2,3,7); ==> [t^2*x^3*y^7] u(Uy)+u(Ut)-u(Ux); ==> [-1,1,1] u(U)*234-dx*dt*dy*3*u(Uy); ==> [0,(-3*dt*dx*dy),0,234]