User Portlet

Here is code that should be adaptable for the purpose at hand. It returns a list {total time, x intersection, y intersection, new x direction, new y direction}. Hence it is suitable for nesting. [mcode]solveBilliard[{t0_, x0_, y0_, dx_, dy_}] := ...