Integration of a list
- To: mathgroup at yoda.physics.unc.edu
- Subject: Integration of a list
- From: bappadit at ecn.purdue.edu (Banerjee Bappaditya)
- Date: Mon, 11 Jan 93 19:59:44 EST
Dear Mathgroupers,
I have been trying to Integrate a moderate sized list from -Pi to Pi
and it is taking forever.....the individual components of the expression
are quite simple, so I cannot figure out what is taking so long...
test =
{0, ((-0.25*I*E^(-3.*I*Omgn*Pi)*x1*y1)/(Omgn^3*Pi) +
(0.75*I*E^(-1.*I*Omgn*Pi)*x1*y1)/(Omgn^3*Pi) +
(-0.75*I*E^(1.*I*Omgn*Pi)*x1*y1)/(Omgn^3*Pi) +
(0.25*I*E^(3.*I*Omgn*Pi)*x1*y1)/(Omgn^3*Pi) -
(3*E^(-I*Omgn*t)*(I + -I*E^(2*I*Omgn*t))*
((E^(-I*Omgn*t)*(1 + E^(2*I*Omgn*t))*x1)/2 +
(E^(-I*Omgn*t)*(I + -I*E^(2*I*Omgn*t))*y1)/(2*Omgn))^2)/(2*Omgn))*
((0.04166666666666666*I*E^(-3.*I*Omgn*Pi)*t*x1^3)/(Omgn*Pi) +
(0.375*I*E^(-1.*I*Omgn*Pi)*t*x1^3)/(Omgn*Pi) +
(-0.375*I*E^(1.*I*Omgn*Pi)*t*x1^3)/(Omgn*Pi) +
(-0.04166666666666666*I*E^(3.*I*Omgn*Pi)*t*x1^3)/(Omgn*Pi) +
(-0.125*I*E^(-3.*I*Omgn*Pi)*t*x1*y1^2)/(Omgn^3*Pi) +
(0.375*I*E^(-1.*I*Omgn*Pi)*t*x1*y1^2)/(Omgn^3*Pi) +
(-0.375*I*E^(1.*I*Omgn*Pi)*t*x1*y1^2)/(Omgn^3*Pi) +
(0.125*I*E^(3.*I*Omgn*Pi)*t*x1*y1^2)/(Omgn^3*Pi) -
E^(I*Omgn*t)*(((-3*I)/8*x1^3)/Omgn - (3*x1^2*y1)/(8*Omgn^2) +
((-3*I)/8*x1*y1^2)/Omgn^3 - (3*y1^3)/(8*Omgn^4)) +
(E^(-3*I*Omgn*t)*((-I/8*x1^3)/Omgn + (3*x1^2*y1)/(8*Omgn^2) +
((3*I)/8*x1*y1^2)/Omgn^3 - y1^3/(8*Omgn^4)))/3 -
(E^(3*I*Omgn*t)*((-I/8*x1^3)/Omgn - (3*x1^2*y1)/(8*Omgn^2) +
((3*I)/8*x1*y1^2)/Omgn^3 + y1^3/(8*Omgn^4)))/3 +
E^(-I*Omgn*t)*(((-3*I)/8*x1^3)/Omgn + (3*x1^2*y1)/(8*Omgn^2) +
((-3*I)/8*x1*y1^2)/Omgn^3 + (3*y1^3)/(8*Omgn^4)) +
(0.5*F*t*Sin[-1.*Omg*Pi])/(Omg*Pi) - (0.5*F*t*Sin[1.*Omg*Pi])/(Omg*Pi) +
(F*Sin[Omg*t])/Omg)}
Tperiod = 2*Pi
Mean[func_List,t_]:=
(1.0/Tperiod)*Integrate[func,{t,-Tperiod/2.0,Tperiod/2.0}]
Mean[test,t] is taking forever. Any ideas ?
Thanks,
bappa.
Bappaditya Banerjee
bappadit at mn.ecn.purdue.edu
Ray W. Herrick Laboratories
Purdue University
West Lafayette, IN 47907
work : (317) 494 2132
(317) 494 2147
fax : (317) 494 0787
home : (317) 743 3982