CREATE OR REPLACE FUNCTION azymut3(pkt1 geometry, pkt2 geometry) RETURNS FLOAT LANGUAGE plpgsql AS $$
	DECLARE

	azym FLOAT;
	pkt1x FLOAT;
	pkt2x FLOAT;
	pkt1y FLOAT;
	pkt2y FLOAT;
   BEGIN

	pkt1x :=ST_X(pkt1);
	pkt2x :=ST_X(pkt2);
	pkt1y :=ST_Y(pkt1);
	pkt2y :=ST_Y(pkt2);
	azym := ATN2(COS(pkt1x) * SIN(pkt2x) - SIN(pkt1x) * COS(pkt2x)* COS(pkt2y - pkt1y), SIN(pkt2y - pkt1y) * COS(pkt2x));

	 RETURN @azym ;
   END
   $$ ;
    SELECT azymut3(ST_GeomFromEWKT('POINT(1 2)'), ST_GeomFromEWKT('POINT(3 4)'));