逐次逼近最小二乘法的非线性回归程序汇编源码转Delphi源码，愿意重金答谢(100)

有关于数学模型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 即为最终结果。现请教那位高人赐教上述这种算法代码？