
$Oc           @   s^   d  d l  m Z d  d l Td  d l Td   Z d   Z d   Z d   Z d   Z d   Z	 d	 S(
   i(   t   division(   t   *c         C   s   t  |  d  d S(   Ni   gQ63E?(   t   exp(   t   x(    (    sD   /homepages/25/d321765456/htdocs/pub_http_internet/python/normal_m.pyt	   normalpdf%   s    c         C   s   d d t  |  t d   S(   Ng      ?g      ?i   (   t   erft   sqrt(   R   (    (    sD   /homepages/25/d321765456/htdocs/pub_http_internet/python/normal_m.pyt	   normalcdf*   s    c         C   s   t  d  t d |  d  S(   Ni   i   (   R   t   inverf(   R   (    (    sD   /homepages/25/d321765456/htdocs/pub_http_internet/python/normal_m.pyt   invnorm4   s    c         C   s   d t  |   S(   Ng      ?(   t   erfc(   R   (    (    sD   /homepages/25/d321765456/htdocs/pub_http_internet/python/normal_m.pyR      s    c         C   s   d } d } d } d } d } d } d } d } d	 }	 d
 }
 d } t  |   } | d k r^ | Sd d d | } | t | d | | | | | | | | | | | | | | | | |	 | |
  } |  d k  r d | } n  | S(   Ngś?g5 ?g`yg?gƸ?gꪂIǿg#v?g9)gS?gޅ1Ogv(?g      ?g        g      ?i   i    g       @(   t   absR   (   R   t   a1t   a2t   a3t   a4t   a5t   a6t   a7t   a8t   a9t   a10t   rvt   zt   t(    (    sD   /homepages/25/d321765456/htdocs/pub_http_internet/python/normal_m.pyR
      s(    [c         C   sw  d } d } d } t  |   | k r, | |  St  |   d k rg t j d I|  Id IJt j d  n  | t  |   } t d |  } d	 | } x t d
 |  D] } d t |  } t  |   | }	 t  |	  | k  r |  d
 k  r | S| Sn  | | }
 | |	 |
 9} | } | | 7} t  | |  | k  r |  d
 k  r@| S| Sq q Wt j d t |   d IJt j d  d  S(   Ni2   g+=gjﴑ[?g      ?s   inverf domain error: xs   not in (-1, 1).i   g?g?i    s   inverf(s   ): did not converge.(   R   t   syst   stderrt   exitR   t   rangeR
   t   str(   R   t   maxitst   tolt   sqrtpio2t   erfit   y0t   derfit   itert   y1t   dy1t   dy0(    (    sD   /homepages/25/d321765456/htdocs/pub_http_internet/python/normal_m.pyR      s8    


N(
   t
   __future__R    t   matht   kerlutilR   R   R	   R   R
   R   (    (    (    sD   /homepages/25/d321765456/htdocs/pub_http_internet/python/normal_m.pyt   <module>    s   

		
			"