%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
%  cratio.m - forms the ratio of products of complex numbers
%             and displays the result in both polar and
%             rectangular format.
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
disp(' ');
disp('  RATIO OF COMPLEX NUMBER PRODUCTS'); disp(' ');
disp('  Form: 1-polar, 2-rectangular'); disp(' ');
num=input('How many numerator numbers? '); disp(' ');
for i=1:num;
   F=input(['Form of ',num2str(i),' = ']);
   if (F~=1 & F~=2); disp('INVALID FORM'); end
   if F==1
      M=input(['Mag  ',num2str(i),' = ']);
      A=input(['Deg  ',num2str(i),' = '])*pi/180;
      N(i)=M*exp(j*A); disp(' ');
   else
      R=input(['Real ',num2str(i),' = ']);
      I=input(['Imag ',num2str(i),' = ']);
      N(i)=R+j*I; disp(' ');
   end
end
NP=N(1);  for k=2:num; NP=NP*N(k); end
den=input('How many denominator numbers? '); disp(' ');
for i=1:den;
   F=input(['Form of ',num2str(i),' = ']);
   if (F~=1 & F~=2); disp('INVALID FORM'); end
   if F==1
      M=input(['Mag  ',num2str(i),' = ']);
      A=input(['Deg  ',num2str(i),' = '])*pi/180;
      D(i)=M*exp(j*A); disp(' ');
   else
      R=input(['Real ',num2str(i),' = ']);
      I=input(['Imag ',num2str(i),' = ']);
      D(i)=R+j*I; disp(' ');
   end
end
DP=D(1);  for k=2:den; DP=DP*D(k); end; RAT=NP/DP;
disp(' '); disp(['RATIO = ' num2str(real(RAT)) ' +j '...
      num2str(imag(RAT)) ' = ' num2str(abs(RAT)) '|_'...
      num2str(angle(RAT)*180/pi) 'deg']);