Функцию не просчитывает, что не так?

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

Модератор: Admin

Полищук Андрей
Сообщения: 55
Зарегистрирован: Чт ноя 09, 2006 2:40 pm
Откуда: Ижевск

Функцию не просчитывает, что не так?

Сообщение Полищук Андрей » Чт апр 12, 2007 12:04 pm

Следующая ситуация:
Получил характеристич. полином, это неважно - просто функция - HarUr.
Ввожу уравнение с целью далее работать с ним.
HarUr3:= HarUr =0;
Делаю алгебраическую подстановку, чтобы возможно было получить производную (т.е. вместо omega исп. omega(k))
HarUr3:= collect(algsubs(omega2/(2*i) = omega3(k),expand(HarUr3)), k);
Далее преобразования согласованы для HarUr и HarUr3, плюс беру производную:
diff(HarUr3,k);
изолирую производную
GrVel:=rhs(isolate(%,diff(omega3(k),k))):
подстановка omega3(k)-> omega3
GrVel:= algsubs(omega3(k) = omega3, GrVel);
Это групповая скорость, а HarUr -- дисперсионное (частотное) уравнение. Его решаю относительно частот, подставдяю в GrVel --- НЕ РАБОТАЕТ. Причем GrVel получена в аналитическом виде.
Подставляю - это фрагмент
omega3:= Otv[NN][3];
это результат для конкретного числа
omega3 := .9854059595
Получаю серию ответов, в которых нет ни одной не переменной с неприсвоенным значением и получаю после evalf(GrVel) нечто вроде:
.255016263810^7sin(alpha0)(1. + m)... и alpha0 и m присвоены конкретные значения.

В чем причина?

Irmantas
Сообщения: 589
Зарегистрирован: Ср апр 27, 2005 8:49 pm
Откуда: Литва, Каунас

Re: Функцию не просчитывает, что не так?

Сообщение Irmantas » Чт апр 12, 2007 12:19 pm

Ничего неясно.
Пришлите сам мапле файл