neoclassics_calc_D111 Function

    function neoclassics_calc_D111()
        !! Calculates the integrated radial transport coefficients (index 1)
        !! It uses Gauss laguerre integration
        !! https://en.wikipedia.org/wiki/Gauss%E2%80%93Laguerre_quadrature
        !
        ! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
        use const_and_precisions, only: pi

        real(dp),dimension(4) :: D111, neoclassics_calc_D111

        real(dp),dimension(no_roots) :: xi,wi

        xi = roots
        wi = weights

        D111(1) = sum(2.0d0/sqrt(pi) * D11_mono(1,:) * xi**(1.0d0-0.5d0) * wi)
        D111(2) = sum(2.0d0/sqrt(pi) * D11_mono(2,:) * xi**(1.0d0-0.5d0) * wi)
        D111(3) = sum(2.0d0/sqrt(pi) * D11_mono(3,:) * xi**(1.0d0-0.5d0) * wi)
        D111(4) = sum(2.0d0/sqrt(pi) * D11_mono(4,:) * xi**(1.0d0-0.5d0) * wi)

        neoclassics_calc_D111 = D111

    end function neoclassics_calc_D111