Rules automatically generating the classical shape functions and finit
e difference patterns are developed. Finite difference solutions of La
place equation, Fourier equation, and the classical second-order wave
equation are demonstrated by using Mathematica.