begin x = ispan(-50,50,1) ; 101 points (-50,-49,...,49,50) y = sin(0.0628*x) wks = gsn_open_wks("png","xy1d") ; Open a png file called "xy1d.png". res = True ; Create a variable to hold plot resources. res@xyLineColor = "Blue" ; Line color (case insensitive) res@xyLineThicknessF = 5 ; 5x as thick as default res@trYMinF = min(y) ; Set axes limits for res@trYMaxF = max(y) ; both Y and res@trXMinF = min(x) ; X axes. res@trXMaxF = max(x) plot = gsn_csm_xy(wks,x,y,res) ; Call the gsn_csm function for drawing ; the curve. end