2次元のﾎﾟｱｿﾝ方程式

方程式 d2f/dx2 + d2f/dy2 = c

境界値　　f = 0

範囲 　 x= -1 〜 1 , y= 0 〜 1

定数 c= -1 (x<0) , c = 0 (x=0) , c = 1 (x>0)

離散式 -4*f(i,j)+f(i-1,j)+f(i+1,j)+f(i,j-1)+f(i,j+1)= dx*dx*c


これをvisual Fortranを使ってやるのですが全くわかりません。。。
誰か出来る人いませんか？？？？
あとDOSで開いた後立体的に見たいんですが、可能ですか？？？？
よかったら貼り付けてください！！！