% PHD_PQvT  List P/Q ratio for uplink field

% this program was not used to produce any graphs in the thesis

format long
format compact
D=5;
T=0;

theta = atan(D);
Q0 = 2*cos(theta)^2 - sin(theta)^2;
P0 = 3*cos(theta)*sin(theta);

ratio0 = P0/Q0

    disp('T ratio P')
for T=0:1:8
    P = abs( PQintegral(D,T,0,1) );
    Q = abs( PQintegral(D,T,0,3) );
    ratio = P/Q;
    disp([T ratio P])
end

x=0:0.01:20;
plot(abs(sommerfeld(x,D,T,0,1)),'-')
