function g=guassdist(x,c,b)
%計算格林矩陣
% x為np×ni的輸入矩陣
% c為學習后的ni×m中心矩陣
% b為方差
%g為np×m的隱含層輸出矩陣
[np,ni]=size(x);
[ni,m]=size(c);
q=zeros(np,m);
p=zeros(ni,m);
 if ni==1
    for i=1:np
      p=x(i,:)'*ones(1,m);
      q(i,:)=abs(p-c);
    end
  else
     for i=1:np
     p=x(i,:)'*ones(1,m);
     q(i,:)=sum((p-c).^2);
     end
  end
  g=exp(-q/(b*b));