Технический форум
Вернуться   Технический форум > Программирование > Форум программистов > Помощь студентам


Ответ
 
Опции темы Опции просмотра
Старый 07.10.2012, 17:18   #1 (permalink)
ayPinki
Member
 
Регистрация: 27.11.2011
Сообщений: 29
Сказал(а) спасибо: 0
Поблагодарили 0 раз(а) в 0 сообщениях
Репутация: 10
По умолчанию Решить систему уравнений методом Ньютона, Pascal

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

Старый 07.10.2012, 17:18
Helpmaster
Member
 
Аватар для Helpmaster
 
Регистрация: 08.03.2016
Сообщений: 0

Возможно, что вы еще не читали аналогичные по содержанию топики

Pascal, Delphi. Шифрование методом Льюиса
Помогите решить задачу в Pascal ABC
Помогите с заданием: решить систему линейных алгебраических уравнений
Решить задачку по Pascal

Старый 07.10.2012, 19:28   #2 (permalink)
Vladimir_S
Специалист
 
Аватар для Vladimir_S
 
Регистрация: 27.08.2008
Адрес: Санкт-Петербург
Сообщений: 26,360
Сказал(а) спасибо: 289
Поблагодарили 507 раз(а) в 165 сообщениях
Репутация: 91953
По умолчанию

Цитата:
Сообщение от ayPinki Посмотреть сообщение
Решить систему методом Ньютона с точностью e= 0,00001.
или дайте хотя бы пример такой же задачи с комментариями, я сам переделаю.
stm_1.jpg stm_2.jpg
Теперь сама программа:
Код:
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.
И результат:
stm_3.jpg
__________________
With Mozilla Firefox - straight to communism!
Vladimir_S вне форума   Ответить с цитированием
Старый 07.10.2012, 20:03   #3 (permalink)
ayPinki
Member
 
Регистрация: 27.11.2011
Сообщений: 29
Сказал(а) спасибо: 0
Поблагодарили 0 раз(а) в 0 сообщениях
Репутация: 10
По умолчанию

сразу моя программа? спс. будет разбираться.
ayPinki вне форума   Ответить с цитированием
Старый 07.10.2012, 20:05   #4 (permalink)
Vladimir_S
Специалист
 
Аватар для Vladimir_S
 
Регистрация: 27.08.2008
Адрес: Санкт-Петербург
Сообщений: 26,360
Сказал(а) спасибо: 289
Поблагодарили 507 раз(а) в 165 сообщениях
Репутация: 91953
По умолчанию

Цитата:
Сообщение от ayPinki Посмотреть сообщение
спс. будет разбираться.
Успехов!_____
__________________
With Mozilla Firefox - straight to communism!
Vladimir_S вне форума   Ответить с цитированием
Старый 08.10.2012, 20:32   #5 (permalink)
ayPinki
Member
 
Регистрация: 27.11.2011
Сообщений: 29
Сказал(а) спасибо: 0
Поблагодарили 0 раз(а) в 0 сообщениях
Репутация: 10
По умолчанию

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

Яндекс

Member
 
Регистрация: 31.10.2006
Сообщений: 40200
Записей в дневнике: 0
Сказал(а) спасибо: 0
Поблагодарили 0 раз(а) в 0 сообщениях
Репутация: 55070
Старый 08.10.2012, 21:06   #6 (permalink)
ayPinki
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.
ayPinki вне форума   Ответить с цитированием
Старый 08.10.2012, 21:22   #7 (permalink)
ayPinki
Member
 
Регистрация: 27.11.2011
Сообщений: 29
Сказал(а) спасибо: 0
Поблагодарили 0 раз(а) в 0 сообщениях
Репутация: 10
По умолчанию

не получается "внести" вторую функцию в график , help ((
p/s/ и да, там ... -1...
Миниатюры
777.png  
ayPinki вне форума   Ответить с цитированием
Старый 09.10.2012, 11:57   #8 (permalink)
Vladimir_S
Специалист
 
Аватар для Vladimir_S
 
Регистрация: 27.08.2008
Адрес: Санкт-Петербург
Сообщений: 26,360
Сказал(а) спасибо: 289
Поблагодарили 507 раз(а) в 165 сообщениях
Репутация: 91953
По умолчанию

Да... посмотрел (извините, вчера не было возможности ответить). И зачем понадобилось функции переделывать? Получилась чушь несусветная и работать это не будет.
Прямо руки опускаются помогать таким умникам. Тогда уж пишите программы сами!
__________________
With Mozilla Firefox - straight to communism!
Vladimir_S вне форума   Ответить с цитированием
Старый 09.10.2012, 12:11   #9 (permalink)
ayPinki
Member
 
Регистрация: 27.11.2011
Сообщений: 29
Сказал(а) спасибо: 0
Поблагодарили 0 раз(а) в 0 сообщениях
Репутация: 10
По умолчанию

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

Яндекс

Member
 
Регистрация: 31.10.2006
Сообщений: 40200
Записей в дневнике: 0
Сказал(а) спасибо: 0
Поблагодарили 0 раз(а) в 0 сообщениях
Репутация: 55070
Ответ

Опции темы
Опции просмотра

Ваши права в разделе
Вы не можете создавать новые темы
Вы не можете отвечать в темах
Вы не можете прикреплять вложения
Вы не можете редактировать свои сообщения

BB коды Вкл.
Смайлы Вкл.
[IMG] код Выкл.
HTML код Выкл.
Trackbacks are Вкл.
Pingbacks are Вкл.
Refbacks are Выкл.




Часовой пояс GMT +4, время: 05:27.

Powered by vBulletin® Version 6.2.5.
Copyright ©2000 - 2014, Jelsoft Enterprises Ltd.