# Question:plot output of an integral

## Question:plot output of an integral

Maple 2018

how this integration can be computed and plotted? tnx in advance

 > restart:Digits:=20:
 > T := 10^(-6); delta := 10; m := 1; alpha:= 75/10; v := 65; mu := 8;
 (1)
 > z := sqrt((omega^2 - 4*m^2)/alpha^2 - v^2/alpha^2*(4 *delta/alpha*tan(theta)^2 - v^4/alpha^4*tan(theta)^4))
 (2)
 > a := sec(theta)/sqrt(2)*sqrt((2 *delta)/alpha - v^2/alpha^2*tan(theta)^2 + z)
 (3)
 > f1 := 1/(exp(( omega/2 - mu)/T) + 1); f2 := 1/(exp((- (omega/2) - mu)/T) + 1);
 (4)
 > A1 := evalf(-(f1 - f2)*abs(omega/(4*a*alpha*cos(theta)^2*(-delta+ a^2 *alpha* cos(theta)^2) + 2*a*v^2*sin(theta)^2 + 1/10000)))
 (5)
 > b := sec(theta)/sqrt(2)*sqrt((2*delta)/alpha - v^2/alpha^2*tan(theta)^2 - z)
 (6)
 > f3 := 1/(exp(( omega/2 - mu)/T) + 1); f4 := 1/(exp((- (omega/2) - mu)/T) + 1);
 (7)
 > A2 := evalf(-(f3 - f2)*abs(omega/(4*b*alpha*cos(theta)^2*(-delta+ b^2 *alpha* cos(theta)^2) + 2*b*v^2*sin(theta)^2 + 1/10000)))
 (8)
 > plot(int(A1+A2,theta=0..2*Pi),omega=0..50);
 >