Error in Ode23 line 114

8 ビュー (過去 30 日間)
JYOTI PRAKASH BEHERA
JYOTI PRAKASH BEHERA 2020 年 11 月 9 日
回答済み: Stephan 2020 年 11 月 9 日
Error in odearguments(FcnHandlesUsed, solver_name, ode, tspan, y0, options, varargin);
error in [tsol,hsol] =ode23(@(t,h) TnkODE(t,h), [0 2000],[5 5]);
Please Help
The code::
[tsol,hsol] =ode23(@(t,h) TnkODE(t,h), [0 2000],[5 5]);
function Diff=TnkODE(t,h)
h1=h(1); h2=h(2); % given values F0=3 Beta1=0.75, Beta2=1 A1=A2=0.5 are used futher in the equation
Diff(1)=6- 1.5*(h1-h2)^0.5 ; % h1'= F0/A1 - [(4/A1) (h1-h2)^.5]
Diff(2)=1.5*(h1-h2)^0.5 - 2*(h2)^0.5; % h2'= [(Beta1/A2) (h1-h2)^.5 ]- (Beta2/A2)h2^0.5
end

採用された回答

Stephan
Stephan 2020 年 11 月 9 日
[tsol,hsol] =ode23(@(t,h) TnkODE(t,h), [0 2000],[5 5]);
plot(tsol,hsol)
function Diff=TnkODE(~,h)
h1=h(1); h2=h(2); % given values F0=3 Beta1=0.75, Beta2=1 A1=A2=0.5 are used futher in the equation
Diff = zeros(2,1);
Diff(1)=6- 1.5*(h1-h2).^0.5 ; % h1'= F0/A1 - [(4/A1) (h1-h2)^.5]
Diff(2)=1.5*(h1-h2).^0.5 - 2*(h2).^0.5; % h2'= [(Beta1/A2) (h1-h2)^.5 ]- (Beta2/A2)h2^0.5
end

その他の回答 (0 件)

カテゴリ

Help Center および File ExchangeNumerical Integration and Differential Equations についてさらに検索

製品


リリース

R2020b

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by