有关于数学模型Y*Exp(x/T0)=A0+B*Exp[(1/T0-1/T)*x]求解算法？

有关于数学模型Y*Exp(x/T0)=A0+B*Exp[(1/T0-1/T)*x]求解算法吗？Y 因变量，x自变量，A0，B，T0，T待求系数！对于式中所描述的函数曲线，可采用逐次逼近最小二乘法的非线性回归程序进行拟合。对测量时间较短的一组（Xi,Yi）值按上式Y*Exp(x/T0)=A0+B*Exp[(1/T0-1/T)*x]进行拟合， 确定A0，B，T0，T待求系数拟合及的具体步骤如下：
设T0为某个初始值，步长为dT0 ，则分别在T01＝T0+dT0,T01＝T0+2*dT0,T01＝T0+3*dT0对上式进行拟合，并算出每次的剩余平方和S1，S2 ，S3 。比较S1，S2 ，S3  ， 使 T0沿着S 变小的方向增大或减小，直到某个T0 时的S 大于T0+dT0 时的值为止，并设为新的初值，并设新的步长dT0=dT0/2,再重复上述步骤，直到满足给定精度e>(Sn-Sn-1)为止，此时A0，B，T0，T 即为最终结果。
这几弄得我头都大啦！希望那位高人能给我点算法提示，有类似相关算法代码好就更好？
1stOpt中的BFGS准牛顿法以上函数拟后效果很好！但我想自己搞一个delphi程序进行拟合！