Расчет иллюстрации Лоренцева сокращения в программе Maple.

> restart: with(plots):

> Lorenz[x] := ((x,t,V) -> (x+V*t)/ sqrt( 1 - V^2/c^2 ));

> Lorenz[t] := ((x,t,V) -> (t+V*x/c^2)/ sqrt( 1 - V^2/c^2 ));

> x1[stay] := 0;

> t1[stay] := alpha;

> x1[move] := Lorenz[x](x1[stay],t1[stay],V);

> t1[move] := Lorenz[t](x1[stay],t1[stay],V);

> x2[stay] := 500000;

> t2[stay] := alpha;

> x2[move] := Lorenz[x](x2[stay],t2[stay],V);

> t2[move] := Lorenz[t](x2[stay],t2[stay],V);

> alpha[max]:=3;

> c:=300000;

> V:=100000;

> Point1[stay] := plot([x1[stay],t1[stay],alpha=-alpha[max]..alpha[max]],-c*alpha[max]..c*alpha[max],-alpha[max]..alpha[max],color=red,axes=BOXED):

> Point1[move] := plot([x1[move],t1[move],alpha=-alpha[max]..alpha[max]],-c*alpha[max]..c*alpha[max],-alpha[max]..alpha[max],color=red,axes=BOXED):

> Point2[stay] := plot([x2[stay],t2[stay],alpha=-alpha[max]..alpha[max]],-c*alpha[max]..c*alpha[max],-alpha[max]..alpha[max],color=blue,axes=BOXED):

> Point2[move] := plot([x2[move],t2[move],alpha=-alpha[max]..alpha[max]],-c*alpha[max]..c*alpha[max],-alpha[max]..alpha[max],color=blue,axes=BOXED):

> Cone1:=plot(x/c,x=-c*alpha[max]..c*alpha[max],color=yellow):

> Cone2:=plot(-x/c,x=-c*alpha[max]..c*alpha[max],color=yellow):

> ZeroLine:=plot([alpha,0,alpha=-c*alpha[max]..c*alpha[max]],axes=BOXED,color=green):

> display(Point1[stay], Point1[move], Point2[stay], Point2[move], Cone1, Cone2,ZeroLine);