3楼: Originally posted by 月只蓝 at 2015-02-04 20:35:45
以矩阵A=为例，代码如下（以下代码需完全复制进一个新建的m文件中，运行即可）：

function solve_nonlq
clear all;clc
x0=2;
options=optimset('TolFun',1e-9,'TolX',1e-9, 'MaxIter',500,'MaxFunEvals',1000 ...

function solve_nonlq
clear all;clc
x0=2;
options=optimset('TolFun',1e-9,'TolX',1e-9, 'MaxIter',500,'MaxFunEvals',1000);
[x fval]=fsolve(@obj,x0,options)


function y=obj(x)
A=[-2 1 3 ;3 1 7;besselj(1,x) 2 9];
y=det(A);