помогите сделать матрицу на фортане

Program TEST
    integer, parameter:: n=30
    integer::a(n,n)
    real(8)::p=1.0
    do i=1,n
        do j=1,n
            if (i==j) then
                a(i,j)=n-i+1
                p=p*a(i,j)
            else
                a(i,j)=0
            endif
        end do
    end do
    print *,p
end