[R_EQU,FL]=shape(PNr)

function [R_EQU,FL]=shape(PNr)

% function [R_EQU,FL]=shape(PNr)
%
% equatorial radius and flattening of the planets
% algorithms by Montenbruck/Pfleger: "Astronomy with the Personal Computer"
%
% 05.05.04 M.Penzkofer

% check input parameter
if (PNr < 1 | PNr > 9)
   disp('Enter number between 1 and 9!');
   return
end

switch PNr
	case 1	% Mercury
      R_EQU =  2439.0;
      FL = 0.0;
   case 2	% Venus
      R_EQU =  6051.0;
      FL = 0.0;
   case 3	% Earth
      R_EQU =  6378.14;
      FL = 0.00335281;
   case 4	% Mars
      R_EQU =  3393.4;
      FL = 0.0051865;
   case 5	% Jove
      R_EQU = 71398.0;
      FL = 0.0648088;
   case 6	% Saturn
      R_EQU = 60000.0;
      FL = 0.1076209;
   case 7	% Uranus
      R_EQU = 25400.0;
      FL = 0.030;
   case 8	% Neptune
      R_EQU = 24300.0;
      FL = 0.0259;
   case 9	% Pluto
      R_EQU =  1500.0;
      FL = 0.0;
end