Non linear differential equation of 2 degrees of freedom

Equation of motion;
d2Xa/dt2=(-Ka(Xa-Xb)-Ca(dXa/dt-dXb/dt))/Ma
d2Xb/dt2=(( Ka(Xa-Xb)+Ca(dXa/dt-dXb/dt) -Kb(Xb-Xc)-Cb(dXb/dt-dXc/dt))/Mb
0= Xc(t)=X0/△t*t
△t = Xc(t)=X0
Procedure for solution;
d2Xa/dt2=((-Ka(Xa-Xb)-Ca(dXa/dt-dXb/dt))/Ma
d2Xb/dt2=(( Ka(Xa-Xb)+Ca(dXa/dt-dXb/dt) -Kb(Xb-Xc)-Cb(dXb/dt-dXc/dt))/Mb
dXa(t+dt)/dt=dXa(t)/dt+d2Xa(t)/dt2*dt
dXb(t+dt)/dt=dXb(t)/dt+d2Xb(t)/dt2*dt
Xa(t)=Xa(t)+dXa(t)/dt*dt
Xb(t)=Xb(t)+dXb(t)/dt*dt

※=+$D$5*B9
※1=@IF(+D9<$E$5,+$F$5/$E$5*D9,$F$5)
※2=(-$I$5*(J9-K9)-$K$5*(H9-I9))/$G$5
※3=($I$5*(J9-K9)+$K$5*(H9-I9)-$J$5(K9-E9)-$L$5*(I9-(E10-E9)/$D$5))/$H$5
※4=+H9+F9*$D$5
※5=+I9+G9*$D$5
※6=+J9+H9*$D$5
※7=+K9+I9*$D$5
Acceleratoin of Ma


お名前をどうぞ(Your name,please):
E-mail addressもどうぞ(please)
一言コメント(comment)を:

頭を休めるために、時々御立ち寄り下さい。花の名前も覚えられます。
きっとデートのときに役に立ちます。
庭の花 Garden flower

私のささやかなコレクション
世界の栞 ; Bookmark

腰痛、痔、水虫でお悩みの方、そして酒の好きな人,是非お寄りください。
風邪の温水うがい治療法も紹介します
私の健康法