高斯拟合算法求助

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