b'\n \n \n
 
\n  
sfduffing (4.0)
index
user/yliu/Mduffing.c
\n Duffing differential equation solved by 4th order Runge-Kutta method. \n

\n \n \n \n \n \n
 
\n Synopsis
       sfduffing < in.rsf > out.rsf sfile=sfile.rsf gamma=0.75 omega=1 kxi=1 x0=0 y0=0 pow1=1 pow2=3 verb=n ricker=n
Duffing equation: x\'\' + 0.5 x\' - x + x^3 = gamma cos(omega t) + kxi input(t)
\n\n

\n \n \n \n \n \n
 
\n Parameters
       \n \n \n
\n  
float gamma=0.75
\tstrength of external force
\n
\n \n\n \n \n
\n  
float kxi=1
\tadjustment for input signal
\n
\n \n\n \n \n
\n  
float omega=1
\tangular frequence of external force
\n
\n \n\n \n \n
\n  
int pow1=1
\tpower of first non-linear restitution term
\n
\n \n\n \n \n
\n  
int pow2=3
\tpower of second non-linear restitution term
\n
\n \n\n \n \n
\n  
bool ricker=n [y/n]
\tif y need extenal input for external force
\n
\n \n\n \n \n
\n  
string sfile=
\tauxiliary input file name
\n
\n \n\n \n \n
\n  
bool verb=n [y/n]
\tverbosity flag
\n
\n \n\n \n \n
\n  
float x0=0
\tinitial value of x0
\n
\n \n\n \n \n
\n  
float y0=0
\tinitial value of y0
\n
\n \n
'