%Script file to show linear regression using the EDM:
X=randn(30,2);
Y=exp(-(X(:,1).^2+X(:,2).^2)./4);

A=edm(X,X);
alpha=inv(A)*Y;

Z=A*alpha;
plot3(X(:,1),X(:,2),Y,'r*',X(:,1),X(:,2),Z,'b^');
hold on

%Create a 3-d surface
[X1,Y1]=meshgrid(-3:0.125:3);
[m,n]=size(X1);
for j=1:m
    for k=1:n
        Z(j,k)=0;
        for t=1:30
            Z(j,k)=Z(j,k)+alpha(t)*norm([X1(j,k) Y1(j,k)]-X(t,:));
        end
    end
end
meshc(X1,Y1,Z);
hold off