t=pi/2;
p=pi/2;
a=sin(t)*cos(p);
b=sin(t)*sin(p);
c=cos(t);

figure(1)
[fx,fy,fz]=hopfmap(a,b,c);
plot3(fx,fy,fz,'b','LineWidth',2)
axis equal
hold on 
title('Hopf Links for \theta = \pi/2 and \phi=\pm \pi/2')
t=pi/2;
p=3*pi/2;
a=sin(t)*cos(p);
b=sin(t)*sin(p);
c=cos(t);
[fx,fy,fz]=hopfmap(a,b,c);
plot3(fx,fy,fz,'r','LineWidth',2)

t=pi/2;
for k=0:20
    p=pi/2*(1+k/10);
    a=sin(t)*cos(p);
    b=sin(t)*sin(p);
    c=cos(t);
    [fx,fy,fz]=hopfmap(a,b,c);
    plot3(fx,fy,fz,'k','LineWidth',1)
end
hold off

t=0;
p=0;
a=sin(t)*cos(p);
b=sin(t)*sin(p);
c=cos(t);

figure(2)
[fx,fy,fz]=hopfmap(a,b,c);
plot3(fx,fy,fz,'b','LineWidth',2)
axis equal
hold on 


t=pi;
p=0;
a=sin(t)*cos(p);
b=sin(t)*sin(p);
c=cos(t);
[fx,fy,fz]=hopfmap(a,b,c);
plot3(fx,fy,fz,'r','LineWidth',2)

p=0;
for k=0:20
    t=pi*k/20;
    a=sin(t)*cos(p);
    b=sin(t)*sin(p);
    c=cos(t);
    [fx,fy,fz]=hopfmap(a,b,c);
    plot3(fx,fy,fz,'k','LineWidth',1)
end
hold off
title('Hopf Links for \phi = 0 and \theta=0, \pi')

figure(3)
t=pi/6;
for k=0:40
    p=pi*k/20;
    a=sin(t)*cos(p);
    b=sin(t)*sin(p);
    c=cos(t);
    [fx,fy,fz]=hopfmap(a,b,c);
    plot3(fx,fy,fz,'k','LineWidth',1)
    hold on
end
hold off
title('Hopf Fibration for \theta = \pi/6 and \phi=0, 2\pi')


figure(4)
cc=['k','b','r','y'];

for k=1:4
    t=k*pi/5;
for j=1:200;
    p=j*pi/25;
    a=sin(t)*cos(p);
    b=sin(t)*sin(p);
    c=cos(t);
    [fx,fy,fz]=hopfmap(a,b,c);
     plot3(fx,fy,fz,cc(k),'LineWidth',1)
    hold on
end
end
hold off
axis equal
title('Hopf Fibration for \theta = k\pi/5')

