clear all
close all

epsilon=5;

% [X,Y]=ndgrid(-1:0.01:1);
% XY=epsilon*[X(:) Y(:)];
% NXY= sqrt(sum(XY.^2,2));

X=linspace(-1,1,200);
NXY=abs(X);
Z=exp(-NXY); % Matérn-gaussiana
Z=(1+NXY).*exp(-NXY); % Matérn-2

plot(X,Z)
% n=max(size(X));
% Zj=reshape(Z,n,n);
% H=surf(X,Y,Zj);
% %axis square
% set(H,'FaceColor','interp','EdgeColor','none')