| 查看: 1016 | 回复: 0 | |||
[交流]
Matlab求解微分方程bvp4c
|
|
solinit = bvpinit(linspace(0,14,10),[1;0]); sol=bvp4c(@ODEfun,@BCfun,solinit); x = linspace(0,14,127); y = deval(sol,x); plot(x,y(2,) function dydx=ODEfun(x,y) a=pi*14.701^4/64; b=pi*(24^4-14.701^4)/64; c=3.54365802196686; d=5.94058826800217; e=2/14.701/pi; f=14.701/2; dydx(1,1)=y(2); dydx(2,1)=(1/a+1/b)/(c+d)*f^2/e*(cos(x/(f*e)))^2*y(1)-(tan(x/(f*e)))/f/e*y(2); function res=BCfun(ya,yb) res=[ya(1); yb(1)-1]; 哪位大神给看看我这个微分方程的数值解是多少。用bvp4c解出来的总感觉不对 |
回复此楼» 猜你喜欢
A期刊撤稿
已经有3人回复
-
职称评审没过，求安慰
已经有34人回复
-
垃圾破二本职称评审标准
已经有17人回复
-
回收溶剂求助
已经有6人回复
-
投稿Elsevier的Neoplasia杂志，到最后选publishing options时页面空白，不能完成投稿
已经有22人回复
-
申请26博士
已经有5人回复
-
EST投稿状态问题
已经有7人回复
-
毕业后当辅导员了，天天各种学生超烦
已经有4人回复
-
求助文献
已经有3人回复
-
投稿返修后收到这样的回复，还有希望吗
已经有8人回复