Calculation for an object moving towards a galaxy

This page gives the details of the calculation for an object with no internal forces moving from the origin towards a comoving galaxy in an expanding universe, as described in my paper Expanding Space: Just say No at quantropy.org. The calculation is done using the Euler Math Toolbox. 
The file fxspace.e contains the function

    togalaxy(t0,t1,v,chig,w)

which calculates when the two ends of the object reach the galaxy at
χg and returns tf0,tf1,vf0,lf,lf/l

>load "fxspace.e"

Now run the calculation for five different values of w:
-1/2,-1/3,0,1/3,2/3

>togalaxy(1,1.005,0.2,0.1,-1/2)

[ 2.82753179947  2.85311215555  0.0509861813527  0.00130424467423
1.30424467423 ]

>togalaxy(1,1.005,0.2,0.1,-1/3)

[ 1.98345631632  1.99337359791  0.102372660939  0.00101525850479
1.01525850479 ]

>togalaxy(1,1.005,0.2,0.1,0)

[ 1.72215757485  1.72906425559  0.140660781554  0.00097149910982
0.97149910982 ]

>togalaxy(1,1.005,0.2,0.1,1/3)

[ 1.64521106023  1.65139549762  0.157163820363  0.000971969805636
0.971969805636 ]

>togalaxy(1,1.005,0.2,0.1,2/3)

[ 1.60800659179  1.61386610643  0.166447891354  0.000975303856894
0.975303856894 ]