SUB draw_axes(xmin,xmax,ymin,ymax)
    ! number of tick marks
    LET ntick = 10
    ! dx distance between tick marks on x axis
    LET dx = (xmax - xmin)/ntick
    ! dy distance between tick marks on y axis
    LET dy = (ymax - ymin)/ntick
    ! include margin in window statement
    SET WINDOW xmin - dx,xmax + dx,ymin - dy,ymax + dy
    LET x0 = max(0,xmin)
    LET y0 = max(0,ymin)
    IF ymin*ymax < 0 then
       LET y0 = 0
    ELSE
       LET y0 = ymin
    END IF
    PLOT LINES: xmin,y0;xmax,y0   ! horizontal axis
    PLOT LINES: x0,ymin;x0,ymax   ! vertical axis
    LET tx = 0.1*dy               ! size of tick mark
    LET ty = 0.1*dx
    FOR itick = 0 to ntick
        LET x = xmin + itick*dx
        PLOT LINES: x,y0 - tx; x,y0 + tx    ! draw ticks on x axis
        LET y = ymin + itick*dy
        PLOT LINES: x0 - ty,y; x0 + ty,y    ! draw ticks on y axis
    NEXT itick
END SUB