B žµs^¤ã@sHddlmZddlZddlZddlmZddde  ej ¡fdd„Z dS)éNi'éc Cs*dd„}t |¡}d}xXt|ƒD]L}|t ¡}|t ¡}|||ƒkrP|d7}||||d} | ||<q Wtjdd„t|ƒDƒ|ddt|ƒd d tjd|gt |¡t |¡gd t|ƒd d t  d ¡t  d¡t  ¡t  t |¡|dt |¡|d¡t  dt|ƒd¡dS)NcSs dt tj¡t |d ¡S)Nr)ÚnpÚsqrtÚpiÚexp)Úx©rúC/Users/vegardjervell/Documents/6_semester/matmod/matmod10/mc_erf.pyÚerf_funcszplot..erf_funcrécSsg|] }|d‘qS)r r)Ú.0Úirrr ú szplot..ÚgreenzCalculated value of erf(ú))ÚcolorÚlabelzExact value of erf()rzValue of errorfunctionzNumber of trialséz.) using Monte Carlo method vs. the exact value)rÚzerosÚrangeÚrandomÚpltÚplotÚstrÚssÚerfÚylabelÚxlabelÚlegendÚylimÚtitle) ÚNÚbÚy_maxr ÚIntegralZN_Ir rÚyÚIrrr rs"     ,0  (r) Úmatplotlib.pyplotÚpyplotrÚnumpyrrZ scipy.specialÚspecialrrrrrrrr Ús