class Program
    {
        static void Main(string[] args)
        {
            double[] y = new double[] { 2, 2, 3, 3};
         
            double[,] x = new double[y.Length, 1];
           
            for (int i = 0; i < y.Length; i++)
            {
                x[i, 0] = (double)(i);
            }

            double[] s = new double[] {1, 1};       // scale
            double[] c = new double[] {2.2, 1};    // start values
            double[] bndl = new double[] {0.0, 0.0};     // lower bound                                                                                                                                                     
            double[] bndu = new double[] {4.0, 4.0};     // upper bound                                                                                                                                                             
            double epsx = 0.0001;
            int maxits = 0;
            int info;
            alglib.lsfitstate state;
            alglib.lsfitreport rep;
            double diffstep = 0.01;

            alglib.lsfitcreatef(x, y, c, diffstep, out state);         
            alglib.lsfitsetbc(state, bndl, bndu);
            alglib.lsfitsetcond(state, epsx, maxits);
            alglib.lsfitsetscale(state, s);
            alglib.lsfitfit(state, function_cx_2_func, null, null);
            alglib.lsfitresults(state, out info, out c, out rep);
           
            int iterationscount = rep.iterationscount;

            System.Console.WriteLine("{0}", info);
       System.Console.WriteLine("iterationscount = " + iterationscount.ToString());
            System.Console.WriteLine("{0}", alglib.ap.format(c, 7));
            System.Console.ReadKey();
        }

        public static void function_cx_2_func(double[] c, double[] x, ref double func, object obj)
        {
            // vectors a and b
           
            [b]double[] a = new double[] { 1, 1, 0, 0};
            double[] b = new double[] { 0, 0, 1, 1};[/b]
           
            func = func + c[0] * a[(int)x[0]] + c[1] * b[(int)x[0]];
        }
    }

class Program
    {
        static void Main(string[] args)
        {
            //
            // This example demonstrates minimization of F(x0,x1) = f0^2+f1^2 + ....  fn^2       
            //
            // using "V" mode of the Levenberg-Marquardt optimizer.
            //
            // Optimization algorithm uses:
            // * function vector f[] = {f1,f2,...,fn^2}
            //
            // No other information (Jacobian, gradient, etc.) is needed.
            //

            double[] x = new double[] { 0, 0 };
            double epsx = 0.0000000001;
            int maxits = 100;                     
            alglib.minlmstate state;
            alglib.minlmreport rep;
            double[] bndl = new double[] { 0, 0 };               
            double[] bndu = new double[] { +4, +4 };   

            alglib.minlmcreatev(4, x, 0.0001, out state);
            alglib.minlmsetbc(state, bndl, bndu);           
            alglib.minlmsetcond(state, epsx, maxits);             
            alglib.minlmoptimize(state, function2_fvec, null, null);
            alglib.minlmresults(state, out x, out rep);

            System.Console.WriteLine("{0}", alglib.ap.format(x, 7));
            System.Console.WriteLine("{0}", alglib.ap.format(rep.iterationscount, 7));
            System.Console.ReadLine();
           
        }
       
        public static void function2_fvec(double[] x, double[] fi, object obj)
        {
            double[] array1 = new double[] { 1, 1, 0, 0 };
            double[] array2 = new double[] { 0, 0, 1, 1};
            double[] array3 = new double[] { 2, 2, 3, 3};       

            for (int i = 0; i < array1.Length; i++)
            {
                fi[i] = System.Math.Pow(array3[i] - x[0] * array1[i] - x[1] * array2[i], 2);
            }
        }
    }