Can someone help me generate numerical values for this sample problem using bvp4c solver?
3 visualizzazioni (ultimi 30 giorni)
Mostra commenti meno recenti
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 Commenti
Risposta accettata
Rajeev
il 6 Mar 2023
Spostato: Image Analyst
il 6 Mar 2023
Can you try adding the line:
y_values(i) = interp1(sol.x, sol.y(1,:), g);
in the loop?
Più risposte (0)
Vedere anche
Categorie
Scopri di più su Boundary Value Problems in Help Center e File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!