Wow ! that's great. thanks very much. I'm an ameture user, therefore "CodeTools" is new to me, also I don't understand some parameters inside the procedure for EnergyDensity.
"epsilon" is machine epsilon, right ?
But what is _dO( or is it zero ?)1amc method ?
Aso when I run your worksheet , my version can't work with "Usage", it says : "Error, Usage is not a command in the CodeTools package."
Please note that I use Maple 13, and you made your worksheet using newer version, therefore your code is essentially useless for me because with CodeTools package and usage command I have to find the roots. Please give me a command that can work with Maple 13.
Anyway, thanks again.