
r\                 @   s  d  d l  Z d  d l j Z d  d l Z e j j e j j e	   Z
 e e
  d Z d Z d Z d Z d Z d Z d Z d	 Z d
 Z d Z d Z d Z d d   Z e d  Z e d  Z e d Z e d Z e d Z e d Z e d Z e d Z e e e e Z e e e e e j e d e j  e   Z! e e e e  e e j e d e j  e    e e  e e!  e j  e  Z" e j# e" e! d  Z$ e j% d  e& e"  d  Z' e$ d  e' e$ d Z( e$ d  e d e e e j d  Z) d S)    Ng?ga2U0*#?g{Gz?gMb@?gV-?g&- @g_LJ?gMbU@giUMu>g      ?g?c             C   s   t  |  d  8 } t j d d   | j   D  } t j |  } Wd  QRXt j d d   | D  } xH t t |   D]4 } t j d d   | | d d   D  | | <qy W| S)Nrc             S   s   g  |  ] } | j     q S )split).0liner   r   T/Users/vegardjervell/Documents/4_semester/termo/reduksjonspotensiale_lab/plotting.py
<listcomp>   s   	 zread_data.<locals>.<listcomp>c             S   s   g  |  ] } d   q S)Nr   )r   r   r   r   r   r       s   	 c             S   s   g  |  ] } t  |   q Sr   )float)r   xr   r   r   r   "   s   	    )opennpZarray	readlinesZ	transposerangelen)filenamefilelinesdatair   r   r   	read_data   s    "2r   zdata.txtr         i  g[rPL?gfffffq@
   )*Znumpyr   Zmatplotlib.pyplotZpyplotZpltospathdirnameabspath__file__hereprintVs_Vc0s_c0rohRAFs_V0s_Ts_ps_Ur   Zdata_matrixZ	V0_pointsZT_pointsZE_pointsZp_pointsZm_pointslogZsqrtZE_merketZx_pointsZpolyfitZcoeffZlinspacemaxZx_lineZy_lineTr   r   r   r   <module>   sF   







0$

