如何用matlab画超越方程组的图像,【求助】超越方程组的fsolve解法

【求助】超越方程组的fsolve解法

我用matlab的fsolve函数求解三元超越方程组，但是得到的解和文献上不同，困扰几个月了，在其他的论坛也求助了，但是一直都找不到解决的办法，现在课题全卡住了，现将代码和公式附上，请高手帮忙看看,万分感谢！sx1,sx2,tc是未知数，在程序运行时部分解输出的exitflag<0。

%%%%%%%%%%主程序%%%%%%%%%%%%%%%%%%%%%%

clear all;clc;close all;

global js;

global q;

Tc =[];

Sx1=[];

Sx2=[];

%q_par = [0.9];

%q_par = [1.065];

q_par = [2.0];

js_par = 2.0:-0.1:0;

n0= [ unifrnd(0,0.4) unifrnd(0,0.4) 1.5];

for q = q_par

js = js_par(1);

options=optimset('display','on');

[n(1,:),fval,exitflag]=fsolve(@myfun,n0,options)

for i = 2:length(js_par)

js = js_par(i);

[n(i,:)