static void Main() { newtonRaphson SL = new newtonRaphson(0.0001, 0.01, 2, 2); GR = new GoldenRatio(); for (SL.k = 0; ; SL.k++) { Console.WriteLine(); Console.WriteLine("--------------------------------------"); Console.WriteLine("\nточка:" + SL.next.x + " " + SL.next.y); Console.WriteLine("итерация: " + SL.k); Console.WriteLine("градиент: " + SL.grad.x+" "+SL.grad.y); Console.WriteLine("норма градиента: " + SL.normV(SL.grad)); if (SL.normV(SL.grad) <= SL.e1) { SL.answer = SL.next; break; } if ((SL.m.Invert().Det()) > 0) { SL.d = -((SL.m.Invert()) * SL.grad); } else { SL.d.x = -SL.grad.x; SL.d.y = -SL.grad.y; } SL.step = GR.Calculate(SL.next, SL.d); Console.WriteLine("шаг: " + SL.step); //if(k>1000) { break;} SL.next.x = SL.first.x + SL.step * SL.d.x; SL.next.y = SL.first.y + SL.step * SL.d.y; // Console.WriteLine(k); if (SL.normV(new Point(SL.next.x - SL.first.x, SL.next.y - SL.first.y)) < SL.e2) { if (Math.Abs(SL.F(SL.next.x, SL.next.y) - SL.F(SL.first.x, SL.first.y)) < SL.e2) { SL.answer = SL.next; break; } } SL.first.x = SL.next.x; SL.first.y = SL.next.y; SL.grad = SL.gradF(SL.next); } SL.answer.x = SL.next.x; SL.answer.y = SL.next.y; Console.WriteLine(); Console.WriteLine("--------------------------------------"); Console.WriteLine("--------------------------------------"); Console.WriteLine("\nОтвет:" + SL.answer.x + " " + SL.answer.y); Console.WriteLine("итерация: " + SL.k); Console.WriteLine("норма градиента: " +SL.normV(SL.grad)); Console.Read(); }
newtonRaphson(double e1, double e2, double x, double y) { this.e1 = e1; this.e2 = e2; first = new Point(x, y); next = new Point(first.x, first.y); grad = gradF(next); answer = new Point(); d = new Point(-grad.x, -grad.y); preGrad = new Point(grad.x, grad.y); preD = new Point(d.x, d.y); GR = new GoldenRatio(); m = new Matrix(16, -1, -1, 2); }
static void Main() { newtonRaphson SL = new newtonRaphson(0.0001, 0.01, 2, 2); GR = new GoldenRatio(); for (SL.k = 0; ; SL.k++) { Console.WriteLine(); Console.WriteLine("--------------------------------------"); Console.WriteLine("\nточка:" + SL.next.x + " " + SL.next.y); Console.WriteLine("итерация: " + SL.k); Console.WriteLine("градиент: " + SL.grad.x + " " + SL.grad.y); Console.WriteLine("норма градиента: " + SL.normV(SL.grad)); if (SL.normV(SL.grad) <= SL.e1) { SL.answer = SL.next; break; } if ((SL.m.Invert().Det()) > 0) { SL.d = -((SL.m.Invert()) * SL.grad); } else { SL.d.x = -SL.grad.x; SL.d.y = -SL.grad.y; } SL.step = GR.Calculate(SL.next, SL.d); Console.WriteLine("шаг: " + SL.step); //if(k>1000) { break;} SL.next.x = SL.first.x + SL.step * SL.d.x; SL.next.y = SL.first.y + SL.step * SL.d.y; // Console.WriteLine(k); if (SL.normV(new Point(SL.next.x - SL.first.x, SL.next.y - SL.first.y)) < SL.e2) { if (Math.Abs(SL.F(SL.next.x, SL.next.y) - SL.F(SL.first.x, SL.first.y)) < SL.e2) { SL.answer = SL.next; break; } } SL.first.x = SL.next.x; SL.first.y = SL.next.y; SL.grad = SL.gradF(SL.next); } SL.answer.x = SL.next.x; SL.answer.y = SL.next.y; Console.WriteLine(); Console.WriteLine("--------------------------------------"); Console.WriteLine("--------------------------------------"); Console.WriteLine("\nОтвет:" + SL.answer.x + " " + SL.answer.y); Console.WriteLine("итерация: " + SL.k); Console.WriteLine("норма градиента: " + SL.normV(SL.grad)); Console.Read(); }