清楚的
信谊x1x2x3x4
xx = (X1, X2, X3, X4)
X1=xx(1);
d =长度(xx);
term1 = (X1-1) ^ 2;
金额= 0;
为2 = 2: d
席=XX(II);
xold=xx(ii-1);
ii * (2*xi^2 - xold)^2;
Sum = Sum + new;
结束
xc1之前(1)= 1;
xc2 (1) = 1;
xc3 (1) = 1;
xc4 (1) = 1;
呃= 10 ^ (2)
yb4 = (term1 +总和)
dyb4_x1 = diff (yb4 X1)
dyb4_x2=diff(yb4,x2)
dyb4_x3=diff(yb4,x3)
dyb4_x4 = diff (yb4 X4)
阁下=[潜艇(dyb4_x1 (X1, X2, X3, X4) [xc1之前(1)xc2 (1) xc3 (1) xc4(1)])潜艇(dyb4_x2, (X1, X2, X3, X4) [xc1之前(1)xc2 (1) xc3 (1) xc4(1)])潜艇(dyb4_x3, (X1, X2, X3, X4) [xc1之前(1)xc2 (1) xc3 (1) xc4(1)])潜艇(dyb4_x4, (X1, X2, X3, X4) [xc1之前(1)xc2 (1) xc3 (1) xc4 (1))))
S4 = -(阁下)
i = 1;
虽然norm(J4)>err
当= [xc1之前(我),xc2(我),xc3(我),xc4 (i)] ';
信谊h;
gr=subs(yb4[X1,X2,X3,X4],[xc1(i)+h*S4(1),xc2(i)+h*S4(2),xc3(i)+h*S4(3),xc4(i)+h*S4(4)])
h dgr_h = diff (gr)
h = vpasolve (dgr_h = = 0 h[无穷,正])
i=i+1
xc1之前(i + 1) =预告(1)+ h * S4 (1)
xc2 (i + 1) =预告(2)+ h * S4 (2)
xc3(i+1)=I4(3)+h*S4(3)
xc4(i+1)=I4(3)+h*S4(4)
结束
错误在迪克森(第38行)xc1之前(i + 1) =预告(1)+ h * S4 (1)