コード例 #1
0
    public static double p04_q(ref typeMethods.r8ErrorData data, ref typeMethods.r8ErrorcData cdata, int m, double[] c, double[] w)
//****************************************************************************80
//
//  Purpose:
//
//    P04_Q evaluates the integral for problem p04.
//
//  Licensing:
//
//    This code is distributed under the GNU LGPL license.
//
//  Modified:
//
//    29 August 2012
//
//  Author:
//
//    John Burkardt
//
//  Reference:
//
//    Alan Genz,
//    A Package for Testing Multiple Integration Subroutines,
//    in Numerical Integration: Recent Developments, Software
//    and Applications,
//    edited by Patrick Keast and Graeme Fairweather,
//    Reidel, 1987, pages 337-340,
//    ISBN: 9027725144,
//    LC: QA299.3.N38.
//
//  Parameters:
//
//    Input, int M, the dimension of the argument.
//
//    Input, double C[M], W[M], the problem parameters.
//
//    Output, double P04_Q, the integral.
//
    {
        double q = 1.0;

        for (int i = 0; i < m; i++)
        {
            q = q * Math.Sqrt(Math.PI)
                * (typeMethods.r8_error(ref data, ref cdata, c[i] * (1.0 - w[i]))
                   + typeMethods.r8_error(ref data, ref cdata, c[i] * w[i]))
                / (2.0 * c[i]);
        }

        return(q);
    }
コード例 #2
0
    public static double p00_q(ref typeMethods.r8ErrorData data, ref typeMethods.r8ErrorcData cdata, int prob, int m, double[] c, double[] w)
//****************************************************************************80
//
//  Purpose:
//
//    P00_Q returns the integral of any function.
//
//  Licensing:
//
//    This code is distributed under the GNU LGPL license.
//
//  Modified:
//
//    29 August 2012
//
//  Author:
//
//    John Burkardt
//
//  Parameters:
//
//    Input, int PROB, the index of the function.
//
//    Input, int M, the spatial dimension.
//
//    Input, double C[M], W[M], the problem parameters.
//
//    Output, double P00_Q, the integral.
//
    {
        double q;

        switch (prob)
        {
        case 1:
            q = p01_q(m, c, w);
            break;

        case 2:
            q = p02_q(m, c, w);
            break;

        case 3:
            q = p03_q(m, c, w);
            break;

        case 4:
            q = p04_q(ref data, ref cdata, m, c, w);
            break;

        case 5:
            q = p05_q(m, c, w);
            break;

        case 6:
            q = p06_q(m, c, w);
            break;

        default:
            Console.WriteLine("");
            Console.WriteLine("P00_Q - Fatal error!");
            Console.WriteLine("  Illegal function index PROB = " + prob + "");
            return(1);
        }

        return(q);
    }