Thank you very much
Your great hints are useful for me
But I have more question (not only about lsode)
I will be grateful, if you can help me again
Please download the attached file
The value of C is assumed to be zero for relative error checking ( as you recommended)
1- It is reasonable that by assigning a little value to C, the little tolerance of C for little amount of absolute maximum step size in interval [0,1], should be ineffective to the answer. By taking the value C=Float(1,-15) and C=1E-3, the answer is 0.229491 but for C=Float(1.-14) the answer is 0.221455. Is it show the instability of solution?
2- I use the Implode and Explode commands in line 14. Is there a better way for obtaining the expression from 'RootOf' without using the String Tools?
3- Can I solve the ode in line 12 directly by maple commands?