clc;clear all;

dataMatrix=[1;2;4;5;6];y=[1;1;-1;-1;1];C=50;

[n,d]=size(dataMatrix);

alphas=zeros(n,1);K=zeros(n,n);

for i=1:n

   for j=1:n   K(i,j)=  exp(-  norm(  dataMatrix(i,:)-dataMatrix(j,:) ,2)/(2*2)  )  ;H(i,j)=K(i,j).*y(i).*y(j); end

end

c = -ones(n,1);

Aeq = y.';

beq = 0;

lb = zeros(n,1);

ub = ones(n,1)*C;

alphas0 = zeros(n,1);

alphas_final=quadprog(H,c,[],[],Aeq,beq,lb,ub,alphas0);

for k=1:n

   if(abs(  alphas_final(k)  )>0.00001) break;end

end

sum=0;

for g=1:n  sum=sum+K(g,k).*y(g).*alphas_final(g);end

b_final=1/y(k)-sum;  dataTest=dataMatrix;

for i=1:n

   sum=0;

   for g=1:n K(g,i)=exp(-  norm(  dataMatrix(g,:)-dataTest(i,:) ,2)/(2*2)  )  ;sum=sum+K(g,i).*y(g).*alphas_final(g);end

  if ((sum+b_final)>0) z(i)=1; else z(i)=-1; end;

end