A simple method has been developed in order to get a contour plot in an element with given nodal values. The contour lines are obtained by solving a differential equation obtained from the shape function using symbolic computation. This method can be applied to any kind of shape (interpolation) function of a two-dimensional element.