유한차분법으로 열 방정식 풀기

 

 

열방정식의 marching solution

 

 

program heat

integer i,j,k !==========반복문 변수


real a(1000,1000)

a(1000,1)=0

a(1,1)=0


do i=2,500,1

a(i,1)=0.1*i



end do

do i=501,999,1

a(i,1)=(1000-i)*0.1

end do


do i=2,999,1


do j=1,100,1

a(i,j+1)=a(i,j)+0.01*(a(i+1,j)-2*a(i,j)+a(i-1,j))




end do





end do


do i=1,500,1

a(1,k+1)=a(1,k)

a(1000,k+1)=a(1000,k)



end do



open(unit=3,file='ther.txt',status='replace')



do i=1,1000,1

do j=1,100,1

write(3,*) i,j,a(i,j)
end do

end do

close(3)














end program heat