Wolfram Cloud Document

ClearAll["Global`*"];

eq11a=x'[t]y[t];

eq11b=y'[t]lam(x-lamx^3+lam^2x^5)(x'[t])^2-alfgam(1+lamx^2)x'[t]-w0^2x-delx^3-xix^5+gam(1+lamx^2)(fCos[wt]+gCos[wwt])/.{x'[t]y[t],xx[t]};

val={alf.2,w0I,ww9.842,w.5,f.05,lam0,beta1,m01.5,gam1/m0,delbetagam-lamw0^2,xibetagamlam+(lamw0)^2};

eq11b2=eq11b//.val;

sol[gg_?NumericQ,t1_?NumericQ]:=Module[{},eq11b3=eq11b2/.ggg;NDSolve[{eq11a,eq11b3,x[0]1,y[0]1},{x,y},{t,0,t1}]]

n=1;w=0.5;T=(2*N[Pi])/w;Ts=n*T;

qQnum[g_]:=Module[{},Ts=n*T;lsg=sol[g,Ts];fx=x/.lsg[[1,1]];fy=y/.lsg[[1,2]];qQs=NIntegrate[fx[t](Sin[wt]/.val),{t,0,Ts}]2/Ts;qQc=NIntegrate[fx[t](Cos[wt]/.val),{t,0,Ts}]2/Ts;qnum=Sqrt[qQs^2+qQc^2]]

tt=Table[{g,qQnum[g]},{g,0,Ts}]

$Aborted

ListLinePlot[tt]

ReplaceAll::argx:ReplaceAll called with 2 arguments; 1 argument is expected.



NIntegrate::nlim:t = T is not a valid limit of integration.



ReplaceAll::argx:ReplaceAll called with 2 arguments; 1 argument is expected.



NIntegrate::nlim:t = T is not a valid limit of integration.



ReplaceAll::argx:ReplaceAll called with 2 arguments; 1 argument is expected.



General::stop:Further output of ReplaceAllargx will be suppressed during this calculation.



NIntegrate::nlim:t = T is not a valid limit of integration.



General::stop:Further output of NIntegratenlim will be suppressed during this calculation.

