As for the plotting part, there is a points command that will do what you want. You'll want to zip up the points first, given its syntax.
Pts = zip(X,Y) # this should make a list of pairs from X and Y