/*
 * 示例程序Sample5_9: NLEquations類的求非線性方程組一組實根的擬牛頓法
 */

using System;
using CSharpAlgorithm.Algorithm;

namespace CSharpAlgorithm.Sample
{
	class Class1
	{
        // 建立NLEquation的子類，在其中重載函數func
        class NLEq : NLEquations
        {
            protected override double Func(double[] x, double[] y)
            {
                y[0]=x[0]*x[0]+x[1]*x[1]+x[2]*x[2]-1.0;
                y[1]=2.0*x[0]*x[0]+x[1]*x[1]-4.0*x[2];
                y[2]=3.0*x[0]*x[0]-4.0*x[1]+x[2]*x[2];

                return 0.0;
            }
        }
        
        [STAThread]
		static void Main(string[] args)
		{
            // 求解
            NLEq nleq = new NLEq();
            double[] x={1.0,1.0,1.0};
            double t=0.1; 
            double h=0.1;
            double eps = 0.0000001;
            if (nleq.GetRootsetNewton(3, x, t, h, 100, eps))
            {
                for (int i=0; i<3; ++i)
                {
                    Console.WriteLine(x[i]);
                }
            }
            else
                Console.WriteLine("求解失敗");
        }
	}
}