문제풀이--

 

↑

program relaxation


 integer d
real a(700)
real b(2)

real c(3)
integer i, j, k
 
 a=0
 
open(unit=3,file='in.txt',status='old')

read(3,*) b(:)

read(3,*) c(:)

d=c(1)

read(3,*) a(2:d)


close(3)


do j=1,10000,1

a(1)=b(1)

a(d+1)=b(2)

do i=2,d,1


a(i)=a(i)+(a(i-1)+a(i+1)-2*a(i))/2.0

end do


end do


open(unit=4,file='red.txt',status='replace')



do i= 1,d,1

write(4,*) c(2)+i*((c(3)-c(2))/d) , a(i)



end do
 
 
 close(4)
 
 
 
 end program relaxation