How do i merge graphs in bvp4c format ?
1 次查看(过去 30 天)
显示 更早的评论
What are a various ways that a loop can be created ? Also what is the use of 'global' ?
0 个评论
回答(1 个)
Stephan
2019-7-14
hold on
for k = 1:4;
Ha = k;
init = bvpinit(linspace(-1,1,20),[0 0]);
sol=bvp4c(@(x,y)rhs_bvp(x,y,Ha),@bc_bvp,init);
x = linspace(-1,1,20);
BS(:,:,k) =deval(sol,x);
plot(BS(1,:,k),x, '--');
xlim([-1 1])
end
hold off
function rhs=rhs_bvp(x,y,Ha)
rhs = [y(2); ((Ha^2+(1/3))/(1+0.3))*y(1)];
end
function bc=bc_bvp(yl,yr)
bc = [yl(1)-(0.1*1.3)*yl(2)-1; yr(1)+(0.1*1.3)*yr(2)+1];
end
0 个评论
另请参阅
产品
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!