</pre><pre name="code" class="cpp">/*函数名称:三角单元,线性插值,计算单元系数矩阵ke[3][3]////////////////////////
*xco[3]: double,input,store the x coordination,xi,xj,xm
*yco[3]: double,input,store the y coordination,yi,yj,ym
*ke[3][3]: double,output,store the element coefficient matrix
*/
void uke1 (double xco[], double yco[], double ke[][3])
{ int ce=0;
double a[3], b[3];
a[0] = yco[1] - yco[2];
a[1] = yco[2] - yco[......
阅读全文