parameter n=100
	real tv(n,n,n),tu(n,n,n)
	open(10,file='soojus.in')
	open(20,file='soojus.dat')
	read (10,*) nt,alpha,iviz
	do i=1,n
	do j=1,n
	do k=1,n
	tv(i,j,k)=0.
	if (i.ge.35.and.i.le.65.and.j.ge.35.and.j.le.65
     *  .and.k.ge.35.and.k.le.65) then
	tv(i,j,k)=100.
	else
	tv(i,j,k)=0.
	endif
	tu(i,j,k)=0.
	enddo
	enddo
	enddo
	
	do l=1,nt
	
	do i=2,n-1
	do j=2,n-1
	do k=2,n-1
	if (i.ge.35.and.i.le.65.and.j.ge.35.and.j.le.65
     *  .and.k.ge.35.and.k.le.65) then
    	tu(i,j,k)=100.
	else
	tu(i,j,k)=tv(i,j,k)+alpha*(tv(i+1,j,k)+tv(i-1,j,k)+tv(i,j+1,k)+
     *  tv(i,j-1,k)+tv(i,j,k-1)+tv(i,j,k+1)-6.*tv(i,j,k))
	endif
	enddo
	enddo
	enddo
	
	do i=1,n
	do j=1,n
	do k=1,n
	tv(i,j,k)=tu(i,j,k)
	enddo
	enddo
	enddo
	
	if ((l/iviz)*iviz.eq.l) then
	write (20,100) tu
	endif
	
100	format (100 f15.6)
	
	enddo
	stop
	end