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