B
    SûB^ù  ã               @   sH   d dl mZ d dlZddd„Zdd„ Zeddd	d
fdd„Zddd„ZdS )é    Nçš™™™™™¹?é   c             C   s   | ||  S )N© )ÚtÚTÚrÚT_sr   r   úN/Users/vegardjervell/Documents/6_semester/matmod/matmod4/diff_cooling_paths.pyÚdT_dt   s    r
   c             C   s€   || |ƒ}|| |d  ||| d  ƒ}|| |d  ||| d  ƒ}|| | |||  ƒ}||d |d|  d|  |   S )Né   é   r   )ÚxÚyÚhZdy_dxÚk1Úk2Úk3Zk4r   r   r	   Úrk4_step_func	   s
    
r   éZ   éA   é   c             C   sÌ   |}d}|g}|g}x:||| krN| |||t ƒ}| |¡ ||7 }| |¡ qW ||8 }| |¡ | |¡ || }	d}
||	g}ddg}x6|	|kr¾| |
|	|t ƒ}	| |	¡ |
|7 }
| |
¡ qŠW ||||fS )Nr   )r
   Úappend)Znum_step_funcÚT_0ÚdtÚT_coldÚcream_coolingZlate_TZlate_tÚlate_T_listÚlate_t_listZearly_TZearly_tÚearly_T_listÚearly_t_listr   r   r	   Úoptimum_cooling_time   s,    




r    éP   ç{®Gáz”?c             C   s\   t | |||d\}}}}tj||ddd tj||ddd t d¡ t d¡ t ¡  d S )	N)r   r   r   r   ÚbluezEarly cream)ÚcolorÚlabelÚredz
Late creamz
Time [min]u
   Temp [Â°C])r    ÚpltÚplotÚxlabelÚylabelÚlegend)r   r   r   r   r   r   r   r   r   r   r	   Úplot_optimum_cooling8   s    

r,   )r   r   )r!   r"   r   r   )	Úmatplotlib.pyplotÚpyplotr'   ÚnumpyÚnpr
   r   r    r,   r   r   r   r	   Ú<module>   s
   
#