% ns=[1 3 5 6 10];
% for k=1:length(ns)
%     c=MLS_Fred_II_sym(K,f,fii,ns(k),lambda,a,b);
%     y=@(x) fii(x,1:ns(k))*c;
%     eps(k)=norm(double(subs(y,x,testGrid))-y_exact(testGrid))/norm(y_exact(testGrid));
% end
% plot(ns,eps);
% xlabel('n');
% ylabel('eps');