有100组数据. x与y的关系是一个非线性函数（一条曲线），y=(a+x/b)*exp((-x/b)*c)，现在想要求出每一个id的未知系数a,b,c。每个id有5个分散的x, 以及对应的y。并且求出曲线下的面积，不用梯形法则。

用proc iml怎么实现啊，自定义公式不太会用。
proc nlin已经试过了，很多id的收敛失败了，并且求出的系数好夸张，看着根本不对，重新带入x,y后，算出来的y也对不上。