Здравствуйте!
Создал новую тему (хотя похожие есть), т.к. к имеющимся уж очень давно в последний раз обращались! Решение о правильности моего поступка, по всей видимости, остается за модераторами. Очень нужна помощь в этом вопросе, поскольку обратиться больше не к кому. Вопрос в общем-то не сложный, но странный
Есть система ДУ с условиями:
Код:
sys:={14.52*a[3](t)+18.60*a[1](t)+2.35*(diff(a[3](t), t)), 18.60*a[3](t)+6.28*(diff(a[1](t), t))+33.61*a[1](t)};
ics:={a[1](0) = .63, a[3](1) = .21}
Решая ее в Maple 13 и 15 получаю
Код:
{a[1](t) = Float(undefined)*exp(-.9051675582*t)+Float(undefined)*exp(-10.60750253*t), a[3](t)=Float(undefined)*exp(-.9051675582*t)+Float(undefined)*exp(-10.60750253*t)}
где
Float(undefined) -- насколько я понимаю, не определенное вещественное число. Возникает такое, вроде как, когда пытаешься найти значение функции, которое не принадлежит ее области определения. Однако, решая систему без условий получаем
Код:
dsold := evalf(dsolve(sys);assign(dsold):
{a[1](t) = _C1*exp(-.9051675582*t)+_C2*exp(-10.60750253*t), a[3](t) = -1.501489565*_C1*exp(-.9051675582*t)+1.776014119*_C2*exp(-10.60750253*t)}
далее подставляем
соответственно в
Код:
raz := subs(t = 0, a[1](t)):dva := subs(t = 1, a[3](t)):
solve({dva = .21, raz = .63}, {_C1, _C2});
{_C1 = -.35, _C2 = .99}
находим
. Как это объяснить? И почему Maple не может решить эту систему с условиями сразу, найдя эти константы?
Задача у меня именно такая: одни условия задаются в начале временного промежутка, другие -- в конце. Система решается "в лоб" мейплом, если задавать
, но только в этом случае. Если оба условия задавать в
-- тоже не решается. С чем это связано и почему так происходит? Может я все-таки что-то делаю не правильно?