program laplace real dimension 15 15 V1 real dimension 15 15 V2 do 15

program laplace

real, dimension(15,15) :: V1
real, dimension(15,15) :: V2

do I=1,15
	do J=1,15
		V1(I,J) = 1000*random(L)
	end do
end do

do M=1,50
  do I=1,15
	do J=1,15
		V1(I,J) = (V1(I-1,J)+V1(I+1,J)+V1(I,J-1)+V1(I,J+1))/4
	end do
  end do
end do

do I=1,15
	do J=1,15
		write (*,*) V1(I,J) , " "
	end do
end do

end program laplace

function random(L)
integer :: L
L = mod(1366*L+150889,714025)
Random=Real(L)/Real(7142024)
return
end