function res=primigauss(a,b)
    res=integrate('exp(0.5*x^2)','x',a,b)
endfunction