% x + y + z + 1 = 0 => c=[1;1;1], d=1
c = [ 1 ; 1 ; 1 ] ;
d = 1 ;
[q0,v1,v2] = piano_esplicito(c,d) ;

[c1,d1] = piano_canonico(q0,v1,v2) ;



q0 = [ 1 ; 2 ; 3 ] ;
v1 = [ 1 ; 0 ; 0 ] ; 
v2 = [ 0 ; 1 ; 0 ] ; 
[c,d] = piano_canonico(q0,v1,v2) ;
[q0,v1,v2] = piano_esplicito(c,d) ;