% Input: P0, vettore colonna 3x1 del punto nello spazio
%        q0, vettore colonna 3x1 di un punto nel piano
%        v1 e v2, vettori colonna 3x1 delle direzioni lungo il piano  
% Output: dist, distanza di P0 dal piano
function dist = distanza_punto_piano_esplicito(P0, q0, v1, v2)
  V = [ v1 v2 ] ; 
  q = @(lambda) q0 + V*lambda ; % funzione piano

  % Metodo A
  A = V' * V ;       % matrice quadrata dei coefficienti
  b = V' * (P0-q0) ; % vettore dei termini noti
  lambda = A\b ;     % soluzione, lambda
  dist   = norm( P0 - q(lambda) ) ;
end