非线性方程组牛顿迭代法matlab,matlab实现牛顿迭代法求解非线性方程组

已知非线性方程组如下

3*x1-cos(x2*x3)-1/2=0

x1^2-81*(x2+^2+sin(x3)+=0

exp(-x1*x2)+20*x3+(10*pi-3)/3=0

求解要求精度达到————————————————————————————————首先建立函数fun

储存方程组编程如下将保存到工作路径中:

function f=fun(x);

%定义非线性方程组如下

%变量x1 x2 x3

%函数f1 f2 f3

syms x1 x2 x3

f1=3*x1-cos(x2*x3)-1/2;

f2=x1^2-81*(x2+^2+sin(x3)+;

f3=exp(-x1*x2)+20*x3+(10*pi-3)/3;

f=[f1 f2 f3]; ————————————————————————————————建立函数dfun

用来求方程组的雅克比矩阵将保存到工作路径中:

function df=dfun(x);

%用来求解方程组的雅克比矩阵储存在dfun中

f=fun(x);

df=[diff(f,'x1');diff(f,'x2');diff(f,'x3')];

df=conj(df');————————————————————————————————编程牛顿法求解非线性方程组将保存到工作路径中: