real tv(100,100),tu(100,100),alpha data tv,tu/20000*0./ read *, nt open (10,file="s_2d_t=0.out") alpha=0.05 c******* alg tingimused *********************** do i=1,100 tv(i,1)=100. enddo c********************************************** do it=1,nt do i=2,99 do j=2,99 sum=alpha*(tv(i+1,j)-2*tv(i,j)+tv(i-1,j)+tv(i,j+1) * -2*tv(i,j)+tv(i,j-1)) tu(i,j)=tv(i,j)+sum enddo enddo c****** raja tingimused ********************** do i=1,100 tu(i,1)=100. enddo do i=2,100 tu(1,i)=0. tu(i,100)=0. tu(100,i)=0. enddo c*********************************************** do i=1,100 do j=1,100 tv(i,j)=tu(i,j) enddo enddo if((it/1000)*1000.eq.k)then write(10,100) tu 100 format(100e15.5) endif enddo stop end