Технический форум - Решить систему уравнений методом Ньютона, Pascal

Технический форум (http://www.tehnari.ru/)

- Помощь студентам (http://www.tehnari.ru/f41/)

Решить систему уравнений методом Ньютона, Pascal

Решить систему методом Ньютона с точностью e= 0,00001.
или дайте хотя бы пример такой же задачи с комментариями, я сам переделаю.

Цитата:

Сообщение от ayPinki (Сообщение 801805)

Вложение 99410 Вложение 99411
Теперь сама программа:

Код:

Const

 e=0.00001;



Var

 x_old, x_new, y, z:real;



Function Fun_1(x:real):real;

begin

 Fun_1:=Sin(x-Sqrt(1-Sqr(x)))-1.5*x;

end;



Function Fun_2(x:real):real;

begin

 Fun_2:=Sin(x+Sqrt(1-Sqr(x)))-1.5*x;

end;



Function Der_1(x:real):real;

begin

 Der_1:=Cos(x-Sqrt(1-Sqr(x)))*(1+x/Sqrt(1-Sqr(x)))-1.5;

end;



Function Der_2(x:real):real;

begin

 Der_2:=Cos(x+Sqrt(1-Sqr(x)))*(1-x/Sqrt(1-Sqr(x)))-1.5;

end;



Begin    

 x_old:=-0.6;

 Repeat

  z:=x_old;

  x_new:=x_old-Fun_1(x_old)/Der_1(x_old);

  x_old:=x_new;

 Until Abs(x_old-z)<e;

 y:=-Sqrt(1-Sqr(x_old));

 Writeln('x1= ',x_old:0:5);

 Writeln('y1= ',y:0:5);

 Writeln('Sin(x1+y1)= ',Sin(x_old+y):0:5);

 Writeln('1.5*x1= ',1.5*x_old:0:5);

 Writeln;

 x_old:=0.6;

 Repeat

  z:=x_old;

  x_new:=x_old-Fun_2(x_old)/Der_2(x_old);

  x_old:=x_new;

 Until Abs(x_old-z)<e;

 y:=Sqrt(1-Sqr(x_old));

 Writeln('x2= ',x_old:0:5);

 Writeln('y2= ',y:0:5);

 Writeln('Sin(x2+y2)= ',Sin(x_old+y):0:5);

 Writeln('1.5*x2= ',1.5*x_old:0:5);

 Readln;

End.

И результат:
Вложение 99412

сразу моя программа? спс. будет разбираться.

Цитата:

Сообщение от ayPinki (Сообщение 801936)

спс. будет разбираться.

Успехов!_____

я не понимаю как в маткаде график этот построить...
подскажите

вот что получилось, я не могу проверить правильность, т.к. не могу построить маткад
помогите пожалуйста, завтра сдавать

Код:

Const

 e=0.00001;



Var

 x_old, x_new, y, z:real;



Function Fun_1(x:real):real;

begin

 Fun_1:=sin(x+y)-1.5*x;

end;



Function Fun_2(x:real):real;

begin

 Fun_2:=x*x+y*y-1;

end;



Function Der_1(x:real):real;

begin

 Der_1:=cos(x+y)-1.5;

end;



Function Der_2(x:real):real;

begin

 Der_2:=2*x+2*y;

end;



Begin

 x_old:=-0.6;

 Repeat

  z:=x_old;

  x_new:=x_old-Fun_1(x_old)/Der_1(x_old);

  x_old:=x_new;

 Until Abs(x_old-z)<e;

 y:=-Sqrt(1-Sqr(x_old));

 Writeln('x1= ',x_old:0:5);

 Writeln('y1= ',y:0:5);

 Writeln('Sin(x1+y1)= ',Sin(x_old+y):0:5);

 Writeln('1.5*x1= ',1.5*x_old:0:5);

 Writeln;

 x_old:=0.6;

 Repeat

  z:=x_old;

  x_new:=x_old-Fun_2(x_old)/Der_2(x_old);

  x_old:=x_new;

 Until Abs(x_old-z)<e;

 y:=Sqrt(1-Sqr(x_old));

 Writeln('x2= ',x_old:0:5);

 Writeln('y2= ',y:0:5);

 Writeln('Sin(x2+y2)= ',Sin(x_old+y):0:5);

 Writeln('1.5*x2= ',1.5*x_old:0:5);

 Readln;

End.

не получается "внести" вторую функцию в график , help ((
p/s/ и да, там ... -1...

Да... посмотрел (извините, вчера не было возможности ответить). И зачем понадобилось функции переделывать? Получилась чушь несусветная и работать это не будет.
Прямо руки опускаются помогать таким умникам. Тогда уж пишите программы сами!

ну я мог ничего совсем ненаписать, я буквально вчера познакомился с Маткадом.
в общем сделал я эту задачу, используя частные производные. ждите следующих))