solving the equation using rk-4 method, but couldn't get the right output, last time I forget to add the program please check this one.
1 visualización (últimos 30 días)
Mostrar comentarios más antiguos
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/1072310/image.png)
i want to solve this, where σ = 0.1, L0 = 0.5, k= 0.01, tc = 70E-9
format long
a = 0;
b = 4;
h = 0.1;
k = 0.01; % critical coupling strength
L = 0.5;
s = 0.1;
tc = 70E-9; % photon life time in cavity
n = (b-a)/h ;
t = a + (0:n)*h;
func1 = @(t,y) [(-(s^2)*k/tc)*sin(y - pi/2) + L*(s^2)/(2*tc)*sin(y + pi/2)/sqrt(1 + cos(y + pi/2))];
for i = 1:n
k1 = func1(y(i));
k2 = feval(func1, y(i)+(k1/2)*h);
k3 = feval(func1, y(i)+(k2/2)*h);
k4 = feval(func1, y(i)+(k3*h));
y(i+1) = y(i) + (h/6)*(k1+(2*(k2+k3))+k4);
end
plot(t,y)
7 comentarios
Walter Roberson
el 21 de Jul. de 2022
You should preallocate y as length n+1. With your current code if you just happen to have a longer y in your workspace you would have a problem.
Respuestas (0)
Ver también
Categorías
Más información sobre Ordinary Differential Equations en Help Center y File Exchange.
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!