program plplotting
    use plplot
    implicit none

    call plinit ()
    call plenv(-2.1, 2.1, -2.1, 2.1, 0, 0)
    call plend
end program