//      g13bd Example Program Text
//      C# version, NAG Copyright 2008
using System;
using NagLibrary;
using System.Globalization;
namespace NagDotNetExamples
{
  public class G13BDE
  {
    static string datafile = "ExampleData/g13bde.d";
    static void Main(String[] args)
    {
      if (args.Length == 1)
      {
        datafile = args[0];
      }
      StartExample();
    }
    public static void StartExample()
    {
      try
      {
        DataReader sr = new DataReader(datafile);
        double r0,   s; int i,  iwa,  nl,  nwds;
        int[] isf = new int[2];
        int[] nna = new int[3];
        int ifail;
        Console.WriteLine("g13bd Example Program Results");
        //      Skip heading in data file
        sr.Reset();
        sr.Reset();
        nl = int.Parse(sr.Next());
        sr.Reset();
        r0 = double.Parse(sr.Next(), CultureInfo.InvariantCulture);
        double[] r = new double[nl];
        if (nl > 0)
        {
          sr.Reset();
          for (i = 1; i <= nl; i++)
          {
            r[i - 1] = double.Parse(sr.Next(), CultureInfo.InvariantCulture);
          }
          sr.Reset();
          for (i = 1; i <= 3; i++)
          {
            nna[i - 1] = int.Parse(sr.Next());
          }
          sr.Reset();
          s = double.Parse(sr.Next(), CultureInfo.InvariantCulture);
          nwds = nna[1] + nna[2] + 1;
          iwa = nna[2] * (nna[2] + 1);
          double[] wds = new double[nwds];
          if (nwds > 0 )
          {
            //
            G13.g13bd(r0, r, nl, nna, s, nwds, wds, isf, out ifail);
            //
            Console.WriteLine("");
            Console.WriteLine("  {0}{1,4}{2,4}", "Success/failure indicator", isf[0], isf[1]);
            Console.WriteLine("");
            Console.Write(" {0}", "Transfer function model B, Q, P =");
            for (i = 1; i <= 3; i++)
            {
              Console.Write(" {0}", nna[i - 1]);
            }
            Console.WriteLine(""); ;
            Console.WriteLine("");
            Console.WriteLine(" {0}", "Parameter initial estimates\r\n");
            for (i = 1; i <= nwds; i++)
            {
              Console.Write(" {0:0.0000}", wds[i - 1]);
            }
            Console.WriteLine(""); ;
          }
        }
        //
      }
      catch (Exception e)
      {
        Console.WriteLine(e.Message);
        Console.WriteLine("Exception Raised");
      }
    }
  }
}