Dear All,
I have the below equation:
sol = NDSolve[{q'[t] == (a q[t]^b newefr[t] - q[t]) + (a q[t]^b (newefr[t] - newefr1[t])), q[1] == 0.001}, q, {t, 1, 456}]
Then I want to extract the value of q for each time step by using below command:
Evaluate[q[t] /. sol]
For example for t=430:
Evaluate[q[10] /. sol] result: {-8.60556*10^-13 - 1.18404*10^-16 i}
So may I know why does the solution contain imaginary number? Thanks in advance!
Best Regards, Intan