Labfans是一个针对大学生、工程师和科研工作者的技术社区。 | 论坛首页 | 联系我们(Contact Us) |
![]() |
![]() |
#1 |
初级会员
注册日期: 2010-12-20
年龄: 35
帖子: 1
声望力: 0 ![]() |
![]()
for var1=0:0.001:15
var2=var1/R1; [x,y]=ode23('fun0088',[0,60],[var1;var2]); pp=(-p0.*(1+R1)+sqrt(p0^2.*(1+R1)^2+4.*R1.*y(:,1).*(p0+R2.*y(:,1)+p0.*R2)))/(2.*R1); ppp=(-p0.*(1+R2)+sqrt(p0^2.*(1+R2)^2+4.*R2.*pp.*(p0+R1.*pp+p0.*R1)))/(2.*R2); pppp=ppp*R2; if abs(mean((y(:,2)-pppp)))<1; break end end |
![]() |
![]() |