#include #include int main(int argc, char *argv[]) { if (argc!=3) { fprintf(stderr, "%s m n\n", argv[0]); return 0; } int m = atoi(argv[1]), n = atoi(argv[2]); int numElements = 2*m*n; int numNodes = (m+1)*(n+1); int numInteriorNodes = (m-1)*(n-1); int element[numElements][3]; double node[numNodes][2]; int p[numNodes], q[numNodes]; int cInt = 0; int cBound = numNodes-1; double hx = 1./m, hy = 1./n; int i, j, cNode, cElement=0; for (i=0; i<=m; ++i) { for (j=0; j<=n; ++j) { cNode = j+i*(m+1); node[cNode][0] = i*hx; node[cNode][1] = j*hy; if ((i>0)&&(i0)&&(j