Дискуссионный математический форумМатематический форум
Математический форум Math Help Planet

Обсуждение и решение задач по математике, физике, химии, экономике

Теоретический раздел
Часовой пояс: UTC + 3 часа [ Летнее время ]
новый онлайн-сервис
число, сумма и дата прописью

Часовой пояс: UTC + 3 часа [ Летнее время ]




Начать новую тему Ответить на тему  [ 1 сообщение ] 
Автор Сообщение
 Заголовок сообщения: Решение сиcтемы уравнений методом секущих2
СообщениеДобавлено: 02 июн 2021, 20:25 
Не в сети
Продвинутый
Зарегистрирован:
24 дек 2019, 00:16
Сообщений: 58
Cпасибо сказано: 18
Спасибо получено:
0 раз в 0 сообщении
Очков репутации: 1

Добавить очки репутацииУменьшить очки репутации
Доброго времени суток) Нужна помощь в понимании готовой процедуры решения системы нелинейных уравнений методом секущих:
M2 := unapply( (Mnew(X1,X2,X3)), X1,X2,X3);
Q2 := unapply( (Qnew(X1,X2,X3)), X1,X2,X3);

sol1 := proc(x11,x13,x2)
global X30;
local m1,q1,m2,q2,m3,q3, X11,X12,X13,X2,i;
# print(` sol1:start `);
X11, X13, X2 := x11, x13, x2;
m1 := evalf(M2(X11,X2,X30));
q1 := evalf(Q2(X11,X2,X30));
# print(X11, X2, m1,q1);
m3 := evalf(M2(X13,X2,X30));
q3 := evalf(Q2(X13,X2,X30));
# print(X13, X2, m3,q3);
m2 := m1;
i := 1;
while abs(m2) > 1e-6 and i < 20 do
X12 := -(X11*m3 - X13*m1)/(m1-m3);
m2 := evalf(M2(X12,X2,X30));
q2 := evalf(Q2(X12,X2,X30));
# print(X12, X2, m2,q2);
X13 := X12;
m3 := m2;
i := i + 1;
end do;
# print(` sol1:end `);
return(X12, X2, m2, q2);
end proc:
sol2 := proc(x11,x13,x21,x23)
global X30;
local m1,q1,m2,q2,m3,q3, X11,X12,X13,X21,X22,X23,i, dX1;
X11, X13, X21, X23 := x11, x13, x21, x23;
dX1 := (x13-x11);
X12, X21, m1, q1 := sol1(X11, X13, X21);
# print(X12, X21, m1, q1);
X12, X23, m3, q3 := sol1(X11, X13, X23);
# print(X12, X23, m3, q3);
q2 := q1;
i := 1;
while abs(q2) > 1e-6 and i < 20 do
X22 := -(X21*q3 - X23*q1)/(q1-q3);
X12, X22, m2, q2 := sol1(X11, X13, X22);
# print(X12, X22, m2, q2, X11,X13,X21,X23);
X11 := X12 - dX1/2;
X13 := X12 + dX1/2;
X23 := X22;
q3 := q2;
i := i + 1;
end do;
return(X12, X22, m2, q2);
end proc:

Вызываю в виде, например:
X1 := .72;
X2 := .952;
X30 := -2;
X1,X2,mm,qq := sol2(X1,X1*1.2,X2,X2*1.1);

опишите, пожалуйста, работу sol1 и sol2

Вернуться к началу
 Профиль  
Cпасибо сказано 
Показать сообщения за:  Поле сортировки  
Начать новую тему Ответить на тему      Страница 1 из 1 [ 1 сообщение ]

 Похожие темы   Автор   Ответы   Просмотры   Последнее сообщение 
Решение системы уравнений методом Гаусса

в форуме Линейная и Абстрактная алгебра

Alinmora

3

408

13 дек 2015, 19:53

Решение линейных уравнений методом Гаусса (4х5)

в форуме Линейная и Абстрактная алгебра

Alien_Ted

0

308

13 янв 2016, 18:56

Решение системы уравнений методом Гаусса

в форуме Линейная и Абстрактная алгебра

Kostya151

3

408

04 ноя 2015, 11:38

Решение системы уравнений методом сложения, умножения и тд

в форуме Линейная и Абстрактная алгебра

Oleg9

1

284

13 окт 2015, 14:39

Решение системы нелинейных уравнений методом ньютона

в форуме Линейная и Абстрактная алгебра

mathematic

1

402

15 фев 2018, 12:19

Система линейных уравнений, решение методом обратной м-цы

в форуме Линейная и Абстрактная алгебра

milageras

10

1086

07 окт 2015, 19:09

Решение дифференциальных уравнений н порядка методом Эйлера

в форуме Численные методы

jonygibson

10

890

02 фев 2015, 10:17

Решение системы нелинейных уравнений геометрическим методом

в форуме Интересные задачи участников форума MHP

AGN

6

174

09 авг 2023, 19:08

Mathcad 14 - Решение систем нелинейных уравнений методом Нью

в форуме Численные методы

Ilya_Sergeev

0

259

07 июл 2021, 17:22

Найти общее решение системы линейных уравнений методом Гаусс

в форуме Линейная и Абстрактная алгебра

Elena116

11

1142

20 окт 2014, 12:44


Часовой пояс: UTC + 3 часа [ Летнее время ]



Кто сейчас на конференции

Сейчас этот форум просматривают: нет зарегистрированных пользователей и гости: 7


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

Найти:
Перейти:  

Яндекс.Метрика

Copyright © 2010-2023 MathHelpPlanet.com. All rights reserved