境界フラグをたてる

この時点で境界上にある節点と、そうでない節点を分けておくとあとが楽です。

int *bn=malloc(sizeof(int)*non);//border nodes 1:境界 0:内点
for(i=0;i<non;i++)
 bn[i]=i/(n+1)%n==0||i%(n+1)%n==0;