Can someone help me generate numerical values for this sample problem using bvp4c solver?
1 回表示 (過去 30 日間)
古いコメントを表示
y'' + y = 0
Boundary conditions are y(0)=1; y(pi/2)=0.
Compute y(pi/3) and y(pi/4) using MATLAB Bvp4c solver.
Can someone help me generate numerical values for this sample problem using bvp4c solver?
I want code command to generate the values of y(pi/3) and y(pi/4).
Trial code
function sampletrial
clear all;
close all;
clc;
gva=[pi/3 pi/4];
for i=1:2
g=gva(i);
figure(1);
solinit=bvpinit(linspace(0,pi/2,50),[1 0]);
sol=bvp4c(@shootode5,@shootbc5,solinit);
eta=sol.x;
f=sol.y;
plot(eta,(f(1,:)),'--','linewidth',2);
hold on
end
function dydx=shootode5(eta,f)
dydx=[f(2)
-f(1)];
end
function res=shootbc5(fa,fb)
res=[fa(1)-1
fb(1)];
end
0 件のコメント
採用された回答
Rajeev
2023 年 3 月 6 日
移動済み: Image Analyst
2023 年 3 月 6 日
Can you try adding the line:
y_values(i) = interp1(sol.x, sol.y(1,:), g);
in the loop?
その他の回答 (0 件)
参考
カテゴリ
Help Center および File Exchange で Boundary Value Problems についてさらに検索
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!