function gewdbi = rphv(ant,freq,elev,ax)
% gewdbi = rphv(ant,freq,elev,ax)
% Horizontal and vertical radiation diagram
%
% ant  :  antenna description.  There are three options:
%         MATLAB expression  :   description in MADL
%         sring :   file with description in NEC-syntax
%         {string; p} : string = file with description in AAL-syntax; 
%                       p = antenna parameters
%
% freq :  frequency in MHz, in case of a column vector the diagrams will be 
%         drawn for all frequncies in different figures, in case of a row 
%         vector all in one figure. A matrix is possible too.
%
% elev :  elevation angle in degrees for the horizontal diagram (default = 0).
%         If ground parameters are used, there is zero radiation at elevation 0.
%         In the case of HF-antennas  elev = 15  should be used.
%         If a vector of elevations is given all the curves will bedrawn into 
%         the same diagram. This feature is not possible together with 
%         freq being a vector or a matrix.
%
% ax   :  vector of limits of the diagram in dBi:
%         ax = [max_forward  max_backward  max_sidelobe].
%         If ax is not used or if ax = 0 two full diagrams will be presented.
%         If ax has length 1 and if ax~=0 the axes will appropriately chosen.
%
% gewdbi: gain over isotropic radiator (dB) 
%
% rphv erzeugt eine Inputdatei ant_temp.dat fr NEC2, ruft NEC2 auf 
% und sucht in der Ausgabetextdatei ant_temp.txt von NEC2 nach dem 
% Richtdiagramm.
%
% Erfolgt eine Zuweisung von gewdbi an eine Variable, so wird kein 
% Bild gezeichnet. 
% Andernfalls wird das Richtdiagramm grafisch dargestellt.
% Falls das Bild des Vertikaldiagramms zu klein ausfllt, liegt dies 
% an einer Nullstelle, die im 3D-Bild besser zu sehen ist.
% Das Horizontaldiagramm wird bei Bercksichtigung von Bodeneffekten zu
% Null. Dann ist elev zu vergrssern (oder ein 3D-Diagramm zu berechnen).
%
% Das Richtiagramm wird in einem 2 Grad-Raster berechnet. 
% Siehe auch  rphv1  mit 1 Grad-Raster.
