function myerfvalue=myownerf(x)
g=inline('exp(-z^2)');
myerfvalue=2/pi^0.5*trapezoidal_rule(g,0,x,1000);