# solving a second order linear differential equation

24 ビュー (過去 30 日間)
Eliraz Nahum 2018 年 9 月 23 日
Answered: Eliraz Nahum 2018 年 9 月 23 日
hello everybody, I was trying to solve a simple pendulum second order linear differential equation of the form y''=-(g/l)*sin(y) while using the ode45 function. since it's a second order equation I understood that I have to manipulate the problem, so it will fit the ode45.
The mathematical manipulation I did is described in the attached picture.
g=9.8;
l=0.5;
t_span=[0 30];
F=@(t,q) [q(2);-(g/l)*sin(q(2))];
[t,q]=ode45(@(t,q) F(t,q),t_span,teta0)

#### 0 件のコメント

サインイン to comment.

### 採用された回答

Mischa Kim 2018 年 9 月 23 日

Almost there:
g = 9.8;
l = 0.5;
t_span = [0 30];
F = @(t,q) [q(2);-(g/l)*sin(q(1))]; % check your derivation
This is a second order DE so you need two initial conditions, one for teta and one for tetad.

#### 0 件のコメント

サインイン to comment.

Eliraz Nahum 2018 年 9 月 23 日
thank you very much... sometimes it's not about the code, but the mathematics :-)

#### 0 件のコメント

サインイン to comment.

サインイン してこの質問に回答します。

R2018a

Translated by