p=300;
q=115;
s=321;
disp("Part a");
pf=p/s;
disp("power factor of the combination is"); disp(pf);disp("lagging");
disp("Part b");
deg=acos(pf)*180/%pi;
disp("the phase angle (in deg) between the applied voltage and current is"); disp(deg);