%gbf.m
function[res]=gbf(i,m,s) 
res=exp(-((i-m)*(i-m)')/(2*s^2));
