function result=gauss_f(n,a)

for i=1:2*n+1
    
      b=i-(n+1);
       
       result(i)=exp(-(b*b)/(2*a*a))/sqrt(2*pi*a*a);
 
   
end
cc=0;

for i=1:2*n+1
    
      cc=cc+result(i );
      
 
end   
 result=result/sqrt(cc);
%result=result/sqrt(2*n+1);