
function I=integral_evaluation(x,y,f_pts,rep_r,r_weights,theta_weights,degree_GL)

Lr=degree_GL+1;
Ltheta=Lr;
 
integrand=rep_r.*f_pts;       % INTEGRAND IN POLAR COORDINATES.
integrand=reshape(integrand,Lr,Ltheta);

I_inner=r_weights'*integrand; % INNER INTEGRAL IN "[0,1]", I.E.
                              % RADIAL COMPONENT.

I=I_inner*theta_weights; 