라플라스 방정식 풀기-유한차분법과 gauss seidel 법

 

라플라스 방정식의 해

 

program mt

integer i, j, k

real a(100,100)!----------- a(A,b,c) 에서 A,b,c, 는 각각 x,y, interation 횟수이다.

do i=2,99,1
do j=2,99,1

a(i,j)=0.0


end do
end do

!=================================초기조건

do i=1,100,1



a(1,i)=0.0
a(100,i)=50




end do

do i=2,99,1

a(i,1)=0.0
a(i,100)=50



end do


!=====================================================================================초기조건




!---------------------------초기조건(위)


!----------------------------------------- <프로그램 메인 부분 수치해법>




do k=2,9000,1
do i=2,99,1
do j=2,99,1



a(i,j)=(a(i-1,j)+a(i,j-1)+a(i+1,j)+a(i,j+1))/4.0

!----------------interation method gauss seidel-----------------------



end do


end do


end do

!===================================================================================================



do i=1,100,1



a(1,i)=0.0
a(100,i)=50




end do





do i=2,99,1




a(i,1)=0.0
a(i,100)=50


end do




!---------------------------------------프로그램 메인 부분 수치해법


!==========================================================================>


open(unit=1,file='laa.txt',status='replace')
!============= 파일에 값 저장 lap.txt 파일임

do i=1,100,1

do j=1,100,1


write(1,*) i,j,a(i,j)

end do

end do
close(1)


end program mt

4 모서리에서 경계조건이 주어진 라플라스 방정식 풀기 -포트란 코드