本人想通过实现一维高斯拟合y=a*exp(-(x-b)^/c^2)，得到a，b，c的值，但是不知道如何去做。但是不能对它来求对数，因为y有为0的值，请教高手如何做，该如何实现？其中我使用过最小二乘法加上牛顿迭代和高斯主消元编写程序，但是由于奇异的原因不能迭代下去。
    高手帮帮忙吧！

目前我的解法是通过拟牛顿算法来实现这个过程。

数值计算类的东西忘了，不过这种程序教材上就应该有吧