07.10.2012, 17:18
|
#1 (permalink) |
|
Member
Регистрация: 27.11.2011
Сообщений: 29
Сказал(а) спасибо: 0
Поблагодарили 0 раз(а) в 0 сообщениях
Репутация: 10
|
Решить систему уравнений методом Ньютона, Pascal
или дайте хотя бы пример такой же задачи с комментариями, я сам переделаю. |
|
|
|
07.10.2012, 17:18
|
|
|
Helpmaster
Member
Регистрация: 08.03.2016
Сообщений: 0
|
Возможно, что вы еще не читали аналогичные по содержанию топики Pascal, Delphi. Шифрование методом Льюиса Помогите решить задачу в Pascal ABC Помогите с заданием: решить систему линейных алгебраических уравнений Решить задачку по Pascal |
|
07.10.2012, 19:28
|
#2 (permalink) | |
|
Специалист
Регистрация: 27.08.2008
Адрес: Санкт-Петербург
Сообщений: 27,807
Сказал(а) спасибо: 340
Поблагодарили 583 раз(а) в 208 сообщениях
Репутация: 113184
|
Цитата:
Теперь сама программа: Код:
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.
|
|
|
|
|
|
08.10.2012, 20:32
|
#5 (permalink) |
|
Member
Регистрация: 27.11.2011
Сообщений: 29
Сказал(а) спасибо: 0
Поблагодарили 0 раз(а) в 0 сообщениях
Репутация: 10
|
я не понимаю как в маткаде график этот построить...
подскажите |
|
|
|
| Ads | |
|
Member
Регистрация: 31.10.2006
Сообщений: 40200
Записей в дневнике: 0
Сказал(а) спасибо: 0
Поблагодарили 0 раз(а) в 0 сообщениях
Репутация: 55070
|
|
|
08.10.2012, 21:06
|
#6 (permalink) |
|
Member
Регистрация: 27.11.2011
Сообщений: 29
Сказал(а) спасибо: 0
Поблагодарили 0 раз(а) в 0 сообщениях
Репутация: 10
|
вот что получилось, я не могу проверить правильность, т.к. не могу построить маткад
помогите пожалуйста, завтра сдавать Код:
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.
|
|
|
|
|
09.10.2012, 11:57
|
#8 (permalink) |
|
Специалист
Регистрация: 27.08.2008
Адрес: Санкт-Петербург
Сообщений: 27,807
Сказал(а) спасибо: 340
Поблагодарили 583 раз(а) в 208 сообщениях
Репутация: 113184
|
Да... посмотрел (извините, вчера не было возможности ответить). И зачем понадобилось функции переделывать? Получилась чушь несусветная и работать это не будет.
Прямо руки опускаются помогать таким умникам. Тогда уж пишите программы сами! |
|
|
|
|
09.10.2012, 12:11
|
#9 (permalink) |
|
Member
Регистрация: 27.11.2011
Сообщений: 29
Сказал(а) спасибо: 0
Поблагодарили 0 раз(а) в 0 сообщениях
Репутация: 10
|
ну я мог ничего совсем ненаписать, я буквально вчера познакомился с Маткадом.
в общем сделал я эту задачу, используя частные производные. ждите следующих)) |
|
|
|
| Ads | |
|
Member
Регистрация: 31.10.2006
Сообщений: 40200
Записей в дневнике: 0
Сказал(а) спасибо: 0
Поблагодарили 0 раз(а) в 0 сообщениях
Репутация: 55070
|
|
 |
|
|
Линейный вид
