Third Order Boundary Value Problem (BVP) with Boundary Coditions only in Terms of Derivatives

11 ビュー (過去 30 日間)
Saeid
Saeid 2020 年 7 月 7 日
コメント済み: Saeid 2020 年 7 月 8 日
I need to solve the boundary value problem (related to flow of fluids in conical ducts):
Y''' + 1 - Y'^2 = 0
with all the three BCs bing only on the derivatives of Y (and none on Y itself):
Y'(0)=Y"(inf)=0
Y'(inf)=1
By formulating the problem exactly as suggested here:
dYdX = @(X,Y) [Y(2); Y(3); Y(2).^2-1]; % Differential equation
res = @(ya,yb) [ya(2); yb(3); yb(2)-1];
SolYinit = bvpinit([0 5], [0; 0; 0]);
Fsol = bvp4c(dYdX, res, SolYinit);
X = Fsol.x;
Y = Fsol.y;
I tried to solve the problem, but received the Jacobian error.
Is it possible that the error is generated because the BCs do not contain Y(...), and is there a way around his problem?
  4 件のコメント
darova
darova 2020 年 7 月 8 日
I mean this
Y'(0)=Y"(inf)=0
Y'(inf)=1
res = @(ya,yb) [ya(2); yb(3); yb(2)-1];
THere are only derivatives. There should be at least one Y(0) or Y(inf) (starting or final position)
Saeid
Saeid 2020 年 7 月 8 日
That's what I suspected, thanks! Then I guess the way the problem is posed in the original literature source is not right.

サインインしてコメントする。

回答 (0 件)

カテゴリ

Help Center および File ExchangeBoundary Value Problems についてさらに検索

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by