// d01ak Example Program Text // C# version, NAG Copyright 2008 using System; using NagLibrary; using System.Globalization; namespace NagDotNetExamples { public class D01AKE { public static int kount = 0; static void Main(String[] args) { StartExample(); } public static void StartExample() { try { const int lw = 800; D01.D01AK_F fD01AK = new D01.D01AK_F(f); double a, abserr, b, epsabs, epsrel, pi, result; double[] w = new double[lw]; int subintvls; int ifail; Console.WriteLine("d01ak Example Program Results"); pi = Math.PI; epsabs = 0.00e0; epsrel = 1.00e-03; a = 0.00e0; b = 2.00e0 * pi; kount = 0; // D01.d01ak(fD01AK, 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 <= 3) { 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("** d01ak 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.Cos(x); return fValue; } } }