%用牛頓-雅可比迭代法求非線性方程組的一個(gè)根
function?[rm]=mulMix（Fx0hleps）
if?nargin==4
????eps=1.0e-4;
end

n?=?length（x0）;
J?=?zeros（nn）;
Fx?=?subs（Ffindsym（F）x0）;
for?i=1:n
????x1?=?x0;
????x1（i）?=?x1（i）+h（i）;
????J（:i）?=?（subs（Ffindsym（F）x1）-Fx）/h（i）;
end
D?=?diag（diag（J））;
C?=?D?-?J;
inD?=?inv（D）;
H?=?inD*C;
Hm?=?eye（nn）;
for?i=1:l-1
????Hm?=?Hm?+?power（Hi）;
end
dr?=?Hm*inD*Fx;
r?=?transpose（x0）-dr;???????????????????????????????
m=1;
tol=1;
while?tol>eps
????x0=r;???
????Fx?=?subs（Ffindsym（F）x0）;
????J?=?zeros（nn）;
????for?i=1:n
????????x1?=?x0;
????????x1（i）?=?x1（i）+h（i）;
????????J（:i）?=?（subs（Ffindsym（F）x1）-Fx）/h（i）;
????end
????D?=?diag（diag（J））;
????C?=?D?-?J;
????inD?=?inv（D）;
????H?=?inD*C;
????Hm?=?eye（nn）;
????for?i=1:l-1
????????Hm?=?Hm?+?power（Hi）;
????end
????dr?=?Hm*inD*Fx;????
????r?=?x0-dr;?????????????????????????????????%核心迭代公式
????tol=norm（r-x0）;
????m=m+1;
????if（m>100000）??????????????????????????????????????????????%迭代步數(shù)控制
????????disp（‘迭代步數(shù)太多，可能不收斂！‘）;
????????return;
????end
end

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----
?????文件????????1119??2017-11-23?20:08??mulMix.m