Сообщение aar » Пн май 26, 2008 9:24 am
Не знаю, насколько это законно. Проверяйте:
> S1 := dsolve({ur1, ur2}) assuming t>=0:
> S2 := dsolve({ur1, ur2}) assuming t<=0:
> S2 := eval(S2, {_C1 = _C3, _C2 = _C4}):
> eval(S1, S2):
> S:= eval(%, t = 0):
Дальше надо смотреть, что получилось в S. Поскольку у меня вышло:
> S[2];
_C4 = -_C1 + _C2
то можно сделать так:
> eval(S[1], S[2]):
> solve(%);
Получается ответ из четырех вариантов, в одном варианте есть lambda = A, в другом lambda = -A.