function [v1,Pw] = proyec_ortog(A,v,G)
%Función que calcula la proyección ortogonal de v sobre W en un espacio
%euclídeo con un producto escalar con matriz de Gram G.
%   Ultano Kindelán Bustelo 09/2022

Pw=A*(A'*G*A)^(-1)*A'*G;
v1=Pw*v;
end