實現主代碼
n = 1;
while (norm(x - x0)>tol)& amp;& amp(n & lt1000)
x0 = x;
x = x0-A \ funm(x0);
p = x-x0;
q = funm(x)-funm(x0);
a = a+(q-a*p)*p'^rm(p)^2;
n = n+1;
結束