load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"

begin

N = 500
x = (1.0 * ispan(1,N,1)) / N

y = new( (/ 2, N /), float )

y(0,:) = x
y(1,:) = x ^ 2

wks = gsn_open_wks("eps", "foo")
res = True
res@xyXStyle = "Log"
res@trXMaxF = 1.0
res@trXMinF = 5.0 * x(0)
res@trYMaxF = 1.0
res@trYMinF = -1.0

p = gsn_xy(wks,x,y,res)

end