neoclassics_calc_D112 Function

    function neoclassics_calc_D112()
        !! Calculates the integrated radial transport coefficients (index 2)
        !! 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) :: D112, neoclassics_calc_D112

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

        xi = roots
        wi = weights

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

        neoclassics_calc_D112 = D112

    end function neoclassics_calc_D112