x1=-5:1.1:5;
x2=-5:1.1:5;
[X1,X2]=meshgrid(x1,x2);
X=[];
[m,n]=size(X1);
for i=1:m;
  for j=1:n;
      X=[X;X1(m,n),X2(m,n)];
  end
end
Y=X(:,1).^2+X(:,2).^2;
Yd=X1.^2+X2.^2;

global p1
p1=10;
ker='rbf';
[nsv, beta, bias] = svr(X,Y,ker);
%svrplot(X,Y,ker,beta,bias)
mx=size(X,1);
Xt=X;
Yt = svroutput(X,Xt,ker,beta,bias);

k=size(Yt,1);
for i=m:-1:1
    for j=n:-1:1
        Ys(i,j)=Yt(k);
        k=k-1
    end
end
figure(2)
mesh(X1,X2,Yd)

figure(3)
mesh(X1,X2,Ys);

figure(4)
mesh(X1,X2,Ys-Yd);