【正文】
j=2 do i=2,m1 w(i,j)=a(i,j)+s*(1s)/2*(a(i1,j)+a(i+1,j)+a(i,j1)+a(i,j+1)4*a(i *,j))+s*s/4*(a(i1,j1)+a(i1,j+1)+a(i+1,j1)+a(i+1,j+1)4*a(i,j)) end do j=n1 do i=2,m1 w(i,j)=a(i,j)+s*(1s)/2*(a(i1,j)+a(i+1,j)+a(i,j1)+a(i,j+1)4*a(i *,j))+s*s/4*(a(i1,j1)+a(i1,j+1)+a(i+1,j1)+a(i+1,j+1)4*a(i,j)) end do i=2 do j=2,n1 w(i,j)=a(i,j)+s*(1s)/2*(a(i1,j)+a(i+1,j)+a(i,j1)+a(i,j+1)4*a(i *,j))+s*s/4*(a(i1,j1)+a(i1,j+1)+a(i+1,j1)+a(i+1,j+1)4*a(i,j)) end do i=m1 do j=2,n1 w(i,j)=a(i,j)+s*(1s)/2*(a(i1,j)+a(i+1,j)+a(i,j1)+a(i,j+1)4*a(i *,j))+s*s/4*(a(i1,j1)+a(i1,j+1)+a(i+1,j1)+a(i+1,j+1)4*a(i,j)) end do do i=1,m do j=1,n a(i,j)=w(i,j) end do end do return end