4楼: Originally posted by 月只蓝 at 2015-09-08 09:56:55
MATLAB代码：
function solve_eqs
clear all;clc
global x
y0=1e-5;
x_data=linspace(1e-10,1e-6,100);
n=length(x_data);
options=optimset('TolFun',1e-12,'TolX',1e-15,'MaxFunEvals',1000);
for i=1:n
...

function solve_eqs
clear all;clc
global x
y0=1e-5;
x_data=linspace(1e-10,1e-6,100);
n=length(x_data);
options=optimset('TolFun',1e-12,'TolX',1e-15,'MaxFunEvals',1000);
for i=1:n
    x=x_data(i);
  
   [ y(i) fval(i)]=fsolve( @(y) log10(256*x^2*y^5*(6*y-x)^4)+( -13*( sqrt(13*y)/(1+sqrt(13*y)) -3.9*y) )+29.259637310505752, y0,options   );
end
[x_data; y;fval]'
figure(1)
subplot(1,2,1)
plot(x_data,y,'b-o');xlabel('x'),ylabel('y'),


subplot(1,2,2)
loglog(x_data,y,'r-');xlabel('log10(x)'),ylabel('log10(y)')