Здравствуйте. Нужно численно решить дифур, некоторые функции которого являются решением другого дифура. Привожу упрощённый пример (функции на самом деле другие):
Код:
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))...
Как быть?