neoclassics_calc_D113 Function

    function neoclassics_calc_D113()
        !! Calculates the integrated radial transport coefficients (index 3)
        !! 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) :: D113, neoclassics_calc_D113

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

        xi = roots
        wi = weights

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

        neoclassics_calc_D113 = D113
    end function neoclassics_calc_D113