節点-座標対応表の作成

int i;//index
int n=3;//分割数
int non=(n+1)*(n+1);//節点数(number of nodes)
double **nbc=malloc(sizeof(double*)*non);//節点座標対応表(nodes by coodinate)
for(i=0;i<non;i++)
 nbc[i]=malloc(sizeof(double)*2);//x,y座標なので2つ
for(i=0;i<non;i++){
 nbc[i][0]=i/(n+1)*1./n;//x座標
 nbc[i][1]=i%(n+1)*1./n;//y座標
}