|
Debug了一会儿，不过还是搞定了。为了方便你阅读，也为了方便debug，我写成了嵌套函数的形式。最外层的函数就是myfun，你只要输入特定的x数值，特定的初始值y0,ma,mb,ms（这些你都自己定，随便；ms是迭代步长，ma和mb是变量m的起点和终点，起点终点越远，步长越小，迭代越花时间，但图形越精确）function output = myfun(x,y0,ma,mb,ms)m=[ma:ms:mb];for i=1:length(m)y(i)=fsolve(@(y) myfun1(x,y)-m(i),y0,optimset('Display','off'));z(i)=myfun2(x,y(i));endoutput.y=y;ou... |
