B
    ¿*™\þ  ã               @   sl   d dl Zd dlZd dlmZ ddd„Zddd„Ze	d e	d e	e	dfd	d
„Z
dd„ Ze	d e	d e	e	dfdd„ZdS )é    Né   Fc             C   s°   t  | ¡} t  |¡}tj|  tj }tj| tj }|r¤tjtj tj }tjtj tj }t  ||| d  | d | || d  | d  ¡}|||  |fS |||  S )Né   )	ÚnpÚarrayÚkÚroh_isoÚMm_isoÚroh_mbÚMm_mbÚs_VÚsqrt)ÚV_isoÚV_mbÚ
usikkerhetÚmol_isoÚmol_mbÚ	s_mol_isoÚs_mol_mbZs_x_iso© r   úW/Users/vegardjervell/Documents/4_semester/termo/tokomponent_faselikevekt_lab/regning.pyÚx_iso   s    

r   c             C   s°   t  | ¡} t  |¡}tj|  tj }tj| tj }|r¤tjtj tj }tjtj tj }t  ||| d  | d | || d  | d  ¡}|||  |fS |||  S )Nr   )	r   r   r   r	   r
   r   r   r   r   )r   r   r   r   r   r   r   Zs_x_mbr   r   r   Úx_mb   s    

r   Úisoc       	   	   C   sô  |dkrxt  tjtj d| dtj   ¡||  tj|   }t  t  tjtj d| dtj   ¡||  tj|   tj tj|d    tj d t  tjtj d| dtj   ¡||  tj|   tjtjtjd    tj d  t  tjtj d| dtj   ¡|tj|   | d  t  tjtj d| dtj   ¡||   tj|d    | d  ¡}||fS |dkrðt  tj	tj d| dtj
   ¡||  tj|   }t  t  tj	tj d| dtj
   ¡||  tj|   tj	 tj|d    tj d t  tj	tj d| dtj
   ¡||  tj|   tj	tjtj
d    tj d  t  tj	tj d| dtj
   ¡|tj|   | d  t  tj	tj d| dtj
   ¡||   tj|d    | d  ¡}||fS d S )Nr   r   r   Úmb)r   Úexpr   Ú	H_vap_isoÚRÚT_b_isoÚp0r   Ús_TZH_vap_mbZT_b_mb)	ÚyÚs_yÚxÚs_xÚTÚpÚ	komponentZkoeffr   r   r   r   Úaktivitetskoeff(   s    
6äH
6äHr'   c             C   s\   t  | ¡} t| ƒt| ƒ }| | }tt|ƒƒt|ƒ }t  t|| d ƒt|ƒ ¡}||fS )Nr   )r   r   ÚsumÚlenÚabsr   )ZpunkterZsnittZavvikZsnitt_avvikZ	std_avvikr   r   r   ÚstandardavvikD   s    
r+   c          	   C   s4  t  t  tjtj d| dtj   ¡||  tj|   tj tj|d    tj d t  tjtj d| dtj   ¡||  tj|   tjtjtjd    tj d  t  tjtj d| dtj   ¡|tj|   | d  t  tjtj d| dtj   ¡||   tj|d    | d  ¡}d S )Nr   r   )	r   r   r   r   r   r   r   r   r   )r    r!   r"   r#   r$   r%   r&   r   r   r   r   Ú usikkerhet_aktivitetskoeffisientQ   s    är,   )r   r   F)r   r   F)Únumpyr   Ú
konstanterr   Úscipy.optimizeÚoptimizeÚscor   r   Úlistr'   r+   r,   r   r   r   r   Ú<module>   s
   

