You may want to express Rate of Climb (ROC) as a function of altitude, and then integrate the expression over the altitude region which you are interested in.
Thank you very much, but I found it. The source of error may come from SFC with the thrust varying with the mach number: how do I take into account this? Thank you.