[1] 授業単元：数値解析
[2] 問題文(含コード&リンク)：下記
[3] 環境 
　[3.1] OS： WindowsXP
　[3.2] Visual Studio .NET
　[3.3] 言語： C
[4] 期限： 4月末
[5] その他の制限： 基本的なCの知識しか習っていません(曖昧ですいません)。
副プログラム(サブルーチン)構造体、モジュールなども習っていません。
今持っている基本的な知識(if文、for文、配列)でこれを解けと言われました…。

ポアソン方程式の5点差分法による求解を考える。
ただし領域は、
0≦x≦1,  0≦y≦1,  とし、
境界条件は
φ[i,0]=φ[i,Ny]=0 (i=0,1,2,...,Nx)
φ[0,j]=φ[Nx,j]=0 (j=0,1,2,...,Ny)
とする。
またx方向の分割数をNxとし、y方向の分割数をNyとし
Nx*Nyが格子点数。
ポアソン方程式は
φ[i+1,j]-2φ[i,j]+φ[i-1,j]/(Δx)^2+φ[i,j+1]-2φ[i,j]+φ[i,j-1]/(Δy)^2=-13π^2sin(3πx)sin(2πy)
で与えられる。
上式をAφ＝bの形に変形し、その連立一次方程式をガウスの消去法で未知数φを求める
プログラムを作成せよ。