thetavals = -1:0.01:1;
A = [0.96058 0.82596 0.576 0.29171 0.12276 0.12276 0.29171 0.576 0.82596 0.96058];
u = sin(thetavals);
syms x U real
fun = A*exp(1j*pi*30*U*x)
Eint = int(fun,x,-1,1)
E0 = limit(Eint, U, 0);
E0idx = find(u == 0);
u(E0idx) = nan;
E = subs(Eint, U, u(:));
E(E0idx,:) = E0;
M = max(E, 2);
Y = abs(E);
G = Y./M;
Edb = 20*log10(G);
plot(thetavals, Edb);



