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