//      d01aj Example Program Text
//      C# version, NAG Copyright 2008
using System;
using NagLibrary;
using System.Globalization;
namespace NagDotNetExamples
{
  public class D01AJE
  {
    public static double pi;
    public static int kount = 0;
    static void Main(String[] args)
    {
      StartExample();
    }
    public static void StartExample()
    {
      try
      {
        const int lw = 800;
        D01.D01AJ_F fD01AJ = new D01.D01AJ_F(f);
        double a,   abserr,   b,   epsabs,   epsrel,   result;
        double[] w = new double[lw]; 
        int subintvls;
        int ifail;
        Console.WriteLine("d01aj Example Program Results");
        pi = Math.PI;
        epsabs = 0.00e0;
        epsrel = 1.00e-04;
        a = 0.00e0;
        b = 2.00e0 * pi;
        kount = 0;
        // 
        D01.d01aj(fD01AJ, a, b, 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("");
          if (ifail != 0)
          {
            Console.WriteLine("  {0}{1,4}", "ifail = ", ifail);
          }
          if (ifail <= 5)
          {
            Console.WriteLine("  {0}{1,9:f5}", "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("** d01aj failed with ifail = {0,5}", ifail);
        }
        // 
      }
      catch (Exception e)
      {
        Console.WriteLine(e.Message);
        Console.WriteLine("Exception Raised");
      }
      ;
    }
    //  
    public static double f(double x)
    {
      double fValue;
      kount = kount + 1;
      fValue = x * Math.Sin(30.00e0 * x) / Math.Sqrt(1.00e0 - (x) * (x) / (4.00e0 * (pi) * (pi)));
      return fValue;
    }
  }
}