B
    ]\U^                 @   s4   d dl Zd dlZd dlZG dd dZd	ddZdS )
    Nc               @   s   e Zd Zd	ddZdd ZdS )
Temp      ?X  2   c             C   sN   |t j | | _| j| | _| j|t j |  | _|| _|| _|| _|| _d S )N)	cT0t_iso_startt_iso_sluttt_endT_isoHRCRt_iso)selfr   r   r   r    r   A/Users/vegardjervell/Documents/6_semester/matmod/matmod7/matte.py__init__   s    zTemp.__init__c             C   sZ   || j k rtj| j|  | _n6|| jk r2| j| _n"|| jkrT| j| j|| j   | _| jS )N)	r   r   r   r   tempr	   r   r
   r   )r   tr   r   r   __call__   s    



zTemp.__call__N)r   r   r   r   )__name__
__module____qualname__r   r   r   r   r   r   r      s   
r         ?X    c                s   t | |||d  fdd}tj|d jd }|tj| j jd  }|tj| j jd  }tjd|  | d|  d }	tjd|  | d|  d }
tjd|  | d|  d }|	|
|gS )N)r   r   r   r   c                s(   t jt j tt j t j |    S )N)r   ZgammaZM0npZexpQZR_gas)r   )r   r   r   <lambda>(   s    zR.<locals>.<lambda>r      g    .A)	r   scipyZ	integrateZquadr   r	   r
   r   ZR0)r   r   r   r   nZ	integrandZintegral_heatingZintegral_isoZintegral_totalZ	R_heatingZR_isoZR_totalr   )r   r   R#   s    r"   )r   r   r   r   r   )numpyr   Zscipy.integrater    Z	constantsr   r   r"   r   r   r   r   <module>   s   