Проекция кривой, заданной параметрически, на внутреннюю поверхность сферы

restart;
x := -10*cos(t)-2*cos(5*t)+15*sin(2*t);
y := -15*cos(2*t)+10*sin(t)-2*sin(5*t);
z := 10*cos(3*t);
r := 30; #радиус сферы
n := 1080; #количество шагов
h := (2*3.141592654)/n; #ширина шага

Lx := [seq(xx[i], i = 1 .. n)];
Ly := [seq(yy[i], i = 1 .. n)];
Lz := [seq(zz[i], i = 1 .. n)];
Ld := [seq(d[i], i = 1 .. n)];

#находим координаты точек
tt := 0;
for i to n+1 do
xx[i] := -10*cos(tt)-2*cos(5*tt)+15*sin(2*tt);
yy[i] := -15*cos(2*tt)+10*sin(tt)-2*sin(5*tt);
zz[i] := 10*cos(3*tt);
tt := tt+h
end do;

# находим расстояния точек от начала координат [0,0,0]
for i to n do
d[i] := sqrt((xx[i]+0)^2+(yy[i]+0)^2+(zz[i]+0)^2)
end do;

# проецируем точки на сферу
for i to n do
xx[i] := xx[i]*r/d[i];
yy[i] := yy[i]*r/d[i];
zz[i] := zz[i]*r/d[i]
end do;

# рисуем
with(plottools);
with(plots);
Sphere := plot3d([r*sin(theta)*cos(phi), r*sin(theta)*sin(phi), r*cos(theta)], theta = 0 .. Pi, phi = 0 .. 2*Pi, style = wireframe);
Curve := spacecurve([x, y, z], t = 0 .. 2*Pi, color = green);
points := point(Lx, Ly, Lz, symbol = circle, color = black, symbolsize = 15);

plots[display](points, Curve, Sphere, axes = box);