Matlab 泊松回归的问题?

大家好，现在我想用Matlab 解决泊松回归的问题。有人能解决泊松回归的问题吗？就像Matlab 中 curve fit中高斯回归一样，不过matlab 没有现成函数，自己编也没编出来，网上查也没找到合适的，一阶就可以了。Matlab中custom function无法使用阶乘。问题如下：问一列x:x1,x2,x3.........,y:y1,y2,y3............... 找到方程y=a*b^((x-c)/d))*exp(-b)/((x-c)/d))!,使相关系数最大。