o
    ª3Ih*•  ã                   @   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snJ ‚t t||ƒ|dƒ|tdd}t	|||iƒ}t
|||g d¢iƒ}|td |tƒ dtd   dt  d ks¥J ‚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srJ ‚| || tdƒ d  |tƒdttƒ ttƒ t|t ƒ tt|ƒttƒ ƒ dt|ƒ  ddfks§J ‚| |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s3J ‚| |d ||  |tƒd|d  td|t ƒ ddfksRJ ‚| t|ƒ|tƒdtd  ddfksfJ ‚| 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s°J ‚| |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s'J ‚| t| t|ƒ ƒ|tƒ|t tt |ƒ ddfksDJ ‚| |tdƒd  td| ƒ ttdƒd ƒ |tƒtd tdƒ d  ddfksoJ ‚| |tdƒd  td| ƒ |tƒdttƒ dtd tdƒd    ddfks™J ‚| || 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s!J ‚| |td|d  ƒ |tƒ|ttd|d  ƒ|tƒ t dfksDJ ‚| |t| |d  ƒ |tƒdd|  tttdt|ƒ  ƒ dttƒ |tdƒd     ddfksyJ ‚| t| | ƒ|tƒdt|ƒ tdt ƒ tddt|ƒ ttƒ ƒ ddfks£J ‚| 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s6J ‚| |d td| ƒ |tƒdtdƒ dt tdƒd   tddtdƒ ttƒ ƒ ddfksgJ ‚| 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s7J ‚| t|d ƒ|tddttdƒ tdƒ td d  ddfksYJ ‚| d!t|| | ƒd  |tƒ\}}}|d"td tt td| ƒ d| td| ƒ    d|d    ttd d|d       ¡ dks¢J ‚| 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sKJ ‚| t|ƒ|tddttƒ t t ddfksdJ ‚| 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s¥J ‚| 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sVJ ‚| t|ƒd |tddttƒt d td d(  t ddfksxJ ‚| t| | ƒt|ƒ |tddt|t ƒ t |t  | dfksJ ‚| t|| ƒ|tƒ||d td   ddfks·J ‚| 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	sJ ‚| t|| ƒd |d  |tƒ|t!d| t ƒ ttd|d  td  d ƒ d  ddfk	sKJ ‚| t|| ƒ|tƒt|d td   ddfk	seJ ‚| t|| ƒd |tƒd|d  td  td|d  td    ddfk	sJ ‚| 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
s0J ‚| | | t|| ƒ t|| ƒ |tddd|d  |d d|d  td   td   ddfk
seJ ‚| |t| | ƒ t|| ƒ |tƒ|| |d |t d   | dfk
sJ ‚| |t| | ƒ t|| ƒ |tƒ||t  |d |t d   | dfk
s·J ‚| 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s'J ‚| td|| ƒ|tdd||d td  tt|d td  ƒ   ddfksPJ ‚| 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s¢J ‚| |td|| ƒ |tƒ||d td  tdƒd   ddfksÅJ ‚| |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s#J ‚| t#d|| ƒ|tƒdt|d  td  ƒ |dfksAJ ‚| t#d|| ƒ|tdd||d  td  tt|d  td  ƒ   |dfkslJ ‚| t#d|| ƒ|tdd|d t|d  td  ƒtt|d  td  ƒ d   |dfksJ ‚| |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s4J ‚| td|| ƒ|tƒttt|d  td  ƒ | ƒt|d  td  ƒ | dfksbJ ‚| 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sµJ ‚| ||	|ƒ ||ƒ|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s.J ‚| t|| ƒ|	|ƒ |tddt(t|	|ƒ|t( | t ƒ t|	|ƒ|t(| t ƒ  d t dfkscJ ‚| 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sªJ ‚| 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s(J ‚| t|	|ƒ|tƒd t|	|ƒ|dt( | t ƒtdt( | ƒ d  t|	|ƒ|dt( | t ƒtdt( | ƒ d   dksjJ ‚| t|| | ƒd |	|ƒ |tdd}t	||	|
iƒ|
tƒd |
dt( | t ƒtdt( | ƒ d  |
dt( | t ƒtdt( | ƒ d  ks¶J ‚| 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s%J ‚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saJ ‚||t|ƒ}| |t&|ƒ  |d t&|d ƒ  d|d  t&|d ƒ   ¡ dksJ ‚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sJ ‚|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 ])}| ||tƒd  ,¡ d| t |tt ƒ t  |tdt ƒ t  ksJ ‚qW| t*dt-|dƒfd.ƒ|tƒd tt*dt-|dƒfd.ƒ|tƒks¢J ‚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seJ ‚| 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s(J ‚| dt|ƒ t&|ƒ |tdddtd d  ddfksFJ ‚| d|d  tdt |d  ƒ t0|d ƒ t&|d0 ƒt&|d ƒ  |tƒd dt tt ƒ kszJ ‚| || | t| ||  ƒ t&|| ƒ |tƒ|t | d  t| t ƒ t|d ƒ | dfks°J ‚| 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s&J ‚| t0|ƒt0|d! ƒ |tƒdtd1t ƒ t dfksCJ ‚| t0|ƒ|t| | ƒ  |tddt|t  | dfkscJ ‚| 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s¨J ‚| |	|ƒ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sJ ‚| t1|||ƒ|tƒd tt1|||ƒ|tƒks2J ‚| t0|d ƒ|tddtt0|d ƒ|tƒksKJ ‚| t0|d d ƒ|tƒtt ƒd t dfkseJ ‚| 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s°J ‚| 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sJ ‚| t&|| ƒ|tƒt| t ƒt ddfksJ ‚| t&|d ƒ|tƒtt ƒt ddfks6J ‚| t&d| d ƒ|tƒtdt ƒt ddfksPJ ‚| t&d| d ƒ|tƒdt ddfksfJ ‚| t&d| d ƒ|tdddt tdt ƒt  ddfks†J ‚| ||ƒt&|| ƒ |tƒt||ƒt&|| ƒ |tƒt dfks©J ‚| t&|| ƒ||ƒ |tƒt||| ƒ|tƒt| t ƒ t dfksÏJ ‚| t&|| ƒ||ƒ |tƒt||ƒ|tƒt dfksìJ ‚| t&| | ƒ||ƒ |tƒt||ƒ|tƒt||| ƒ|tƒt| t ƒ  t dfksJ ‚| t&| | ƒ||ƒ |tƒd2ks/J ‚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svJ ‚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s¼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 s#w   Y  t5ƒ  | ||tƒ|ks7J ‚W d   ƒ n	1 sBw   Y  | ||td3d4|ddfksWJ ‚| ||tdd|ksdJ ‚| ||td3dd5|ksrJ ‚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   Ú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   ÚcondÚx1ÚX1Úy1Úx2ÚY1ÚZ1ÚMtÚMsÚMs_conds© r~   úb/home/air/sanwanet/gpt-api/venv/lib/python3.10/site-packages/sympy/integrals/tests/test_laplace.pyÚtest_laplace_transform"   sh   ÿÿ<þ00&<ÿ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ÿ""ÿÿÿÿ  r€   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sLJ ‚|| ˆ ˆ| ƒ ƒ| ˆˆƒ t| ˆ ƒ kseJ ‚|tˆ ˆƒˆdƒƒˆd  ˆˆƒ kszJ ‚|tˆ ˆƒˆdƒƒˆd ˆˆƒ ksŽJ ‚|dƒdtˆƒ ksšJ ‚|dˆ ƒtˆƒks¦J ‚|dˆd  ƒˆtˆƒ ks¶J ‚|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sBJ ‚||ˆ| ˆ|   ƒ|t|ˆ ƒt|ˆ ƒ  tˆ ||  ƒ tˆƒ ||  kspJ ‚||ˆ ˆ| ˆ|   ƒ||t|ˆ ƒ |t|ˆ ƒ   tˆ ||  ƒ tˆƒ ||  ks¤J ‚|ˆ|ˆ 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s"J ‚|ˆ| |  ƒˆ|d  t|ˆ ƒ tˆƒ t|ƒ ksBJ ‚||ˆ d ƒˆt| ˆ ƒ tˆƒ ksZJ ‚||ˆ d ƒˆd t| ˆ ƒ tˆƒ d ksvJ ‚|ˆd ˆd d  ƒtˆƒ tˆƒ tˆƒ ks’J ‚|ddˆd d   ƒtˆƒ tˆƒ tˆƒ ks®J ‚|||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sYJ ‚||ˆ |d ˆ| d  |d   ƒ|| | t|ˆ | ƒ |t	|ˆ | ƒ   t| ˆ ƒ tˆƒ ||d   ksšJ ‚||ˆ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s,J ‚||ˆ ˆd |d   ƒ|t|ˆ ƒ tˆƒ ksIJ ‚||ˆˆ|   ƒ|tˆƒ| t| ˆ ƒtˆƒ |   kskJ ‚|ˆd ˆd ˆd dˆ  d   ƒdtdˆ ƒ dt	dˆ ƒ  tˆƒ d tdˆ ƒ tˆƒ d ks¤J ‚|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sJ ‚|dt| ˆ ƒ |ˆ |  ˆ ƒd| tˆƒ ttƒ|  |tdƒd  t|ˆ | ƒ tt|ƒtˆƒ t|ƒ ƒ |tdƒd    tˆƒ |t| ƒ  ksnJ ‚|dt|ˆ ƒ|  ƒ| t|d ˆ | ƒ t|tˆƒ t|ƒ ƒ | dttƒt|ƒ tˆƒ   tˆƒ ksªJ ‚||ˆ |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sXJ ‚||d tˆƒ|  ˆtdƒd   ƒd| tˆƒ ttƒ t|d ˆ ƒt|tˆƒ ƒ  d tˆƒ ks‘J ‚||| 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s#J ‚|dtˆƒt|ƒ d  ˆ ƒdˆ d|  t|ˆ ƒ tt|ƒtˆƒ ƒ d|  dtˆƒ ttƒt|ƒ   tˆƒ ksdJ ‚|dtˆƒ| d  tˆƒ ƒd| ˆ t|d ˆ ƒ t|tˆƒ ƒ dtˆƒ ttƒ  tˆƒ ks›J ‚|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sbJ ‚|dtˆƒ tˆƒ| d  ƒ tˆƒtˆƒd|d  ˆ d ttƒ| tˆƒ t|d ˆ ƒ t|tˆƒ ƒ d  d  ttƒ|     ¡ }	|	dks²J ‚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
swJ ‚|ˆ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?   ro   r~   r   rk   ó  s   z+test_inverse_laplace_transform.<locals>.ILTc                    s   t | ˆ ˆddS )NTr   r„   r…   r†   r~   r   ÚILTSö  s   z,test_inverse_laplace_transform.<locals>.ILTSc                    s   t t| ˆˆƒˆˆ iƒS rƒ   )r   r   r…   ©rE   rD   r?   ro   r~   r   Ú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   ÚHalfr   r   r8   r7   )
ÚkÚnrl   rm   rn   Údrk   r‡   r‰   rq   r~   rˆ   r   Útest_inverse_laplace_transformë  sø  ,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sAJ ‚|t| t ƒt t|ƒt| | ƒksWJ ‚|t| t ƒ|t  t|ƒt| | |  ƒt| | ƒ ksyJ ‚||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s%J ‚|t| t ƒtdtd  ƒ t|ƒt|| ƒtd|| ƒ ksHJ ‚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s°J ‚|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sJ ‚|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sªJ ‚|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sXJ ‚tdddtddddfD ]#}	|t|	t ƒt|dd}
t|
tƒr}|
d n|
}
|
j| ks‡J ‚qed S )Nr   )r+   za b c dTr@   zn, rr   zt zrD   rE   c                 S   s   t t| ƒƒ t¡S rƒ   )r   r
   Úrewriter&   )Úexprr~   r~   r   Úsimp_hyp¬  s   z4test_inverse_laplace_transform_old.<locals>.simp_hyprI   rQ   r   r\   rJ   rS   iùÿÿÿrU   rZ   ÚrFrK   Úz)Úextended_real)Ú	imaginaryÚ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+   rk   rl   rm   rn   rŽ   r   r“   ro   r”   rD   rE   r’   rr   r~   r~   r   Ú"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s(J ‚tt|| ƒ| tdd	ttttt	 ƒ d|ƒdt
|ƒtjkfksIJ ‚ttd| ƒ| tdd	ttd ƒt ddfksaJ ‚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sšJ ‚tttd ƒt t| dd	 t¡t| ƒt| ƒ ksµJ ‚ttttd ƒ td  t| dd	 t¡ ¡ td| ƒt| ƒ  t¡ t¡ ¡ ksßJ ‚d S )
Nrq   rl   ÚuT)ÚpolarrQ   rI   r   r   )r   r   r4   r?   r   r5   r.   r   r   r   r   r   rh   r   ri   r,   r   r6   r3   )rq   rl   rž   r~   r~   r   Útest_expintî  s:   4&ÿÿÿ$ÿÿþþ
 ýr    N)ZÚsympy.integrals.laplacer   r   r   r   r   r   r   Úsympy.core.functionr	   r
   Ú
sympy.corer   r   r   r   Úsympy.core.exprtoolsr   Úsympy.core.numbersr   r   r   Úsympy.core.relationalr   Úsympy.core.singletonr   Úsympy.core.symbolr   r   Úsympy.simplify.simplifyr   Ú$sympy.functions.elementary.complexesr   r   Ú&sympy.functions.elementary.exponentialr   r   r   Ú%sympy.functions.elementary.hyperbolicr   r   r    r!   Ú(sympy.functions.elementary.miscellaneousr"   Ú$sympy.functions.elementary.piecewiser#   Ú(sympy.functions.elementary.trigonometricr$   r%   r&   Úsympy.logic.boolalgr'   Ú'sympy.functions.special.gamma_functionsr(   r)   r*   r˜   r+   r,   Ú-sympy.functions.special.singularity_functionsr-   Ú&sympy.functions.special.zeta_functionsr.   Ú'sympy.functions.special.error_functionsr/   r0   r1   r2   r3   r4   r5   r6   Úsympy.functions.special.besselr7   r8   r9   r:   Úsympy.testing.pytestr;   r<   Úsympy.matricesr=   r>   Ú	sympy.abcr?   r€   r   r   r    r~   r~   r~   r   Ú<module>   sJ   $ (
   K
 7
K