function y=prodvect(u);

% x et y sont 2 vecteurs y, leur produit vectoriel
% il y a donc 6 entrées 3 + 3

x=u(1:3);
y=u(4:6);

a1=x(2)*y(3)-x(3)*y(2);

a2=x(3)*y(1)-x(1)*y(3);

a3=x(1)*y(2)-x(2)*y(1);

y=[a1 a2 a3]';

end