在本論文中,我們使用 C++ 程式語言設計一個在二維空間上的 p-version 有限元素法程式。在此設計中,我們使用矩形元素計算以下偏微分方程的數值解[sum^{n}_{l=1} { - igtriangledown cdot (A^{kl} igtriangledown u_{l}) + B^{kl} cdot igtriangledown u_{l} + C^{kl} u_{l} - igtriangledown cdot (D^{kl} u_{l}) } = f^{k} , , k = 1 cdots n]這裡 n 是變數的個數,k 代表第幾條式子,因此上式有 n 個條式子。本論文將詳細說明程式架構與所使用的 hierarchical shape function In this paper, we design two-dimensional p-version finite element using rectangular elements by C++. In design, we solve the following partial differential equations by using rectangular elements[sum^{n}_{l=1} { - igtriangledown cdot (A^{kl} igtriangledown u_{l}) + B^{kl} cdot igtriangledown u_{l} + C^{kl} u_{l} - igtriangledown cdot (D^{kl} u_{l}) } = f^{k} , , k = 1 cdots n]where n is the number of variables, k stands for the number of equations, so we have n equations.This paper will show the struct of program and hierarchical shape function.