Интерполяциооный полином Лангранжа

Привет помогите пожалуйста, мне нужны комментарии к этому коду, а завтра сдать надо кровь из носу, буду благодарен. Нужны просто комментарии к этому коду, что для чего нужно. Это Интерполяционный полином Лангранжа

program interpoliacia
integer i,j,n,k,m
real x(10),y(10),matr(3,12),mas(11,11),lag(11),koef(11)
real polinom(10,10)
write(*,*) 'Vvedite kol-vo tochek n = '
read(*,*) n
do i=1,n,1
write(*,1) i
read(*,*) x(i)
write(*,2) x(i)
read(*,*) y(i)
enddo
1 format('Vvedite x[',i2,'] = ')
2 format('Vvedite y[',F8.3,'] = ')
do i=1,n,1
lag(i)=y(i)
k=1
do j=1,n,1
if (i.NE.j) then
lag(i)=lag(i)/(x(i)-x(j))
polinom(i,k)=x(j)
k=k+1
endif
enddo
enddo
polinom(n,n)=0
do i=1,12,1
do j=1,3,1
matr(j,i)=0
enddo
enddo
do i=1,11,1
do j=1,11,1
mas(i,j)=0
enddo
koef(i)=0
enddo
matr(1,2)=1
matr(2,1)=1
do i=1,n,1
m=2
do j=1,n-1,1
if (i.NE.j) then
matr(3,1)=-polinom(i,m)
if ((i.EQ.1).and.(j.EQ.2)) then
matr(1,3)=-polinom(i,m-1)
endif
if ((i.NE.1).and.(j.EQ.1)) then
matr(1,3)=-polinom(i,m-1)
endif
do k=1,j+1,1
matr(2,k+1)=matr(2,1)*matr(1,k+1)
matr(3,k+1)=matr(3,1)*matr(1,k+1)
enddo
do k=1,j+2,1
matr(1,k+2)=matr(2,k+2)+matr(3,k+1)
enddo
m=m+1
endif
enddo
do j=1,n+1,1
mas(i,j)=matr(1,j+1)*lag(i)
enddo
do k=1,3,1
do j=1,12,1
matr(k,j)=0
enddo
enddo
matr(1,2)=1
matr(2,1)=1
enddo
do i=1,n,1
do j=1,n,1
koef(i)=koef(i)+mas(j,i)
enddo
enddo
do i=1,n,1
write(*,3) koef(i),n-i
enddo
3 format(F8.3,'*x^',i2)
End