10楼: Originally posted by SEEDEN at 2013-04-26 10:35:37
fzero可以求解无解析式解的代数方程的符号解吗？

我个人觉得你单独使用这个命令或此类命令恐怕都不能实现你的目的，当然，你可以再多试试几个命令。

这个代数方程组就算有解，你确定只有一组解？而且具体是哪 ...

2楼: Originally posted by 月只蓝 at 2013-04-25 16:25:30
你用MOL法求解的本身就是一个常微分方程组，原方程组有刚性，ode15s能解出，那么双组份的情形，无非是把方程组的维数提高一倍。

9楼: Originally posted by 小木虫zb at 2013-04-26 09:09:40
我修改了下程序，但是还是有很多错误，
function dydt=fangcheng(t,y)
kf=200;
a=0.363;%空隙率
u=0.361032/a;%床层颗粒间隙中流体流速
c0=7.03;
L=8.31e-2;
N=20;
dz=L/N;
p=436.8/(1-0.363);%颗粒表观密 ...