又看了下问题,补充一下
在我看来函数jiaoxian的自变量应该是u,函数值是交线方程的解,因此对fsolve的调用应该在jiaoxian内部吧。也就是
function F=jiaoxian(u)
...
F=fsolve(...);
...
简洁的来说,交线求解的代码结构应该是
clear all;
...
eq1=...;
eq2=...;
eq3=...;
F=[];
for ...
...
F=[F,fsolve(...)];
...
end
...
一个过程文件就可以了,不用定义函数的吧。
|