DataMuseum.dk

top - metrics - download

    Length: 1536 (0x600)
    Types: TextFile
    Names: »errelltx«

└─⟦621cfb9a2⟧ Bits:30002817 RC8000 Dump tape fra HCØ.  Detaljer om "HC8000" projekt.
    └─⟦0364f57e3⟧ 
        └─⟦787c125fb⟧ »adjprocfile« 
└─⟦00964e8f7⟧ Bits:30007478 RC8000 Dump tape fra HCØ.
    └─⟦b2ec5d50f⟧ 
        └─⟦787c125fb⟧ »adjprocfile« 
            └─⟦this⟧ 

TextFile

;       err_ell_tx            * page 1    9 01 78, 14.03;  

;  err_ell
;  *******

if listing.yes
char 10 12 10

err_ell = set 1

err_ell = algol

external real procedure err_ell(mn, mne, me, a, b, dir);  
________________________________________________________
value                           mn, mne, me;  
real                            mn, mne, me, a, b, dir;  

comment

err_ell             (call and return, real)
The station mean error, assuming a mean error of unit weight = 1

mn, mne, me         (call, real)
The variance and covariance elements of the coordinates.

a, b                (return, real)
The major and minor sime half axes of the error ellipse.

dir                 (return, real)
The direction of the major axis counted from north, clockwize, 
in radians.

Prog.: Knud Poder, NOV 1977.

GI. lib nr  77 037.  

;  

begin

  real sum, dif, az2, cos_2az, sin_2az, nd2;  
   long azd;

  sum := 0.5*(mn + me);  
  dif := 0.5*(mn - me);  
   az2 := arg(dif,mne);
azd := az2/rg;
nd2 := dif*cos(az2) + mne*sin(az2); <*neg 2. deriv/2*>
if nd2 < 0 then
begin
nd2 := -nd2;
azd := azd + (extend 1) shift 47;
end;
a := sqrt(sum + nd2);
b := sqrt(sum - nd2);
azd := azd//2;
dir := azd*rg;
  err_ell := sqrt(sum);  

end err_ell;  

end

\f



;       err_ell_tx            * page 2    9 01 78, 14.03;  
 
if ok.no
mode warning.yes
 
if warning.yes
(mode 0.yes
message err_ell not ok
lookup err_ell)

end

finis
▶EOF◀