function hypercube(u,v)
  % Skapar en ON-bas u0, v0 för planet som u,v spänner upp
  u0= u/norm(u,2);
  v0= v-(v'*u0)*u0;
  v0=v0/norm(v0,2);
  % Skapar matrisen P för 4D till 2D projektion
  P=[1;0]*u0'+ [0;1]*v0'
  clf
  hold on
  % Itererar över hörn (x,y,z,w) i hyperkuben
  for x=0:1
      for y= 0:1
          for z=0:1
              for w=0:1
                  % Itererar över de 4 närliggande hörnen
                  for i=1:4
                      a=[x;y;z;w]; % hörnkoordinater
                      b=a; 
                      b(i)= 1-b(i); % grannhörnets koordinater
                      c=[P*a P*b]; % dessa två hörns 2D projektioner
                      plot(c(1,:), c(2,:))
                  end
              end
          end
      end
  end
end

