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(i,j),X2(i,j)];
  end
end
Yd=X1+X2;
Y=X(:,1)+X(:,2);

global p1
p1=1.2;
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);