Add xyfit plot functionality

Started to add more functionality for plotting xy data and models, so it
comes closer to the plot_data routine (also with its qualifiers) and displays
residuals of data with x- and y-error better.