O15-059
> restart;
> with(plots):with(linalg):
> s1:=x^2+y^2+z^2-r^2;

                               2    2    2
                        s1 := x  + y  + z  - 1

> s2:=x^2+y^2-r*x;

                                 2    2
                          s2 := x  + y  - x

> r:=1:implicitplot3d({s1,s2},x=-1..1,y=-1..1,
> z=-2..2);

> # Symétries : yOz et xOy
> n1:=vector(3,[diff(s1,x),diff(s1,y),diff(s1,z)]);n2:=vector(3,[diff(s2,x),diff(s2,y),diff(s2,z)]);# Normales aux 2 surfaces

                        n1 := [2 x, 2 y, 2 z]


                       n2 := [2 x - 1, 2 y, 0]

> t:=simplify(crossprod(n1,n2));# Point singulier :t=0 en un point de la courbe : (R,0,0)

                  t := [-4 z y, 2 z (2 x - 1), 2 y]

> # Paramétrage de (C) : (R*(cos(t))^2,R*cos(t)*sin(t),eps*R*sin(t))