load (draw)$

aphi : 45;    /* Winkel in grad */
ar : 3;
ax : ar * cos (aphi * 2*%pi/360);
ay : ar * sin (aphi * 2*%pi/360);
bphi : 14;    /* Winkel in grad */
br : 0.6;
bx : br * cos (bphi * 2*%pi/360);
by : br * sin (bphi * 2*%pi/360);

draw (
   gr2d (
     user_preamble = ["set noborder"],
     xtics = false,
     ytics = false,

     points_joined=true,
     point_type = none,
     line_width = 1,
     color = black,
     points ([[-4,0], [4,0]]),
     points ([[0,-4], [0,4]]),

     color = black,
     line_width = 2,
     transparent = true,

     ellipse (0,0,0.7,0.7,0,45),
     label_alignment = right,
     label (["45", bx, by]),
     ellipse (bx+0.01,by+0.10,0.03,0.03,0,360),
     label (["180", -4.3, 0]),
     ellipse (-4.3+0.01,0+0.10,0.03,0.03,0,360),
     label (["0", 4.5, 0]),
     ellipse (4.5+0.01,0+0.10,0.03,0.03,0,360),
     label_alignment = center,
     label (["90", 0, 4.3]),
     ellipse (0.15+0.01,4.3+0.10,0.03,0.03,0,360),
     label (
          ["1", 1.15, -0.15],
          ["2", 2.15, -0.15],
          ["3", 3.15, -0.15],
          ["4", 4.15, -0.15]
     ),

     label_alignment = left,
     ellipse (0,0,1,1,0,360),
     ellipse (0,0,2,2,0,360),
     ellipse (0,0,3,3,0,360),
     ellipse (0,0,4,4,0,360),
     point_size = 2,
     point_type = 7,
     points ([[0,0], [ax,ay]]),
     label (["(3|45)", ax + 0.2, ay + 0.0]),
     xrange = [-5, 5],
     yrange = [-5, 5]
   ),

   file_name = "koordpolar",
   dimensions = [600,600],
   terminal = png
);