先说两个问题: 1、你的kk在个函数中都作为形式参数传递进去,没必要也不应该设为global(而且较新的版本可能会不支持这种用法); 2、你的函数myf返回3x1的向量,dmyf返回3x3的矩阵,下面这个式子肯定是没法计算的:r=x0-myf(x0,kk)/dmyf(x0,kk); 到底怎么改,我不知道你的用意,只能先提醒你一下,等你把这两个问题搞清楚了,如果还有其它问题,就继续追问吧。
请到我的百度空间,标题为 To lyuhxm