この作品は、たとえ営利目的であっても、許可を得ずに複製、改変・翻案、配布、上演・演奏することが出来ます。
http://creativecommons.org/publicdomain/zero/1.0/deed.enCC0Creative Commons Zero, Public Domain Dedicationfalsefalse
Scilab source
This media was created with Scilab, a free open-source software.
Here is a listing of the Scilab source used to create this file.
clear;//Parametersk=0.1;B=12.0;// Forced Duffing equation, x=x(1), y=x(2)functiondx=DUFF(t,x)dx(1)=x(2);dx(2)=-k*x(2)-x(1)^3+B*cos(t);endfunction// The time range and step for calculationTin=0;Tout=160;Tstep=0.01;t=Tin:Tstep:Tout;// The initial valuesx0=3;y0=0.1;X0=[x0;y0];X=ode(X0,0,t,DUFF);scf;plot(t(:),X(1,:))xlabel('t')ylabel('x')scf;plot(t(:),X(2,:))xlabel('t')ylabel('y')