[Maple] Численно решить дифур с численно заданной функцией

[srm]

Здравствуйте. Нужно численно решить дифур, некоторые функции которого являются решением другого дифура. Привожу упрощённый пример (функции на самом деле другие):

sol1 := dsolve({diff(x(t), t) = x(t), x(0) = 0.0}, numeric);
func := q -> op(2, op(2, sol1(q)));
dsolve({diff(z(t), t) = func(t) * z(t), z(0) = 0.1}, numeric);

Понятное дело, этот код не работает, ибо пытается вызвать op(2, sol1(t))...
Как быть?

[Andrew Gubarev]

Решать систему

>sys:= {diff(x(t), t) = x(t), diff(z(t), t) = x(t) * z(t), x(0)=0.1, z(0)=0.1}:
> f:= dsolve (sys, numeric, range=0..3):
> plots[odeplot](f, [t, z(t)]); # вывод на экран решения

(Изменил одно из начальных условий, поскольку оно было для меня не очень показательным.)

[srm]

Andrew Gubarev, нет, такой вариант не подходит. Я уже писал, что реальный пример существенно сложней и там такое сделать не получится. ответ на свой вопрос я уже получил:

нужно задать дополнительный параметр в dsolve:

output = listprocedure

потом задавать функцию-решение ДУ как

f := eval(x(t), sol1)

.

-- Пн апр 23, 2012 18:14:11 --

Если кому-то понадобится, подробный ответ здесь http://www.mapleprimes.com/questions/133218-How-Can-I-Solve-DE-With-Function-Which?submit=133303#comment133303