//      d01aq Example Program Text
//      C# version, NAG Copyright 2008
using System;
using NagLibrary;
using System.Globalization;
namespace NagDotNetExamples
{
  public class D01AQE
  {
    public static int kount = 0;
    static void Main(String[] args)
    {
      StartExample();
    }
    public static void StartExample()
    {
      try
      {
        const int lw = 800;
        D01.D01AQ_G gD01AQ = new D01.D01AQ_G(g);
        double a,  abserr,  b,  c,  epsabs,  epsrel,  result;
        double[] w = new double[lw];
        int subintvls;
        int ifail;
        Console.WriteLine("d01aq Example Program Results");
        epsabs = 0.00e0;
        epsrel = 1.00e-04;
        a = -1.00e0;
        b = 1.00e0;
        c = 0.50e0;
        kount = 0;
        // 
        D01.d01aq(gD01AQ, a, b, c, epsabs, epsrel, out result, out abserr, w, out subintvls,
                  out ifail);
        // 
        if (ifail >= 0)
        {
          Console.WriteLine("");
          Console.WriteLine("  {0}{1,10:f4}", "a      - lower limit of integration = ", a);
          Console.WriteLine("  {0}{1,10:f4}", "b     - upper limit of integration = ", b);
          Console.WriteLine("  {0}{1,9:e2}", "epsabs - absolute accuracy requested = ", epsabs);
          Console.WriteLine("  {0}{1,9:e2}", "epsrel - relative accuracy requested = ", epsrel);
          Console.WriteLine("  {0}{1,9:e2}", "C      - parameter in the weight function = ", c);
          Console.WriteLine("");
          if (ifail != 0)
          {
            Console.WriteLine("  {0}{1,4}", "ifail = ", ifail);
          }
          if (ifail <= 3)
          {
            Console.WriteLine("  {0}{1,9:f2}", "result - approximation to the integral = ", result);
            Console.WriteLine("  {0}{1,9:e2}", "abserr - estimate of the absolute error = ", abserr);
            Console.WriteLine("  {0}{1,4}", "kount  - number of function evaluations = ", kount);
            Console.WriteLine("  {0}{1,4}", "subintvls  - number of subintervals used = ", subintvls);
          }
        }
        else
        {
          Console.WriteLine("");
          Console.WriteLine("** d01aq failed with ifail = {0,5}", ifail);
        }
        // 
      }
      catch (Exception e)
      {
        Console.WriteLine(e.Message);
        Console.WriteLine("Exception Raised");
      }
      ;
    }
    // 
    public static double g(double x)
    {
      double gValue = 0.0;
      double aa = 0.0;
      kount = kount + 1;
      aa = 0.010e0;
      gValue = 1.00e0 / ((x) * (x) + (aa) * (aa));
      return gValue;
    }
  }
}