程序名稱
　　　 shapef(s,t,xl,xsj,shp)――四結點四邊形等參元插值函數的計算
　 功　能
　　　本程序用以計算四結點四邊形等參元的Jacobi矩陣，Jacobi行列式的值及導數之間的變換。
　 使用說明
　　　 子程序語句 subroutine shapef(s,t,xl,xsj,shp)
　　　 參數說明
　　　　◎ s，t 實變量，輸入參數，自然坐標下所求點的座標值。
　　　　◎ xl 實變量，輸入參數，結點坐標。其中xl(1,i)表示x方向坐標，xl(2,i)表示y方向坐標。
　　　　◎ xsj 實變量，輸出參數，所求點Jacobi行列式的值。
　　　　◎ shp 實變量，輸出參數。其中shp（1，i）是插值函數對x的偏導，shp（2，i）是插值函數對y的偏導，
　　　　　 shp（3，i）是插值函數。
　 方法簡介
　　　四結點四邊形自然坐標下的插值函數：
　　　　　　　　 　　　
　　　由等參元的概念得：
　　　則導數之間的變換為：
　　　　　　　　
　 程序（提供下載）