Ñò
¥3sNc        	   @   s  d  d d „  ƒ  YZ  e d j o÷ d d k Z e i d ƒ d d k Te  ƒ  Z d Ge i Gd Ge i Gd	 Ge i GHx( e	 e i d
 ƒ D] Z
 e i e
 ƒ q‚ We e i e i d d ƒe e i e i d d ƒe e i e i d d ƒe d ƒ e d ƒ e d ƒ e e ƒ e ƒ  n d S(   t	   Excitablec           B   s>   e  Z d  „  Z d „  Z d „  Z d „  Z d „  Z d „  Z RS(   c         C   s»   d |  _  d |  _ t |  i |  i  ƒ |  _ |  i d g |  _ |  i d g |  _ |  i d g |  _ |  i d g |  _ |  i d g |  _ d |  _	 d |  _
 d |  _ d |  _ d |  _ d  S(	   Ng{®Gáz„?i    i    g        g{®Gáz´?gffffffæ?gš™™™™™é?g      ð?(   t   timestept   maxtimet   intt
   meshpointst   vt   wt   wdispt   Iextt   timet   dampt   at   bt   ct   tau(   t   self(    (    s'   Parameters/LayerHysteresis/Excitable.pyt   __init__   s    				 	 	 c         C   s&   |  i  | d |  i | | d | S(   Ni   i   (   R   R
   (   R   R   R   R   (    (    s'   Parameters/LayerHysteresis/Excitable.pyt   dvdtstep   s    c         C   s(   |  i  |  i | |  i | |  i | ƒ S(   N(   R   R   R   R   (   R   t   point(    (    s'   Parameters/LayerHysteresis/Excitable.pyt   dvdt   s    c         C   s+   |  i  |  i | |  i |  i | |  i S(   N(   R   R   R   R   R   (   R   R   R   R   (    (    s'   Parameters/LayerHysteresis/Excitable.pyt   dwdtstep   s    c         C   s(   |  i  |  i | |  i | |  i | ƒ S(   N(   R   R   R   R   (   R   R   (    (    s'   Parameters/LayerHysteresis/Excitable.pyt   dwdt   s    c         C   s¯   | d j  o | d |  i  | <n |  i | |  i |  i | d <|  i | |  i | ƒ |  i | d <|  i | |  i | ƒ |  i | d <|  i | d d |  i | d <d  S(   NiX  g-Cëâ6?i   g…ëQ¸…Ë?(   R   R	   R   R   R   R   R   R   (   R   R   (    (    s'   Parameters/LayerHysteresis/Excitable.pyt   StepForward   s     %%(   t   __name__t
   __module__R   R   R   R   R   R   (    (    (    s'   Parameters/LayerHysteresis/Excitable.pyR       s   					t   __main__iÿÿÿÿNt   MacOSX(   t   *s	   Timestep:s     Max time:s      Meshpoints:i   t	   linewidthg      ð?s   time (s)t   ODEs   Excitable system(    (   R    R   t
   matplotlibt   uset   pylabt   excitedR   R   R   t   ranget   iR   t   plotR	   R   R   R   t   xlabelt   ylabelt   titlet   gridt   Truet   show(    (    (    s'   Parameters/LayerHysteresis/Excitable.pyt   <module>   s"   
	" 



