LstSquApprox_Hibert.m

function [fintP,axk,Snx,Sqa_err]=LstSquApprox_Hibert(f,w,n)
% Hibert空间最佳平方逼近计算程序,
syms t

   i=0:n; Funx=t.^i;
   d=int(w*f*Funx,0,1);

   Hibert=hilb(n+1);

   digits(5);
   fintP=vpa(d);
   axk=Hibert\fintP';
   Snx=Funx*axk;

   Sqa_err=abs(vpa(int(f*f,0,1))-vpa(d*axk));