非線形自由振動 初期条件 とする解析解 clear; close; t=0:0.1:10 omega=1 e=0.1 for a=[1,5,10,11,12] x = a*cos(omega*t)+... e*a^3/(32*omega^2)*(cos(3*omega*t)-cos(omega*t))+... e^2*a^5./(1024*omega^4)*(cos(5*omega*t)-cos(omega*t))+... e^3*a^7/…