2, 3]];
                var Bx = [p1x, ds*d1y, p2x, ds*d2y]
                var Px = solved_linear(Ax, Bx, 4);
                var Ay = [
                    [1, 0, 0, 0],
                    [0, 1, 0, 0],
                    [1, 1, 1, 1],
                    [0, 1, 2, 3]];
                var By = [p1y, (-ds)*d1x, p2y, (-ds)*d2x]
                var Py = solved_linear(Ay, By, 4);
                for(var k = 0; k < grid_size; k += 1) {
                    var t1 = (k + 0.01) / grid_size;
                   