! Determinantul LOTKIN = (1/ delta) *(-1)^(n-1)
function Delta(n)
intent(in):: n

delta =1


do i=1,n-1
i2 =i*2

delta =delta*comb(i2,i-1)*comb(i2,i)*(i2+1)
enddo

!delta =delta*(-1)**(n-1)

end

