|
D.14.1.3 scheme
Procedure from library findiff.lib (see findiff_lib).
- Usage:
- scheme([v1,..,vn]); v1,..,vn of type vector
- Return:
- type poly; performs substitutions by the means of Groebner Basis computation of the module generated by the input vectors, then intersects the intermediate result with the suitable component in order to get a finite difference scheme;
- Note:
- works only for a single pde, for the case of a system use matrixsystem
Example:
| LIB "findiff.lib";
list D="Ux","Ut","U";
list P="a";
list V="t","x";
setinitials(V,D,P);
def s1=scheme(u(Ut)+a*u(Ux),backward(Ux,U,x),forward(Ut,U,t));
s1;
==> (-a*dt+dx)/(dx)*x+(a*dt)/(dx)
|
|