численное решение системы ДУ

Форум пользователей пакета Mathcad

Модератор: Admin

aor1120
Сообщения: 30
Зарегистрирован: Вс сен 28, 2008 5:12 pm

численное решение системы ДУ

Сообщение aor1120 » Чт окт 08, 2009 5:30 pm

Всем привет, это опять я и мои ДУ.
Есть система из 4х ДУ, решаю ее численно. Если поставить вектор начальных условий, например, (1,0,1,0), а коэф. с5=1, то все замечательно считается. А при (1,1,1,1) и с5=-1, считать не хочет. Вопрос: почему? Каким интегратором лучше считать такую систему?

PS. При (1,1,1,1) и с5=0.1 на графике разностей I получается какая-то чушь. Что за?

Korobov V I
Сообщения: 1609
Зарегистрирован: Вт янв 21, 2003 5:12 pm

Сообщение Korobov V I » Чт окт 08, 2009 6:32 pm

Попробуйте AdamsBDF (MC14).

aor1120
Сообщения: 30
Зарегистрирован: Вс сен 28, 2008 5:12 pm

Сообщение aor1120 » Чт окт 08, 2009 9:49 pm

Попробовал с помощью AdamsBDF:
Все равно не считает при -0.9<c5<-0.2
Может я с [tol] чего-то не понял?
Последний раз редактировалось aor1120 Ср окт 14, 2009 7:50 pm, всего редактировалось 1 раз.

aor1120
Сообщения: 30
Зарегистрирован: Вс сен 28, 2008 5:12 pm

Сообщение aor1120 » Ср окт 14, 2009 7:45 pm

И все-таки, что ему надо (при c5<-0.2)? Что за ошибка такая "Found a number with a magnitude greater than 10^307 while trying to evaluate this expression"? Т.е. откуда он вообще взял такие значения.