神の言葉 - CGIBOYがサービスやめるからってそれくらいであきらめたりはしない。

■神の言葉（単にメモですが・・・）
{{t -> -(d^2*v1*v2^2 + d^2*v1*v3^2 - v1*v2^2*x0^2 - v1*v3^2*x0^2 -
2*v2^3*x0*y0 - 2*v2*v3^2*x0*y0 + v1^3*y0^2 + 2*v1*v2^2*y0^2 +
v1*v3^2*y0^2 - 2*v2^2*v3*x0*z0 - 2*v3^3*x0*z0 + 2*v1*v2*v3*y0*z0 +
v1^3*z0^2 + v1*v2^2*z0^2 + 2*v1*v3^2*z0^2 +
v1*Sqrt[-4*d^2*(v1^2 + v2^2 + v3^2)*(d^2*v1^2 + v2^2*x0^2 + v3^2*x0^2 -
2*v1*v2*x0*y0 + v1^2*y0^2 - 2*v1*v3*x0*z0 + v1^2*z0^2) +
(d^2*(2*v1^2 + v2^2 + v3^2) + v3^2*x0^2 + v1^2*y0^2 + v3^2*y0^2 -
2*v1*v3*x0*z0 + v1^2*z0^2 - 2*v2*y0*(v1*x0 + v3*z0) +
v2^2*(x0^2 + z0^2))^2])/(2*(v1^2 + v2^2 + v3^2)*
(-(v2^2*x0) + v1*v2*y0 + v3*(-(v3*x0) + v1*z0))),
A -> -(Sqrt[(2*d^2*v1^2 + d^2*v2^2 + d^2*v3^2 + v2^2*x0^2 + v3^2*x0^2 -
2*v1*v2*x0*y0 + v1^2*y0^2 + v3^2*y0^2 - 2*v1*v3*x0*z0 -
2*v2*v3*y0*z0 + v1^2*z0^2 + v2^2*z0^2 -
Sqrt[-4*d^2*(v1^2 + v2^2 + v3^2)*(d^2*v1^2 + v2^2*x0^2 + v3^2*x0^2 -
2*v1*v2*x0*y0 + v1^2*y0^2 - 2*v1*v3*x0*z0 + v1^2*z0^2) +
(d^2*(2*v1^2 + v2^2 + v3^2) + v3^2*x0^2 + v1^2*y0^2 + v3^2*y0^2 -
2*v1*v3*x0*z0 + v1^2*z0^2 - 2*v2*y0*(v1*x0 + v3*z0) +
v2^2*(x0^2 + z0^2))^2])/(v1^2 + v2^2 + v3^2)]/Sqrt[2])},
{t -> -(d^2*v1*v2^2 + d^2*v1*v3^2 - v1*v2^2*x0^2 - v1*v3^2*x0^2 -
2*v2^3*x0*y0 - 2*v2*v3^2*x0*y0 + v1^3*y0^2 + 2*v1*v2^2*y0^2 +
v1*v3^2*y0^2 - 2*v2^2*v3*x0*z0 - 2*v3^3*x0*z0 + 2*v1*v2*v3*y0*z0 +
v1^3*z0^2 + v1*v2^2*z0^2 + 2*v1*v3^2*z0^2 +
v1*Sqrt[-4*d^2*(v1^2 + v2^2 + v3^2)*(d^2*v1^2 + v2^2*x0^2 + v3^2*x0^2 -
2*v1*v2*x0*y0 + v1^2*y0^2 - 2*v1*v3*x0*z0 + v1^2*z0^2) +
(d^2*(2*v1^2 + v2^2 + v3^2) + v3^2*x0^2 + v1^2*y0^2 + v3^2*y0^2 -
2*v1*v3*x0*z0 + v1^2*z0^2 - 2*v2*y0*(v1*x0 + v3*z0) +
v2^2*(x0^2 + z0^2))^2])/(2*(v1^2 + v2^2 + v3^2)*
(-(v2^2*x0) + v1*v2*y0 + v3*(-(v3*x0) + v1*z0))),
A -> Sqrt[(2*d^2*v1^2 + d^2*v2^2 + d^2*v3^2 + v2^2*x0^2 + v3^2*x0^2 -
2*v1*v2*x0*y0 + v1^2*y0^2 + v3^2*y0^2 - 2*v1*v3*x0*z0 -
2*v2*v3*y0*z0 + v1^2*z0^2 + v2^2*z0^2 -
Sqrt[-4*d^2*(v1^2 + v2^2 + v3^2)*(d^2*v1^2 + v2^2*x0^2 + v3^2*x0^2 -
2*v1*v2*x0*y0 + v1^2*y0^2 - 2*v1*v3*x0*z0 + v1^2*z0^2) +
(d^2*(2*v1^2 + v2^2 + v3^2) + v3^2*x0^2 + v1^2*y0^2 + v3^2*y0^2 -
2*v1*v3*x0*z0 + v1^2*z0^2 - 2*v2*y0*(v1*x0 + v3*z0) +
v2^2*(x0^2 + z0^2))^2])/(v1^2 + v2^2 + v3^2)]/Sqrt[2]},
{t -> -(d^2*v1*v2^2 + d^2*v1*v3^2 - v1*v2^2*x0^2 - v1*v3^2*x0^2 -
2*v2^3*x0*y0 - 2*v2*v3^2*x0*y0 + v1^3*y0^2 + 2*v1*v2^2*y0^2 +
v1*v3^2*y0^2 - 2*v2^2*v3*x0*z0 - 2*v3^3*x0*z0 + 2*v1*v2*v3*y0*z0 +
v1^3*z0^2 + v1*v2^2*z0^2 + 2*v1*v3^2*z0^2 -
v1*Sqrt[-4*d^2*(v1^2 + v2^2 + v3^2)*(d^2*v1^2 + v2^2*x0^2 + v3^2*x0^2 -
2*v1*v2*x0*y0 + v1^2*y0^2 - 2*v1*v3*x0*z0 + v1^2*z0^2) +
(d^2*(2*v1^2 + v2^2 + v3^2) + v3^2*x0^2 + v1^2*y0^2 + v3^2*y0^2 -
2*v1*v3*x0*z0 + v1^2*z0^2 - 2*v2*y0*(v1*x0 + v3*z0) +
v2^2*(x0^2 + z0^2))^2])/(2*(v1^2 + v2^2 + v3^2)*
(-(v2^2*x0) + v1*v2*y0 + v3*(-(v3*x0) + v1*z0))),
A -> -(Sqrt[(2*d^2*v1^2 + d^2*v2^2 + d^2*v3^2 + v2^2*x0^2 + v3^2*x0^2 -
2*v1*v2*x0*y0 + v1^2*y0^2 + v3^2*y0^2 - 2*v1*v3*x0*z0 -
2*v2*v3*y0*z0 + v1^2*z0^2 + v2^2*z0^2 +
Sqrt[-4*d^2*(v1^2 + v2^2 + v3^2)*(d^2*v1^2 + v2^2*x0^2 + v3^2*x0^2 -
2*v1*v2*x0*y0 + v1^2*y0^2 - 2*v1*v3*x0*z0 + v1^2*z0^2) +
(d^2*(2*v1^2 + v2^2 + v3^2) + v3^2*x0^2 + v1^2*y0^2 + v3^2*y0^2 -
2*v1*v3*x0*z0 + v1^2*z0^2 - 2*v2*y0*(v1*x0 + v3*z0) +
v2^2*(x0^2 + z0^2))^2])/(v1^2 + v2^2 + v3^2)]/Sqrt[2])},
{t -> -(d^2*v1*v2^2 + d^2*v1*v3^2 - v1*v2^2*x0^2 - v1*v3^2*x0^2 -
2*v2^3*x0*y0 - 2*v2*v3^2*x0*y0 + v1^3*y0^2 + 2*v1*v2^2*y0^2 +
v1*v3^2*y0^2 - 2*v2^2*v3*x0*z0 - 2*v3^3*x0*z0 + 2*v1*v2*v3*y0*z0 +
v1^3*z0^2 + v1*v2^2*z0^2 + 2*v1*v3^2*z0^2 -
v1*Sqrt[-4*d^2*(v1^2 + v2^2 + v3^2)*(d^2*v1^2 + v2^2*x0^2 + v3^2*x0^2 -
2*v1*v2*x0*y0 + v1^2*y0^2 - 2*v1*v3*x0*z0 + v1^2*z0^2) +
(d^2*(2*v1^2 + v2^2 + v3^2) + v3^2*x0^2 + v1^2*y0^2 + v3^2*y0^2 -
2*v1*v3*x0*z0 + v1^2*z0^2 - 2*v2*y0*(v1*x0 + v3*z0) +
v2^2*(x0^2 + z0^2))^2])/(2*(v1^2 + v2^2 + v3^2)*
(-(v2^2*x0) + v1*v2*y0 + v3*(-(v3*x0) + v1*z0))),
A -> Sqrt[(2*d^2*v1^2 + d^2*v2^2 + d^2*v3^2 + v2^2*x0^2 + v3^2*x0^2 -
2*v1*v2*x0*y0 + v1^2*y0^2 + v3^2*y0^2 - 2*v1*v3*x0*z0 -
2*v2*v3*y0*z0 + v1^2*z0^2 + v2^2*z0^2 +
Sqrt[-4*d^2*(v1^2 + v2^2 + v3^2)*(d^2*v1^2 + v2^2*x0^2 + v3^2*x0^2 -
2*v1*v2*x0*y0 + v1^2*y0^2 - 2*v1*v3*x0*z0 + v1^2*z0^2) +
(d^2*(2*v1^2 + v2^2 + v3^2) + v3^2*x0^2 + v1^2*y0^2 + v3^2*y0^2 -
2*v1*v3*x0*z0 + v1^2*z0^2 - 2*v2*y0*(v1*x0 + v3*z0) +
v2^2*(x0^2 + z0^2))^2])/(v1^2 + v2^2 + v3^2)]/Sqrt[2]}}

■プログラム
上の計算をするプログラムを書いた．変数が5個だからN^5回計算しなきゃいけない．並列計算させなきゃ・・・
しっかし，まだ完成していないけど，とんでもないコードになりつつある・・・もっときれいにかけたらいいなぁ・・・