eikonal_raylength : Computes the raylength from S to R using the eikonal equaiton
 
  Call:
    raylength=eikonal_raylength(x,y,v,S,R,tS,doPlot)