clear;
theta = linspace(-90, 90, 100);
phi = 0;
mname = "R";
res = 401;
farfieldfilter(0.5);
f = getdata(mname,"f");
E2=farfield3d(mname,1,res,res);
ux=farfieldux(mname,1,res,res);
uy=farfielduy(mname,1,res,res);

E2_polar = farfieldspherical(E2,ux,uy,theta,phi);
polar(theta/180*pi,E2_polar);
legend("xz plane");
plot(theta,E2_polar);