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