load('dataset10324.mat');

% beta_base = 0.05;
% beta_base = 0.1;
beta = 1.5;
func = @(r) sqrt(1+(beta*r)^2);
% func = @(r) exp(-1.*((beta*r).^2));

x_min = -2;
x_max = 2;
y_min = -3;
y_max = 1;
e_mean_max = 0.1;

phi = zeros(Nt, Nt);

% for b=0:10
% beta = 0.5 + beta_base * b
% func = @(r) exp(-1.*((beta*r).^2));
% func = @(r) sqrt(1+(beta*r)^2);



tic
for i=1:Nt
    for j=i:Nt
        temp = func(norm([xt(i,1), yt(i,1)] - [xt(j,1), yt(j,1)]));
        phi(i,j) = temp;
        phi(j,i) = temp;
    end
end
toc

tic
for i=1:Nt
    for j=1:Nt
        phi(i,j) = func(norm([xt(i,1), yt(i,1)] - [xt(j,1), yt(j,1)]));
    end
end
toc
% phi = phi'+phi;
a = phi\ft;


error = zeros(Nv,1);
f = zeros(Nv, 1);



for i=1:Nv
    for k=1:Nt
        f(i,1) = f(i,1) + a(k,1)*func(norm([xv(i,1),yv(i,1)] - [xt(k,1), yt(k,1)]));
    end
    error(i,1) = abs(fv(i,1) - f(i,1));
end
error_max = max(error)
error_mean = sum(error)/Nv
% end


n = 50;
[X,Y] = meshgrid(linspace(x_min, x_max, n), linspace(y_min, y_max, n));
Z = zeros(n, n);
for i=1:n
    for j=1:n
        for k=1:Nt
            Z(j,i) = Z(j,i) + a(k,1)*func(norm([X(1,i), Y(j,1)] - [xt(k,1), yt(k,1)]));
        end
    end
end

surf(X, Y, Z);