a
    kº”h*•  ã                   @   sð  d dl mZmZmZmZmZmZmZ d dlm	Z	m
Z
 d dlmZmZmZmZ d dlmZ d dlmZmZmZ d dlmZ d dlmZ d dlmZmZ d d	lmZ d d
lm Z m!Z! d dl"m#Z#m$Z$m%Z% d dl&m'Z'm(Z(m)Z)m*Z* d dl+m,Z, d dl-m.Z. d dl/m0Z0m1Z1m2Z2 d dl3m4Z4 d dl5m6Z6m7Z7m8Z8 d dl9m:Z:m;Z; d dl<m=Z= d dl>m?Z? d dl@mAZAmBZBmCZCmDZDmEZEmFZFmGZGmHZH d dlImJZJmKZKmLZLmMZM d dlNmOZOmPZP d dlQmRZRmSZS d dlTmUZU eOdd„ ƒZVeOdd„ ƒZWeOdd„ ƒZXeOd d!„ ƒZYd"S )#é    )Úlaplace_transformÚinverse_laplace_transformÚLaplaceTransformÚInverseLaplaceTransformÚ_laplace_deep_collectÚlaplace_correspondenceÚlaplace_initial_conds)ÚFunctionÚ
expand_mul)Ú
EulerGammaÚSubsÚ
DerivativeÚdiff)Úfactor_terms)ÚIÚooÚpi)ÚEq)ÚS)ÚSymbolÚsymbols©Úsimplify)ÚAbsÚre)ÚexpÚlogÚ	exp_polar)ÚcoshÚsinhÚcothÚasinh)Úsqrt)Ú	Piecewise)ÚatanÚcosÚsin)ÚAnd)Ú
lowergammaÚgammaÚ
uppergamma)Ú
DiracDeltaÚ	Heaviside)ÚSingularityFunction)Úlerchphi)ÚfresnelcÚfresnelsÚerfÚerfcÚEiÚCiÚexpintÚE1)ÚbesseliÚbesseljÚbesselkÚbessely)ÚslowÚwarns_deprecated_sympy)ÚMatrixÚeye)Úsc               
   C   sî4  t } t}tddd\}}}tdddd}tdƒ\}}}tdƒ}	tdƒ}
td	ƒ}td
ƒ}tdƒ}tt|| ||  ƒtdt|| ||  |d  ƒƒ |ƒt|| |d  |||   ƒtdt|| |||   ƒƒ ksÜJ ‚t t||ƒ|dƒ|tdd}t	|||iƒ}t
|||g d¢iƒ}|td |tƒ dtd   dt  d ksLJ ‚dtt|ƒ|tƒ d  ¡ ddtd   ks~J ‚| || |d  |tdƒd   |tƒ|td  dtd   dttƒ dttdƒd     ddfksêJ ‚| |||  |tƒ| t| t ƒ t| t ƒ ddfks,J ‚| dt|| ƒ |tƒttƒtdt ƒ t|t ƒ tt|ƒttƒ ƒ ddfks„J ‚| t|ƒ||  |tƒt t|ƒ t|t ƒ tt|ƒttƒ ƒ ttƒtdt ƒ  ddfksêJ ‚| || tdƒ d  |tƒdttƒ ttƒ t|t ƒ tt|ƒttƒ ƒ dt|ƒ  ddfksTJ ‚| |tdƒd  || d  |tƒt t|ƒ t|t ƒ tt|ƒttƒ ƒ ttƒtdt ƒ  ddfksÆJ ‚| d|t|ƒ |d   |tƒtt|ƒ t|t ƒ tt|ƒttƒ ƒ ddfksJ ‚| || | |tƒt| d  t| t ƒ t|d |t ƒ ddfkslJ ‚| |d ||  |tƒd|d  td|t ƒ ddfksªJ ‚| t|ƒ|tƒdtd  ddfksÒJ ‚| td| ƒ|tƒdtd  ddfksþJ ‚| t|| ƒ|tƒdt|  |dfks*J ‚| t|||  ƒ|tƒt| | ƒ| t  |dfksfJ ‚| |t| | ƒ |tƒ|t d | dfksšJ ‚| |t| ||  ƒ |tƒt|| ƒ|t d  | dfksÞJ ‚| || t| | ƒ |tƒ||t d  | dfksJ ‚| t| t| ƒ ƒ|tƒtt|ƒ|t  ddfksTJ ‚| t| t|ƒ ƒ|tƒ|t tt |ƒ ddfksŽJ ‚| |tdƒd  td| ƒ ttdƒd ƒ |tƒtd tdƒ d  ddfksäJ ‚| |tdƒd  td| ƒ |tƒdttƒ dtd tdƒd    ddfks8J ‚| || t| | ƒ |tƒ|t | d  t|d ƒ | dfks‚J ‚| |t| |d  ƒ |tƒttƒ| ttd d|  ƒ ttdt|ƒ  ƒ dt|ƒ  ddfksìJ ‚| td|d  ƒ|tƒtdƒttƒ ttd d ƒ ttdƒt d ƒ d ddfksHJ ‚| |td|d  ƒ |tƒ|ttd|d  ƒ|tƒ t dfksŽJ ‚| |t| |d  ƒ |tƒdd|  tttdt|ƒ  ƒ dttƒ |tdƒd     ddfksøJ ‚| t| | ƒ|tƒdt|ƒ tdt ƒ tddt|ƒ ttƒ ƒ ddfk	sLJ ‚| t|ƒt| | ƒ |tddttƒt|ƒttƒ dtdƒ   ttd ƒ tdt|ƒ ttƒ ƒ ddfk	sÂJ ‚| t| | ƒt|ƒ |tƒttƒtdt ƒ tdt|ƒ ttƒ ƒ ddfk
sJ ‚| t| | ƒ|t|ƒ  |tƒttƒtd| ƒ tdt|ƒ ttƒ ƒ ddfk
srJ ‚| |d td| ƒ |tƒdtdƒ dt tdƒd   tddtdƒ ttƒ ƒ ddfk
sÔJ ‚| t|| ƒ|tƒ||d  td   |dfks
J ‚| |t|| ƒd  |tƒd|d  | d |d  t td   d| dfks^J ‚| |t|| ƒd  |tddd|d  | td |d  td    d| dfks¶J ‚| td| ƒ|tƒtd td   ddfksæJ ‚| t|| ƒ|tƒt|d  td   |dfksJ ‚| t|| ƒd |tddd|d  td  td|d  td    d| dfkstJ ‚| t|d ƒ|tddttdƒ tdƒ td d  ddfks¸J ‚| d!t|| | ƒd  |tƒ\}}}|d"td tt td| ƒ d| td| ƒ    d|d    ttd d|d       ¡ dksJJ ‚| t|| ƒt|| ƒ |tdd|d |d  td   d| dfks–J ‚| t|| ƒ| |tƒt|t | t  ƒd |dfksÔJ ‚| |tdƒ d  t|| ƒ |tƒttƒ t| t ƒt|t ƒ  |dfks*J ‚| |tdƒd  t|| ƒ |tƒ|t td#ƒ d  | t td#ƒ d   ttd#ƒd ƒ d |dfksœJ ‚| t|ƒ|tddttƒ t t ddfksÎJ ‚| t|| ƒ |tddt|ƒttƒ t t ddfksJ ‚| td||  ƒ|tƒtt| ƒ tt | ƒ t ddfksPJ ‚| t|| ƒ|tddtt|ƒ tt| ƒtt | ƒ  td  ddfks¢J ‚| t|ƒt|ƒ |tddttƒttƒ tdƒ t  ttƒ ddfksðJ ‚| |tdƒd  t|ƒ |tddttƒd$ttƒ td%ƒ dt  d&  dttdƒd    ddfks\J ‚| |d t|ƒ |tddd'd(ttƒ tj tdƒd(   td    ¡ tjks²J ‚| t|ƒd |tddttƒt d td d(  t ddfksöJ ‚| t| | ƒt|ƒ |tddt|t ƒ t |t  | dfks@J ‚| t|| ƒ|tƒ||d td   ddfkstJ ‚| tt|| ƒƒ|tƒ|t tt d|  ƒ |d td   ddfksÀJ ‚| t|| ƒ| |tƒt!|t ƒddfksðJ ‚| t|| ƒd | |tƒtd|d  td  d ƒd ddfks8J ‚| t|| ƒd |d  |tƒ|t!d| t ƒ ttd|d  td  d ƒ d  ddfksœJ ‚| t|| ƒ|tƒt|d td   ddfksÐJ ‚| t|| ƒd |tƒd|d  td  td|d  td    ddfks J ‚| t|| ƒt|| ƒ |tƒd| | t td || d  td || d    ddfks„J ‚| t|| ƒt|| ƒ |tƒ||d  |d  td   td || d  td || d    ddfksöJ ‚| t|| ƒt|| ƒ |tƒt|d |d  td   td || d  td || d    ddfksfJ ‚| | | t|| ƒ t|| ƒ |tddd|d  |d d|d  td   td   ddfksÐJ ‚| |t| | ƒ t|| ƒ |tƒ|| |d |t d   | dfks J ‚| |t| | ƒ t|| ƒ |tƒ||t  |d |t d   | dfkstJ ‚| t|d ƒ|tdd\}}}|dks J ‚|ttdƒ tdƒ td d   dksÎJ ‚| t"|| ƒ|tƒttd d|d   ƒttd|  ƒ t ddfksJ ‚| td|| ƒ|tƒdt|d td  ƒ ddfksTJ ‚| td|| ƒ|tdd||d td  tt|d td  ƒ   ddfks¦J ‚| td|| ƒ|tdd|d t|d td  ƒtt|d td  ƒ d   ddfksJ ‚| |td|| ƒ |tƒt|d td  tdƒd   ddfksJJ ‚| |td|| ƒ |tƒ||d td  tdƒd   ddfksJ ‚| |d td|| ƒ |tƒd|d  |d td  tdƒd   ddfksâJ ‚| td|t|d ||  ƒ ƒ|tddt|tt|d td  ƒ  ƒt|d td  ƒ ddfksLJ ‚| t#d|| ƒ|tƒdt|d  td  ƒ |dfksˆJ ‚| t#d|| ƒ|tdd||d  td  tt|d  td  ƒ   |dfksÞJ ‚| t#d|| ƒ|tdd|d t|d  td  ƒtt|d  td  ƒ d   |dfks@J ‚| |t#d|| ƒ |tƒt|d  td  tdƒd   |dfksˆJ ‚| |t#d|| ƒ |tƒ||d  td  tdƒd   |dfksÐJ ‚| |d t#d|| ƒ |tƒd|d  |d  td  tdƒd   |dfks$J ‚| t$d|| ƒ|tƒdt%t| ƒ tt|d td  ƒ  ddfksnJ ‚| td|| ƒ|tƒttt|d  td  ƒ | ƒt|d  td  ƒ | dfksÊJ ‚| t|| ƒd |tddd)|d  td*|d  d+|d  td   td    ddfks*J ‚| t&|d ƒt|d ƒ |tƒttt ƒ td d  ddfkspJ ‚| ||	|ƒ ||ƒ|t|	|ƒ||ƒ t dfks¦J ‚| |t&|d ƒ |	|d ƒ |tƒ|t|	|d ƒ|tƒ t dfksðJ ‚| |t&|d ƒ |	|d ƒ |tƒ|t|	|ƒ|tƒ tt ƒ t dfks@J ‚| ||	|| ƒ |tƒ|| t|	|ƒ||t ƒ t dfks‚J ‚| t|	|ƒ | ƒ|tƒdt|	|ƒ  t'|	|ƒƒ dfksÂJ ‚| t| | ƒ|	|ƒ |tƒt|	|ƒ||t ƒt dfksJ ‚| t|| ƒ|	|ƒ |tƒt|	|ƒ|| t ƒd t|	|ƒ||t ƒd  t dfks^J ‚| t|| ƒ| |tddd| t |d d|d  td   td   |dfks¶J ‚| t|| ƒ|	|ƒ |tƒt|	|ƒ|| t ƒd t|	|ƒ||t ƒd  t dfksJ ‚| t|| ƒ| |tdddd|t d   dd|t d    |dfksbJ ‚| t|| ƒ|	|ƒ |tddt(t|	|ƒ|t( | t ƒ t|	|ƒ|t(| t ƒ  d t dfksÌJ ‚| t|	|ƒƒ|tƒtt|	|ƒƒ|tƒt dfksJ ‚| t|| ƒ| |tddd| t |d d|d  td   td   ddfksZJ ‚| t|| ƒ|	|ƒ |tƒt|	|ƒ|t( | t ƒd t|	|ƒ|t(| t ƒd  t dfks¾J ‚| t|| ƒ| |tdd|d  td  |d d|d  td   td   ddfk sJ ‚| t|| | ƒd |	|ƒ |tƒ\}}}|t k sVJ ‚| t|	|ƒ|tƒd t|	|ƒ|dt( | t ƒtdt( | ƒ d  t|	|ƒ|dt( | t ƒtdt( | ƒ d   dk sÚJ ‚| t|| | ƒd |	|ƒ |tdd}t	||	|
iƒ|
tƒd |
dt( | t ƒtdt( | ƒ d  |
dt( | t ƒtdt( | ƒ d  k!srJ ‚| t|| ƒd t|| ƒ |tƒ\}}}||k!sªJ ‚| d| dd,|d  |d  d| t  td     d| dd,|d  |d  d| t  td     d| d|d |d  d| t  td     d| d|d |d  d| t  td      ¡ dk"s†J ‚| |d t|d  ƒ |tƒttƒtd  ttd d ƒ ttd ƒ d td  ttƒttd d ƒ ttd ƒ d  ddfk#sJ ‚| ||d  ||  | |	|ƒ |tƒ|t)t|	|ƒ|tƒtdfƒ |t)t|	|ƒ|tƒtƒ  |t|	|ƒ|tƒ  t dfk#s’J ‚| || ||ƒ |tƒd| t)t||ƒ|tƒt|fƒ t dfk#sÚJ ‚t*d|dkfd|dkfd-ƒ}| ||tƒd }|dt tt ƒt  k$s&J ‚||t|ƒ}|t&|ƒt&|d ƒ k$sPJ ‚t*d|dkf||dkfd| |dkfd-ƒ}| ||tƒd  ¡ }|tdt ƒdttƒ  d tdt ƒ td  k$sÈJ ‚||t|ƒ}| |t&|ƒ  |d t&|d ƒ  d|d  t&|d ƒ   ¡ dk%s J ‚t*t|ƒ|dkfd|dkft| ƒdfƒ}| ||tƒd }|tdƒtt ƒ td  dt  tt ƒt  k%s”J ‚||t|ƒ}|tdƒtd| ƒ t&|d ƒ t&|ƒ t&|d ƒ k%sÞJ ‚t*d|dkfd|dkfd-ƒ}| ||tƒd }|t*d|dkfd|dkfd-ƒt k&s4J ‚|t*dt+|dk|dkƒfd.ƒ |t*dt+|dk|dkƒfd.ƒ |t*dt+|dk|dk ƒfd.ƒ |t*dt+|dk|dk ƒfd.ƒ g}|D ]R}| ||tƒd  ,¡ d| t |tt ƒ t  |tdt ƒ t  k&s´J ‚&q´| t*dt-|dƒfd.ƒ|tƒd tt*dt-|dƒfd.ƒ|tƒk'sJJ ‚dt|td  | ƒ }||t|ƒ}| ||tƒd }|d |d   ¡ }|dk's J ‚| t|	|ƒ|dfƒ|tddtt|	|ƒ|tƒ |	dƒ k'sàJ ‚| t|	|ƒ|dfƒ|tddtd t|	|ƒ|tƒ td |	dƒ  tt.t)|	|ƒ|ƒ|dƒ  t.t)|	|ƒ|dfƒ|dƒ k(s`J ‚| t|	|||ƒ|dƒ|tƒtd t|	|||ƒ|tƒ t|	|d|ƒ  t.t)|	|||ƒ|ƒ|dƒ t dfk(sÐJ ‚| d#t|	|ƒ|dfƒ |tddd#t t|	|ƒ|tƒ d#|	dƒ  k)sJ ‚| ||	|| ƒ ||| ƒ |tdd|t|	|ƒ|t| ƒ | t||ƒ|t| ƒ|  k)s~J ‚t|	|ƒ||dd/t/|	|ƒ||dƒk)sªJ ‚| |	|ƒ||ƒ |tddt|	|ƒ||ƒ |tƒk)säJ ‚| ||	|| ƒ |tdd|t|	|ƒ|t| ƒ | k*s"J ‚| dt|ƒ t&|ƒ |tƒdtd  ddfk*sVJ ‚| dt|ƒ t&|ƒ |tdddtd d  ddfk*s’J ‚| d|d  tdt |d  ƒ t0|d ƒ t&|d0 ƒt&|d ƒ  |tƒd dt tt ƒ k*súJ ‚| || | t| ||  ƒ t&|| ƒ |tƒ|t | d  t| t ƒ t|d ƒ | dfk+sfJ ‚| td| ƒd t| | ƒ t&|ƒ d |tddd't| ƒtd d  k+sºJ ‚| t0|ƒ|tƒdt dfk+sÜJ ‚| t0|| ƒ|tƒd| t dfk,sJ ‚| t0|d! ƒ|tƒd!t dfk,s,J ‚| t0|d! ƒ|tƒdt dfk,sRJ ‚| t0|ƒt0|d! ƒ |tƒdtd1t ƒ t dfk,sŒJ ‚| t0|ƒ|t| | ƒ  |tddt|t  | dfk,sÌJ ‚| t| ƒt0|ƒt0|d! ƒ  |tddtd1t d! ƒd t dfk-sJ ‚| |	|ƒt0|d! ƒ |tƒ|	d!ƒtd1t ƒ t dfk-sVJ ‚| |	|ƒt0|| | ƒ |tƒ|	|| ƒt| t | ƒ | t dfk-s¦J ‚| |	|ƒt0|| | ƒ |tƒdt dfk-sØJ ‚| t1||dƒ|tƒd t| t ƒk.sJ ‚| t1||dƒ|tƒd t| t ƒtd  k.s8J ‚| t1|||ƒ|tƒd tt1|||ƒ|tƒk.sjJ ‚| t0|d ƒ|tddtt0|d ƒ|tƒk.sœJ ‚| t0|d d ƒ|tƒtt ƒd t dfk.sÐJ ‚| t0|d|  ƒ|tƒdtt ƒ t dfk/sJ ‚| t0|ƒd t0|d ƒd  |tƒtt0|ƒt0|d ƒ |tƒd tt ƒ dt  ddfk/sfJ ‚| t0d| dt|ƒ  ƒ|tƒtt t|ƒ ƒd t dfk/sªJ ‚| t0d| dt|ƒ  ƒ|tƒtt t|ƒ ƒd t dfk/sîJ ‚| t&|ƒ|tƒdt ddfk0sJ ‚| t&|| ƒ|tƒt| t ƒt ddfk0sDJ ‚| t&|d ƒ|tƒtt ƒt ddfk0srJ ‚| t&d| d ƒ|tƒtdt ƒt ddfk0s¦J ‚| t&d| d ƒ|tƒdt ddfk0sÒJ ‚| t&d| d ƒ|tdddt tdt ƒt  ddfk1sJ ‚| ||ƒt&|| ƒ |tƒt||ƒt&|| ƒ |tƒt dfk1sXJ ‚| t&|| ƒ||ƒ |tƒt||| ƒ|tƒt| t ƒ t dfk1s¤J ‚| t&|| ƒ||ƒ |tƒt||ƒ|tƒt dfk1sÞJ ‚| t&| | ƒ||ƒ |tƒt||ƒ|tƒt||| ƒ|tƒt| t ƒ  t dfk2s<J ‚| t&| | ƒ||ƒ |tƒd2k2sdJ ‚t t2|ƒ|tddttd dt  ƒ t2tt ƒ tdƒttd dt  td  ƒ d  ttd dt  ƒt3tt ƒ  t ddfk2sòJ ‚t t3|ƒ|tddttd dt  ƒt3tt ƒ ttd dt  ƒt2tt ƒ  tdƒttd dt  td  ƒ d  t ddfk3s~J ‚t4t|ƒ|t| ƒ g|t| ƒ t|ƒggƒ}t4dtd  td d gtd d dtd  ggƒ}t5ƒ \ t4dtd  ddftd d ddfgtd d ddfdtd  ddfggƒ}W d   ƒ n1 4sL0    Y  t5ƒ & | ||tƒ|k4stJ ‚W d   ƒ n1 4sŠ0    Y  | ||td3d4|ddfk4s´J ‚| ||tdd|k4sÎJ ‚| ||td3dd5|k4sêJ ‚d S )6Nza, b, cT©ÚpositiveÚnp)ÚintegerrA   zt, w, xÚfÚFÚgÚyÚYé   é   ©Únoconds)rI   é   é   é   é    rM   rN   é   éÿÿÿÿé   é   é   r   éþÿÿÿg      ø?éx   éûÿÿÿiøÿÿÿé   r   éýÿÿÿéüÿÿÿé*   é   é
   iñÿÿÿi   @é.   )rL   r   é   é   é@   é   é	   )r   T)rI   T)Úplaneg      Ð?iÖÿÿÿ)r   r   TF)Úlegacy_matrix)rf   rL   )6r   r   r   r	   r   r   r8   r   r?   r   r   r   Zdoitr   r"   r   r3   r2   r*   r(   r)   r   r9   r   r   r&   r%   r   r   r   ÚZeror   r    r$   r1   r7   r:   r!   r,   r   r   r   r#   r'   Úexpandr   r   r   r+   r-   r0   r/   r=   r<   )ÚLTÚILTÚaÚbÚcrB   ÚtÚwÚxrD   rE   rF   rG   rH   ÚLÚ_re   ZcondÚx1ZX1Úy1Zx2ZY1ZZ1ZMtZMsZMs_conds© ru   úP/var/www/auris/lib/python3.9/site-packages/sympy/integrals/tests/test_laplace.pyÚtest_laplace_transform"   s~    ÿÿ<þÿ22&<ÿ
B8ÿ
@ÿÿ
>ÿÿ
"@ÿÿ
0ÿ
2ÿ
>(,,<4D<::2ÿ
"(ÿ
J>ÿÿ
<ÿ
"ÿ
>ÿÿ
6ÿ
 .ÿþ0ÿ
 0ÿ

8ÿÿ
60ÿ
0ÿ
064ÿ
$ÿ
&2
ÿÿÿÿý
8ÿ>$(ÿ
"@ÿÿ
ÿÿ
D2ÿ
*ÿ
"&ÿÿÿ
&ÿÿ
$ÿ
  ÿ
4,ÿ
0H<ÿ
40ÿ
<ÿ
Dÿÿ
Bÿÿ
,4ÿ
&ÿÿÿBÿ.L:0ÿ
<ÿ
$ÿ
$ÿ
,ÿ
(8ÿ
<4ÿ
@ÿ
HH.ÿ
,ÿ
>ÿ
<ÿ
F6"ÿ
"$ÿ
"ÿ
@ÿ
ÿÿÿ
4ÿ
ÿÿÿ
,ÿ
ÿÿÿÿ
ÿ
4ÿ
ÿÿÿ
:ÿ
*.ÿ.þÿü
(2&ÿÿÿ*82ÿ.þ.ýÿü
6&ÿÿÿ
&ÿþþÿ
 ÿÿ
 *:$ÿÿþ
*:2ÿ
*üPÿÿÿ
&ÿþÿ
*ÿÿÿ
 "ÿ
&ÿÿ
ÿÿ
ÿ
ÿ
4ÿ
,ÿÿÿÿÿ02ÿ
0ÿþÿ"*&&ÿ
$ÿ
(ÿÿ>>ÿ2,4ÿ
ÿ
44"ÿÿÿÿÿ
DD$2.4,ÿ
"ÿ
(ÿÿÿÿ ÿÿÿÿ ÿ &ÿþþþÿ
ÿ&þþþÿ
0ÿ""ÿÿ$6 rw   c            
   	      s¼  t dƒ‰t ddd\} }‰t ddd\}}}}tdƒ‰tdƒ‰ ‡‡fd	d
„}‡‡fdd„}‡ ‡‡‡fdd„}|t| ˆ ƒˆ ˆƒ ƒˆ| ˆ ƒks˜J ‚|| ˆ ˆ| ƒ ƒ| ˆˆƒ t| ˆ ƒ ksÊJ ‚|tˆ ˆƒˆdƒƒˆd  ˆˆƒ ksôJ ‚|tˆ ˆƒˆdƒƒˆd ˆˆƒ ksJ ‚|dƒdtˆƒ ks8J ‚|dˆ ƒtˆƒksRJ ‚|dˆd  ƒˆtˆƒ kstJ ‚|dˆd  ƒˆd tˆƒ d ksžJ ‚|dˆ|  ƒˆ|d  tˆƒ t|ƒ ksÐJ ‚|||ˆ  ƒ|t| ˆ ƒ tˆƒ ks J ‚|ˆ|ˆ  ƒ| t| ˆ ƒ tˆƒ tˆƒ ks:J ‚||ˆ ˆ| d  ƒ|| ˆ t| ˆ ƒ tˆƒ t| ˆ ƒtˆƒ   ks’J ‚||ˆ| ˆ|   ƒ|t|ˆ ƒt|ˆ ƒ  tˆ ||  ƒ tˆƒ ||  ksîJ ‚||ˆ ˆ| ˆ|   ƒ||t|ˆ ƒ |t|ˆ ƒ   tˆ ||  ƒ tˆƒ ||  ksVJ ‚|ˆ|ˆ d  ƒˆ| ˆ d  t| ˆ ƒ tˆƒ d ksœJ ‚|dˆ|ˆ d   ƒ|d  ˆd  d| ˆ  dt|ˆ ƒ  d t| ˆ ƒ tˆƒ d|d   ksJ ‚|dˆ|ˆ |   ƒtˆƒt||ˆ ƒ || t|ƒ  ksRJ ‚|ˆ| |  ƒˆ|d  t|ˆ ƒ tˆƒ t|ƒ ks’J ‚||ˆ d ƒˆt| ˆ ƒ tˆƒ ksÂJ ‚||ˆ d ƒˆd t| ˆ ƒ tˆƒ d ksúJ ‚|ˆd ˆd d  ƒtˆƒ tˆƒ tˆƒ ks2J ‚|ddˆd d   ƒtˆƒ tˆƒ tˆƒ ksjJ ‚|||d ˆd   ƒt|ˆ ƒtˆƒ ksœJ ‚|ˆˆd |d   ƒt	|ˆ ƒtˆƒ ksÎJ ‚|||d |ˆ d   ƒt| ˆ ƒt|ˆ ƒ tˆƒ ksJ ‚||ˆ |d |ˆ d   ƒ|| t| ˆ ƒ t|ˆ ƒ | t| ˆ ƒt	|ˆ ƒ   tˆƒ ks‚J ‚|dˆd ˆd d   ƒˆtˆƒ tˆƒtˆƒ  ksÀJ ‚||ˆ |d ˆ| d  |d   ƒ|| | t|ˆ | ƒ |t	|ˆ | ƒ   t| ˆ ƒ tˆƒ ||d   ksBJ ‚||ˆd  | |d ˆd  d  ƒd|d  | t|ˆ ƒ |d | | |ˆ t	|ˆ ƒ t|ˆ ƒ   tˆƒ d|d   ksÐJ ‚||ˆd |d   ƒ|t
|ˆ ƒ tˆƒ | ks
J ‚||ˆd |d   ƒ|t|ˆ ƒtˆƒ d|  t| ˆ ƒtˆƒ d|    ksfJ ‚||ˆ ˆd |d   ƒ|t|ˆ ƒ tˆƒ ks J ‚||ˆˆ|   ƒ|tˆƒ| t| ˆ ƒtˆƒ |   ksäJ ‚|ˆd ˆd ˆd dˆ  d   ƒdtdˆ ƒ dt	dˆ ƒ  tˆƒ d tdˆ ƒ tˆƒ d k	sVJ ‚|d|ˆd  |ˆ  |  ƒdt| ˆ d|  ƒ tˆtd| | |d  ƒ d|  ƒ tˆƒ td| | |d  ƒ k	sÞJ ‚|dtˆƒ |ˆ |  ƒt|ˆ | ƒtˆƒ tt|ƒtˆƒ t|ƒ ƒ t|ƒt|ƒ  k
sDJ ‚|dt| ˆ ƒ |ˆ |  ˆ ƒd| tˆƒ ttƒ|  |tdƒd  t|ˆ | ƒ tt|ƒtˆƒ t|ƒ ƒ |tdƒd    tˆƒ |t| ƒ  k
sêJ ‚|dt|ˆ ƒ|  ƒ| t|d ˆ | ƒ t|tˆƒ t|ƒ ƒ | dttƒt|ƒ tˆƒ   tˆƒ ksbJ ‚||ˆ |tˆƒ |  ƒt|d ˆ |d  ƒ t|tˆƒ | ƒ d tˆƒ ks¾J ‚||| tˆƒ tˆƒt|ƒ  ˆ|  ƒt|ƒt|ƒ t|ˆ ƒ tt|ƒtˆƒ ƒ |t|ˆ ƒ tt|ƒtˆƒ ƒ  |t|ˆ ƒ  tˆƒ ks\J ‚|dtˆƒ t|ˆ ƒ|  ƒt|d ˆ | ƒtˆƒ t|tˆƒ t|ƒ ƒ t|ƒ ks¾J ‚||d tˆƒ|  ˆtdƒd   ƒd| tˆƒ ttƒ t|d ˆ ƒt|tˆƒ ƒ  d tˆƒ ks0J ‚||| t|ƒ ˆ|  tˆƒ tˆƒt|ƒ  ƒt|ƒt|ˆ ƒ tt|ƒtˆƒ ƒ t|ƒt|ˆ ƒ tt|ƒtˆƒ ƒ  t|ƒt|ˆ ƒ  tˆƒ ksÖJ ‚|dtˆƒt|ƒ d  ƒdt|ƒ tˆƒ ttƒ d| ˆ d tt|ƒtˆƒ ƒd  t|ˆ ƒ  d tˆƒ ksTJ ‚|dtˆƒt|ƒ d  ˆ ƒdˆ d|  t|ˆ ƒ tt|ƒtˆƒ ƒ d|  dtˆƒ ttƒt|ƒ   tˆƒ ksÖJ ‚|dtˆƒ| d  tˆƒ ƒd| ˆ t|d ˆ ƒ t|tˆƒ ƒ dtˆƒ ttƒ  tˆƒ ksDJ ‚|dtˆƒ| d  ƒ| ˆ d|d  ˆ d  t|d ˆ ƒ t|tˆƒ ƒ dtˆƒ |d ˆ d  ttƒ  tˆƒ ksÌJ ‚|tˆƒtˆƒ| d  ƒ dttƒ|d  ˆ dttƒ t|tˆƒ ƒ dt|d  ˆ ƒ |tˆƒ    |d  ˆd  d|d  ˆ d  tj  t|d ˆ ƒ d ttƒ| tˆƒ |d ˆ d  d   tˆƒ t|d  ˆ    ¡ }	|	dksÒJ ‚|dtˆƒ tˆƒ| d  ƒ tˆƒtˆƒd|d  ˆ d ttƒ| tˆƒ t|d ˆ ƒ t|tˆƒ ƒ d  d  ttƒ|     ¡ }	|	dksrJ ‚t|dtˆƒ| d  ƒƒdd|d  ˆtdƒd   d|d  ˆ d  dttƒ  ˆd|d  ˆd  d|d  ˆ  d  t|d ˆ ƒ t|tˆƒ ƒ d   tˆƒ ks,J ‚|tˆƒ| ˆtˆƒ|   ƒdt|d ˆ ƒ t|tˆƒ ƒ d tˆƒ ks‚J ‚|tˆƒ| d ˆtˆƒ| d   ƒdd|d  ˆ t|d ˆ ƒ t|tˆƒ ƒ  dttƒ | tˆƒ  tˆƒ ksJ ‚|tˆƒ| d ˆtˆƒ| d   ƒtˆƒdd|d  ˆd  d|d  ˆ  d  t|d ˆ ƒ t|tˆƒ ƒ dttƒ | tˆƒ d|d  ˆ d   d  ks¶J ‚|tˆ| ƒˆ|  ƒtˆƒt| ˆ ƒtˆƒ ttƒ t|| ƒt| ˆ ƒ tt|| ƒtˆƒ ƒ   ks,J ‚|dtˆ| ƒ ˆ|  ƒtˆƒdt|| ƒ t| ˆ ƒ tt|| ƒtˆƒ ƒ  ksŒJ ‚|dtˆd |d  ƒ ƒtˆƒtd|ˆ ƒ ksÄJ ‚|dtˆd |d  ƒ ƒtˆƒtd|ˆ ƒ ksüJ ‚|ˆd tdˆ ƒ dtˆƒ  d tdˆ ƒ ˆˆd d   ƒ t	ˆƒtˆƒ  dt	ˆd ƒ tˆd ƒ  dt	ˆd ƒ tˆd ƒ  dtˆd ƒ  dtˆd ƒ   ¡ dks¸J ‚d S )Nr?   zk, n, tT©Úrealz
a, b, c, dr@   rD   rE   c                    s   t | ˆ ˆƒS ©N©r   ©rF   ©r?   rn   ru   rv   rj   ó  s    z+test_inverse_laplace_transform.<locals>.ILTc                    s   t | ˆ ˆddS )NTr   r{   r|   r}   ru   rv   ÚILTSö  s    z,test_inverse_laplace_transform.<locals>.ILTSc                    s   t t| ˆˆƒˆˆ iƒS rz   )r   r   r|   ©rE   rD   r?   rn   ru   rv   ÚILTFù  s    ÿz,test_inverse_laplace_transform.<locals>.ILTFrJ   rM   rI   rQ   rS   ra   rV   rX   rN   r^   rd   rT   r   é   )r   r	   r   r   r+   r,   r)   r(   r&   r%   r   r   r"   r1   r   r   r2   ZHalfr   r   r8   r7   )
ÚkÚnrk   rl   rm   Údrj   r~   r€   rp   ru   r   rv   Útest_inverse_laplace_transformë  s"   ,2**"*20:8ÿ
<ÿ
Dÿ
F<ÿÿ
ÿ$ÿ
@088822DHÿ
>&Rÿ
&:
ÿþ
þ:<ÿ
:(ÿ
&Bÿ
@ÿÿÿ
Bÿÿ"Jÿþ
þÿÿJÿÿÿ8ÿÿ*."ÿÿÿÿÿ:ÿÿ&:ÿÿÿ2&&ÿþþÿÿ ÿ
ÿÿþþÿÿ6ÿÿÿÿ*ÿÿÿÿ<"ÿÿÿÿ. ÿ*þþþ&ýüüÿÿÿ "ÿÿÿþ
ýÿÿÿ<:ÿÿÿþÿÿ.ÿÿ&2ÿþÿ,:ÿ*ÿÿÿÿ0ÿÿÿ 4ÿÿ ÿÿ ÿÿDÿÿþþýÿür…   c                  C   s  ddl m}  t}tddd\}}}}tddd\}}tdƒ\}}	td	ƒ}
td
ƒ}dd„ }||tƒt|ƒ}t||
|iƒ|
|ƒks‚J ‚|t| t ƒt t|ƒt| | ƒks®J ‚|t| t ƒ|t  t|ƒt| | |  ƒt| | ƒ ksòJ ‚||t |d |t d   t|ƒt| | ƒt	|| ƒ t|ƒ ks>J ‚|t| t ƒt|  t|ƒ| | |d  t| | ƒ t
|ƒ ksŠJ ‚|t| t ƒttd d ƒ t|ƒt| | ƒtd|| ƒ ksÒJ ‚|dtttd ƒ  t|ƒt|ƒtt|ƒƒ ksJ ‚|t| t ƒt|  t|ƒ|| |d  t|| ƒ t
|ƒ ksPJ ‚|t| t ƒtdtd  ƒ t|ƒt|| ƒtd|| ƒ ks–J ‚t||| tttd |d  ƒ |   ttd |d  ƒ t|ƒ t¡ƒt|ƒt||| ƒ ksþJ ‚||| tttd |d  ƒ |   ttd |d  ƒ t|dd t¡t|ƒt||| ƒ ksfJ ‚|dtttd ƒ  t|ƒt|ƒtt|ƒƒ ksœJ ‚|ttdƒ tddgddggƒ  ¡ t|ƒtt|ƒt|ƒ dgdtd| ƒt|ƒ ggƒksþJ ‚|td |
tƒ t|ƒtt|
tƒt|d ƒ|dfƒks8J ‚|t	tƒt|ƒtt	tƒt|d ƒks`J ‚|ddt  dtd   dtd   t|ƒ| |ƒd| |dƒ  d| |dƒ  d| |dƒ  ksÄJ ‚|dtdt ƒ dtdt ƒ  t|ƒdttdt ƒt|d ƒ d| |d ƒ  ksJ ‚t	tdƒd ƒ}||tdt ƒ t|ƒ|| |d ƒ ksZJ ‚|tdt ƒt|ƒttdt ƒt|d ƒksŠJ ‚tddd}|t|t ƒt|ƒtt|t ƒt|d ƒksÆJ ‚|td |d td   t|ƒ| t|| ƒ t|ƒ | |ƒ ksJ ‚|td |
tƒd|d td     t|ƒ| t|| ƒ t|ƒ | |ƒ tt|
tƒt|d ƒ|dfƒ ks~J ‚|t|t ƒt|ddtt|t ƒt|d ƒdfks¶J ‚tdddtddddfD ]F}	|t|	t ƒt|dd}
t|
tƒr |
d n|
}
|
j| ksÐJ ‚qÐd S )Nr   )r+   za b c dTr@   zn, rrx   zt zrD   rE   c                 S   s   t t| ƒƒ t¡S rz   )r   r
   Úrewriter&   )Úexprru   ru   rv   Úsimp_hyp¬  s    z4test_inverse_laplace_transform_old.<locals>.simp_hyprI   rQ   r   r\   rJ   rS   iùÿÿÿrU   rZ   ÚrFrK   Úz)Zextended_real)Z	imaginaryZzero)Ú'sympy.functions.special.delta_functionsr+   r   r   r	   r?   r   r   r,   r%   r)   r"   r8   r1   r   r†   r7   r>   r=   Úinvr   r   r&   r   Ú
isinstanceÚtupleÚfunc)r+   rj   rk   rl   rm   r„   rƒ   r‰   rn   rŠ   rD   rE   rˆ   rq   ru   ru   rv   Ú"test_inverse_laplace_transform_old¢  s®    ,D" ÿ
&ÿ
$ÿ
6"ÿ
$ÿ
<ÿÿþÿ:ÿÿþÿ6*.ÿ
ÿ
(0ÿÿÿÿ&ÿÿ
.00 ÿ
* ÿÿ
ÿ

ÿr   c                  C   sÈ  t dƒ} t dƒ}t ddd}tt| ƒ| tƒtdtd  ƒ d t ddfksPJ ‚tt|| ƒ| tdd	ttttt	 ƒ d|ƒdt
|ƒtjkfks’J ‚ttd| ƒ| tdd	ttd ƒt ddfksÂJ ‚ttd| ƒ| tdd	tttd ƒ td  ddfksúJ ‚ttdtd  ƒ d t t|ƒ ¡ t|ƒt|ƒ ks6J ‚tttd ƒt t| dd	 t¡t| ƒt| ƒ ksnJ ‚ttttd ƒ td  t| dd	 t¡ ¡ td| ƒt| ƒ  t¡ t¡ ¡ ksÄJ ‚d S )
Nrp   rk   ÚuT)ZpolarrQ   rI   r   r   )r   r   r4   r?   r   r5   r.   r   r   r   r   r   rg   r   rh   r,   r†   r6   r3   )rp   rk   r‘   ru   ru   rv   Útest_expintî  sB    4&ÿÿÿ$ÿ
ÿÿþÿþþ ýÿr’   N)ZZsympy.integrals.laplacer   r   r   r   r   r   r   Zsympy.core.functionr	   r
   Z
sympy.corer   r   r   r   Zsympy.core.exprtoolsr   Zsympy.core.numbersr   r   r   Zsympy.core.relationalr   Zsympy.core.singletonr   Zsympy.core.symbolr   r   Zsympy.simplify.simplifyr   Z$sympy.functions.elementary.complexesr   r   Z&sympy.functions.elementary.exponentialr   r   r   Z%sympy.functions.elementary.hyperbolicr   r   r    r!   Z(sympy.functions.elementary.miscellaneousr"   Z$sympy.functions.elementary.piecewiser#   Z(sympy.functions.elementary.trigonometricr$   r%   r&   Zsympy.logic.boolalgr'   Z'sympy.functions.special.gamma_functionsr(   r)   r*   r‹   r+   r,   Z-sympy.functions.special.singularity_functionsr-   Z&sympy.functions.special.zeta_functionsr.   Z'sympy.functions.special.error_functionsr/   r0   r1   r2   r3   r4   r5   r6   Zsympy.functions.special.besselr7   r8   r9   r:   Zsympy.testing.pytestr;   r<   Zsympy.matricesr=   r>   Z	sympy.abcr?   rw   r…   r   r’   ru   ru   ru   rv   Ú<module>   sH   $(
   K
 7
K