real U2(100,100), U1(100,100), U0(100,100) ! open (10, file="laine.in") open (20,file="laine.out") read *,nt !***** initsialiseerimine alpha=0.05 do i=1,100 do j=1,100 U1(i,j)=0. U0(i,j)=0. enddo enddo !***** alg tingimused ***** do i=48,52 do j=48,52 U1(i,j)=10. U0(i,j)=10. enddo enddo !********************************************** do k=1,nt do i=2,99 do j=2,99 U2(i,j)=alpha*(U1(i+1,j)-4.*U1(i,j)+U1(i-1,j) *+U1(i,j+1)+U1(i,j-1))+2.*U1(i,j)-U0(i,j) enddo enddo !****** raja tingimused ********* do i=1,100 U1(i,100)=0. U0(i,1)=0. U1(100,i)=0. U0(1,i)=0. enddo !****** salvestamine *************** do i=1,100 do j=1,100 U0(i,j)=U1(i,j) U1(i,j)=U2(i,j) enddo enddo !****** kirjutamine ******** if((k/2)*2.eq.k)then write(20,100) U2 100 Format (100F12.5) endif enddo !*********************************************** stop end