Computer generation of difference approximations

An algorithmic approach, based on the method of undetermined coefficients, for generating difference approximation to partial differential operators is presented here. The approach is based on inverting Vandermonde matrices of lower order. Consequently, the computational time required to obtain the symbolic coefficients is very small in all cases. The algorithm is demonstrated by generating multi-parameters approximations for several problems.