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