a=0; b=3.09;sigma=1;x0=0;nn=11
gau(a,b,sigma,x0,nn)
disp(p)