R中关于nls函数执行非线性回归的问题

很是奇怪，我使用R中的nls（）函数和Matlab中的nlinfit()函数对相同的数据集和相同的初值执行拟合，结果却不一样！！！
都是所谓的“高斯牛顿”法，为什么差距这么大呢？
请真正的高手赐教！

请教诸位高手一个问题：
自己编写了函数 f<-function(a,b,c) { }
调用时，现在a和b是必须录入的值，而c可录入也可以不录入，如果不录入c在程序中其值为500，而如果录入，则其值根据录入的值而定，该如何实现c值的编程呢？

第一个问题很难回答你，有两种可能：
1。我们不知道你到底是怎么用的，你R里怎么写的？matlab里怎么写的？会不会是你参数设置不当？我们不得而知……
2。不同软件包结果不同，这已经不是新闻了……基本上，我所有的东西都是自己写，只有自己做的东西最放心……

第二个问题这么写 f <- function(a,b,c=500){}

谢谢啦。
自己编写的东西最放心，的确如此。