B
    ùB^  ã            	   @   sR   d dl Zd dlmZ d dlZddd„Zddd	„Zddd„Z	ddd„Z
ddd„ZdS )é    Nçš™™™™™¹?é   éS   é2   c             C   s2   t  d|d¡}||| t  |  | ¡  }||fS )Nr   gš™™™™™¹?)ÚnpÚarangeÚexp)ÚrÚT_sÚT_0Út_fÚt_0ÚtidÚtemp© r   úJ/Users/vegardjervell/Documents/6_semester/matmod/matmod4/analytic_Euler.pyÚ
Analytical   s    r   c             C   sH   t | ||||d\}}tj||ddd t ¡  t d¡ t d¡ d S )N)r	   r
   r   r   r   r	   ZAnalytic)ÚlabelzTime, [min]u   Temperature, [Â°C])r   ÚpltÚplotÚlegendÚxlabelÚylabel)r	   r
   r   r   r   Úx_pointsÚy_pointsr   r   r   Úplot_Analytical   s
    
r   é
   c             C   s€   t || |  ƒ}t |¡}||d< ||fdd„}	x8tt|ƒd ƒD ]$}
||
 | |	||
 ƒ  ||
d < qBW t |||¡}||fS )Nr   c             S   s   | | |  S )Nr   )ÚTr
   r	   r   r   r   Úf   s    zEuler.<locals>.fé   )Úintr   ÚzerosÚrangeÚlenÚlinspace)Úhr	   r
   r   ÚT_fr   r   ÚNZT_pointsr   ÚnZt_pointsr   r   r   ÚEuler   s    
$r)   c       	   	   C   sL   t | ||||||d\}}tj||ddd t ¡  t d¡ t d¡ d S )N)r	   r
   r   r&   r   r   ÚgZEulers)r   zTime, [min]u   Temperature, [Â°C])r)   r   r   r   r   r   )	r%   r	   r
   r   r&   r   r   r   r   r   r   r   Ú
plot_Euler$   s
    
r+   FTc
          	   C   s`   |rt ||||dd |	r0t| ||||||d |rPt ¡ \}
}tj|
|dd t ||¡ d S )Nr   )r	   r
   r   r   r   )r	   r
   r   r&   r   r   zExperimental data)r   )r   r+   ÚflÚles_real_cupr   ÚscatterZxlim)r%   r	   r
   r   r&   r   r   Zshow_experimental_dataZshow_analyticalZshow_numericr   r   r   r   r   Úplot_Analytical_Euler,   s    r/   )r   r   r   r   r   )r   r   r   r   r   )r   r   r   r   r   r   )r   r   r   r   r   r   )	r   r   r   r   r   r   FTT)Únumpyr   Úmatplotlib.pyplotÚpyplotr   Ú	fillesingr,   r   r   r)   r+   r/   r   r   r   r   Ú<module>   s   



  