Hello, I want to make function Runge-Kutta orde 4 with equation depend on time, but in MATLAB i get some error
2 ビュー (過去 30 日間)
古いコメントを表示
cindyawati cindyawati
2023 年 5 月 19 日
コメント済み: Cris LaPierre
2023 年 5 月 20 日
I want to use this equation with runge-kutta orde 4
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/1388514/image.png)
function M1 =fRK4M1(M1)
%M1 =M1(j)+1./(1+exp(-T(j)));
m1 = (50*M1(j+1)*(1-(M1(j+1)/75))-2*(10^-4)*M1(j+1)*M1(j+1)-M1(j+1)*((5*10^-4).*M2(j+1))-((10-3)*(10^-4)*M1(j+1)*M3(j+1))-((100-10)*0.1*M1(j+1)*O(j+1))-((10^-3)*M1(j+1)));
M1 = dt *m1;
end
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/1388519/image.png)
0 件のコメント
採用された回答
Cris LaPierre
2023 年 5 月 19 日
You have not defined j, so in MATLAB, by default that letter is used to create complex numbers. As the error message states, your index must be positive integers.
j+1
a=1:5;
a(j+1)
2 件のコメント
Cris LaPierre
2023 年 5 月 20 日
Keep in mind that MATLAB is case sensitive. T and t are different variables.
その他の回答 (0 件)
参考
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!