最佳平方逼近的Matlab

最佳平方逼近是先求积得到希尔伯特矩阵，再根据矩阵求解方程。

function [A,B,p]=zjpfbj_zzy(weightfun,primaryfun,fun,a,b)
%      ---wfun：权函数
%      ---phifun：基函数
%      ---fun：逼近函数
%      ---a,b：逼近区间的端点
%      ---n：最佳逼近的次数，默认值为1
% 输出参数：
%      ---A：法方程组的系数矩阵
%      ---b：法方程组的右端向量
%      ---p：最佳平方逼近的系数
n=size(primaryfun,1);
A=zeros(n);
for i=1:n
    for j=1:n
        f1=inline(['(',weightfun,