clc approxerf=myownerf(0.1) exacterf=erf(0.1)