Third Order Boundary Value Problem (BVP) with Boundary Coditions only in Terms of Derivatives
조회 수: 7 (최근 30일)
이전 댓글 표시
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
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)
답변 (0개)
참고 항목
카테고리
Help Center 및 File Exchange에서 Numerical Integration and Differential Equations에 대해 자세히 알아보기
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!