compute_trinomial_coefficients Subroutine

subroutine compute_trinomial_coefficients()

    integer :: i, j, k

    do i = 0, MAX_ORDER
        do j = 0, MAX_ORDER
            do k = 0, MAX_ORDER
                trinomial(i,j,k) = binomial(i,i) * binomial(i+j,j) * binomial(i+j+k,k)
            end do
        end do
    end do

end subroutine compute_trinomial_coefficients