% Valori particolari
P0 = [ 3 ; 3 ; 5 ] ;

q0 = [ 0 ; 4 ; 2 ] ;
v1 = [ 1 ; 3 ; 0 ] ;
v2 = [ 0 ; 2 ; 5 ] ;

V = [ v1 v2 ] ;               % matrice 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_A = norm( P0 - q(lambda) )

% Metodo B
c = cross( v1 , v2 ) ; % prodotto vettoriale, vettore ortognale al piano
e_c = c / norm( c ) ;  % normalizzazione (lunghezza unitaria)
dist_B = abs( dot(P0-q0, e_c) )