function V = gspoly(A,x1,x2)
% Gram-Schmidt ορθοκανονικοποίηση βάσης πολυωνύμων που δίνεται 
% στο διάνυσμα A. Η ορθοκανονικοποιημένη βάση 
% δίνεται στο διάνυσμα V.
[m,n] = size(A);
for k=1:n
    V(k) = A(k);
for j=1:k-1
    R(j,k) = int(V(j)*A(k),x1,x2);
    V(k) = V(k) - R(j,k)*V(j);
end
    R(k,k) = (int(V(k)*V(k),x1,x2))^(1/2);
    V(k) = V(k)/R(k,k);
end