Plot travel path of earthquake to a seismological station

Hi,

I have 4 column vectors: Station latitude, station longitude, earthquake epicentre latitude and earthquake epicentre longitude.

I want to plot lines between each earthquake epicentre and the corresponding station which recorded it on a map I have produced using the worldmap function. For example, if station A detected earthquake 1, plot the line between the station ordinate and the earthquake ordinate.

I can't seem to find a function which effectively plots between two sets of coordinates on a map, so any help would be appreciated.

Thanks