This is my problem:
i've an nonlinear system of six equations and i find the
solution with Runge-Kutta:
for every t>0 i've
x(t), x(t), x(t), xc(t), xc(t), xc(t).
I write this solution in the vector "solution"
Yet i've also another equation ("EQ7") of the type:
EQ7 = f(xc(t),xc(t),xc) with xc,xc,xc presents in "solution".
I desire the graph of
EQ7 = f(t)
I've write this:
> for qtime from 0 to time by 150 do
axes = normal,
title = `without saturation`,
but this "for" cycle in many onerous. Have you another solution?