829 
ANNEX 1 
Program Iglospot(input, output); 
const 
PI - 3.141592717; 
var 
minvis, maxvis, minshad, maxshad : integer; 
sunang, sundir, lookang, lookdir, delta : real; 
idx, i, j : integer; 
lr : char; 
buf : array [0..359] of char; 
function max(a, b : real) : real; 
begin if a > b then max :=■ a else max := b end; 
function deg(a : real) : real; 
begin deg 180 * a / PI end; 
function rad(a : real) : real; 
begin rad pi * a / 180 end; 
function acos(x : real) : real; 
begin 
if x > 0 then 
acos arctan(sqrt(l-x*x)/x) 
else if x < 0 then 
acos PI - acos(-x) 
else 
acos PI / 2 
end; 
function acosd(x : real) : real; 
begin acosd deg(acos(x)) end; 
function tand(x : real) : real; 
begin tand :* sin(rad(x))/cos(rad(x)) end; 
procedure prthdr(var f: text); 
begin 
writeln(f); 
writeln(f, 
' E S W N f ) ; 
writeln(f, 
' 1 2 3 0 r ); 
writeln(f, 
' 901234567890123456789012345012345678'); 
writeln(f, 
' ' ); 
end;