o
    ÏeÝp  ã                   @   s¤  d dl 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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)m*Z*m+Z+ d dl,m-Z-m.Z. d dl/m0Z0m1Z1 d dl2m3Z3 d dl4m5Z5m6Z6m7Z7m8Z8m9Z9m:Z:m;Z;m<Z< d dl=m>Z>m?Z?m@Z@mAZA d dlBmCZCmDZD d dlEmFZFmGZG d dlHmIZI eCdd„ ƒZJeCdd„ ƒZKeCdd„ ƒZLdS )é    )Úlaplace_transformÚinverse_laplace_transformÚLaplaceTransformÚInverseLaplaceTransform)ÚFunctionÚ
expand_mul)Ú
EulerGammaÚSubsÚ
DerivativeÚdiff)Úfactor_terms)ÚIÚooÚpi)ÚS)ÚSymbolÚsymbols©Úsimplify)ÚAbsÚre)ÚexpÚlogÚ	exp_polar)ÚcoshÚsinhÚcothÚasinh)Úsqrt)ÚatanÚcosÚsin)Ú
lowergammaÚgamma)Ú
DiracDeltaÚ	Heaviside)Úlerchphi)ÚfresnelcÚfresnelsÚerfÚerfcÚEiÚCiÚexpintÚE1)ÚbesseliÚbesseljÚbesselkÚbessely)ÚslowÚwarns_deprecated_sympy)ÚMatrixÚeye)Úsc               	   C   s5  t } tddd\}}}tdƒ\}}}tdƒ}tdƒ}dtt|ƒ|tƒ d  ¡ d	dtd   ks2J ‚| || |d  |td
ƒd   |tƒ|td  dtd   dtt	ƒ dttdƒd     ddfksgJ ‚| |||  |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sJ ‚| |tdƒd  || d	  |tƒt	 t|ƒ t|t ƒ tt|ƒttƒ ƒ tt	ƒtdt ƒ  ddfksRJ ‚| 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s¥J ‚| |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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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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sYJ ‚| |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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sAJ ‚| t| | ƒ|t|ƒ  |tƒtt	ƒtd| ƒ tdt|ƒ ttƒ ƒ ddfksnJ ‚| tdt|| ƒ ƒ|tƒdt tt	ƒt|ƒ t|t ƒ tt|ƒtdt ƒ ƒ ttdƒd    ddfks¨J ‚| tdt|| ƒ ƒt|ƒ |tƒt|t ƒtt|ƒtdt ƒ ƒ tt	ƒtdt ƒ  ddfksÜJ ‚| |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sRJ ‚| |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sÝJ ‚| 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sHJ ‚| t|| ƒt|| ƒ |tdd|d|d  td   d| dfksnJ ‚| t|| ƒ| |tƒt|t | t  ƒd |dfksJ ‚| |tdƒ d  t|| ƒ |tƒtt	ƒ t| t ƒt|t ƒ  |dfks¸J ‚| tdt|| ƒ ƒ|tƒtt	ƒt|ƒ t|t ƒ ttdƒd   ddfksâJ ‚| t|ƒtdt|| ƒ ƒ |tddt|ƒ ttd
ƒd   tt	ƒtd  d| t  t|t ƒ tt|ƒtdt ƒ ƒ d  ttdƒd   ddfks3J ‚| tdt|| ƒ ƒt|ƒ |tƒtt	ƒt|t ƒ tt|ƒtdt ƒ ƒ ttƒ ddfkseJ ‚| tt|| ƒƒd t|ƒ |tƒtt	ƒt|t ƒd  dttƒ  ddfksJ ‚| |tdƒd  t|| ƒ |tƒ|t tdƒ d  | t tdƒ d   ttdƒd ƒ d |dfksÈJ ‚| tdt|| ƒ ƒ|tƒtt	ƒt|ƒ t|t ƒ tt|ƒtdt ƒ ƒ ttdƒd   dt  ddfksJ ‚| t|ƒtdt|| ƒ ƒ |tƒtt	ƒ|td   t|t ƒ ttd
ƒd   ddfks2J ‚| tdt|| ƒ ƒt|ƒ |tƒtt	ƒt|t ƒ ttƒ ddfksXJ ‚| tt|| ƒƒd t|ƒ |tƒtt	ƒt|t ƒd  dttƒ  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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	sbJ ‚| |d t|ƒ |tddd!d"ttƒ tj tdƒd"   td    ¡ tjk	sJ ‚| 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	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
sPJ ‚| t|| ƒd |d  |tƒ|td| t ƒ ttd|d  td  d ƒ d  ddfk
s‚J ‚| tdt|| ƒ ƒ|tƒtt	ƒt|ƒ t| t ƒ ttdƒd   ddfk
s­J ‚| tdt|| ƒ ƒ| |tƒt	tt|ƒtdt ƒ ƒ 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dt|| ƒ ƒ |tddtt	ƒ| td   t| t ƒ ttd
ƒd   ddfksGJ ‚| tdt|| ƒ ƒt|ƒ |tƒtt	ƒtdt ƒ t| t ƒ ddfkspJ ‚| 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sJ ‚| | | t|| ƒ t|| ƒ |tddd|d  |d d|d  td   td   ddfksHJ ‚| |t| | ƒ t|| ƒ |tƒ|| |d |t d   | dfkspJ ‚| |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t|| ƒƒ|tƒt|ƒtt|t ƒ  ddfksJ ‚| t|| ƒtt|| ƒƒ |tddt|ƒ ttƒ|t   |dfks2J ‚| tt|| ƒd ƒ|tdddt tt|ƒ ttƒ ƒt  ddfksYJ ‚| tt|| ƒƒ|tddt|ƒ tt|t ƒ  dt  | dfksJ ‚| t|| ƒtt|| ƒƒ |tƒdt|ƒttƒ t  ddfks£J ‚| tt|| ƒd ƒ|tƒtt|ƒ ttƒ ƒ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s
J ‚| td|| ƒ|tdd|d t|d td  ƒtt|d td  ƒ d   ddfks9J ‚| |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dt|| ƒ ƒ|tƒt| t ƒt ddfksÆJ ‚| |tdƒd  tddt|| ƒ ƒ |tƒ|tdƒd  t| t ƒ t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sIJ ‚| td|| ƒ|tdd||d  td  tt|d  td  ƒ   |dfkstJ ‚| 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ƒd  tddt|| ƒ ƒ |tƒ|tdƒd  t|t ƒ td  ddfksFJ ‚| td|| ƒ|tƒdtt| ƒ t	t|d td  ƒ  ddfkskJ ‚| 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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sTJ ‚| |||| ƒ |tƒ|| t||ƒ||t ƒ t dfksuJ ‚| t||ƒ | ƒ|tƒdt||ƒ  t!||ƒƒ dfks•J ‚| t| | ƒ||ƒ |tƒt||ƒ||t ƒt dfksµJ ‚| t| | ƒtt|| ƒd ƒ |tƒ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sjJ ‚| 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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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soJ ‚| 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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sXJ ‚| |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sJ ‚| ||d  ||  | ||ƒ |tƒ|t#t||ƒ|tƒtdfƒ |t#t||ƒ|tƒtƒ  |t||ƒ|tƒ  t df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svJ ‚| 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s9J ‚| dt|ƒ t |ƒ |tdddtd d  ddfksWJ ‚| || | t| ||  ƒ t || ƒ |tƒ|t | d  t| t ƒ t|d ƒ | dfksJ ‚| td| ƒd t| | ƒ t |ƒ d |tddd!t| ƒtd d  ks·J ‚| t(|ƒ|tƒdt dfksÈJ ‚| t(|| ƒ|tƒd| t dfksÝJ ‚| t(|d ƒ|tƒdt dfksðJ ‚| t(|d ƒ|tƒdt dfksJ ‚| t(|ƒt(|d ƒ |tƒdtd(t ƒ t dfks J ‚| t(|ƒ|t| | ƒ  |tddt|t  | dfks@J ‚| t| ƒt(|ƒt(|d ƒ  |tddtd(t d ƒd t dfksfJ ‚| ||ƒt(|d ƒ |tƒ|dƒtd(t ƒ t dfks…J ‚| ||ƒt(|| | ƒ |tƒ||| ƒt| t | ƒ | t dfks­J ‚| ||ƒt(|| | ƒ |tƒdt dfksÆJ ‚| t(|d ƒ|tdd&tt(|d ƒ|tƒksßJ ‚| t(|d d ƒ|tƒtt ƒd t dfksùJ ‚| t(|d|  ƒ|tƒdtt ƒ t dfksJ ‚| t(|ƒd t(|d ƒd  |tƒtt(|ƒt(|d ƒ |tƒd tt ƒ dt  ddfksDJ ‚| t(d| dt|ƒ  ƒ|tƒtt t|ƒ ƒd t dfksfJ ‚| t(d| 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sÊJ ‚| t d| d ƒ|tƒtdt ƒt ddfksäJ ‚| t d| d ƒ|tƒdt ddfksúJ ‚| 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ddttd dt	  ƒ t)tt	 ƒ tdƒttd dt	  t	d  ƒ d  ttd dt	  ƒt*tt	 ƒ  t ddfks„J ‚t t*|ƒ|tddttd dt	  ƒt*tt	 ƒ ttd dt	  ƒt)tt	 ƒ  tdƒttd dt	  t	d  ƒ d  t ddfksÊJ ‚t+t|ƒ|t| ƒ g|t| ƒ t|ƒggƒ}t+dtd  td d gtd d dtd  ggƒ}t,ƒ . t+d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 s1w   Y  t,ƒ  | ||tƒ|ksEJ ‚W d   ƒ n	1 sPw   Y  | ||td)d*|ddfkseJ ‚| ||tdd&|ksrJ ‚| ||td)dd+|ks€J ‚d S ),Nza, b, cT©Zpositivezt, w, xÚfÚgé   é   éÿÿÿÿé   é   é   é   é   r   éþÿÿÿg      ø?éx   éûÿÿÿé   iøÿÿÿé   r   éýÿÿÿéüÿÿÿé*   é   é	   é
   iñÿÿÿi   @é.   )Únocondsr   é   é   é@   é   ©rO   )ÚplaneiÖÿÿÿF)Úlegacy_matrix)rV   rO   )-r   r   r   r   r   r7   Zdoitr   r   r   r+   r*   r"   r#   r   r1   r   r   r!   r    r   r   r)   r   ÚZeror   r   r   r0   r/   r2   r   r%   r   r   r
   r   r	   r   r   r$   r(   r'   r5   r4   )ÚLTÚaÚbÚcÚtÚwÚxr9   r:   ÚLÚ_rU   ZcondZMtZMsZMs_conds© ra   úRD:\Projects\ConvertPro\env\Lib\site-packages\sympy/integrals/tests/test_laplace.pyÚtest_laplace_transform   sH  0&<ÿ@8ÿ@ÿÿ>ÿ
ÿ"@ÿ
ÿ0
ÿ2
ÿ>(,,<4D<2
ÿ"(
ÿJ>ÿ
ÿ<
ÿ"
ÿ>ÿ
ÿ6
ÿ .ÿþ0
ÿ 0
ÿ8ÿÿ
ÿ":ÿ8ÿ
ÿ60
ÿ0
ÿ064
ÿ$
ÿ&2
ÿÿÿ
ý8ÿ>$(
ÿ0
ÿ&>ÿÿÿÿ
ÿ"8
ÿ"(
ÿ"@ÿ
ÿBÿÿ
ÿ"4
ÿ" 
ÿ"(
ÿÿ
ÿD2
ÿ*
ÿ"&ÿÿ
ÿ&ÿ
ÿ$
ÿ  
ÿ4,
ÿ0H<
ÿ2
ÿH40
ÿ&8
ÿ"&
ÿ<
ÿDÿ
ÿBÿ
ÿ,4
ÿ&ÿÿÿBÿ.L<&
ÿ&
ÿ(
ÿ"
ÿB:0
ÿ<
ÿ$
ÿ$
ÿ,
ÿ<,*
ÿ(8
ÿ<4
ÿ@
ÿHH.
ÿ,(
ÿ,
ÿ>
ÿ<
ÿF6"
ÿ"$
ÿ"
ÿ@
ÿ((
ÿÿÿ
ÿ4
ÿÿÿ
ÿ,
ÿÿÿÿ
ÿ
ÿ4
ÿÿÿ
ÿ:
ÿ*.ÿ.þÿ
ü*82ÿ.þ.ý
ü6&ÿÿ
ÿ&ÿþþ
ÿ
ÿ&ÿþ
ÿ*ÿÿ
ÿ "
ÿ&ÿ
ÿÿ
ÿ
ÿ
ÿ4
ÿ02
ÿ0ÿ
þ"*&&
ÿ$
ÿ(
ÿ>>ÿ2
ÿ44"ÿÿÿÿ
ÿDD$2.4,
ÿ"
ÿ &ÿþþþ
ÿÿ&þþþ
ÿ0ÿ""ÿÿÿÿ  rc   c                  C   s¸  ddl m}  t}tddd\}}}}tddd\}}tdƒ\}}	td	ƒ}
d
d„ }|dt|ƒ| |ƒks4J ‚|dt t|ƒt|ƒksBJ ‚|||t  t|ƒ|t| | ƒ t|ƒ ks[J ‚|t|t  t|ƒ| t| | ƒ t|ƒ | |ƒ ksyJ ‚|t|t d  t|dd|| | d  t| | ƒ t|ƒ d ksŸJ ‚|dt|t d   t|dd|d  |d  d| |  dt|| ƒ  d t| | ƒ t|ƒ d|d   ksÛJ ‚|dt|t |   t|ƒt|ƒt||| ƒ || t	|ƒ  ksþJ ‚|t| |  t|ƒ||d  t|| ƒ t|ƒ t	|ƒ ks J ‚||t d t|ƒ|t| | ƒ t|ƒ ks:J ‚||t d t|ƒ|d t| | ƒ t|ƒ d ksXJ ‚|||d td   t|ƒt
|| ƒt|ƒ kssJ ‚|ttd |d   t|ƒt|| ƒt|ƒ ksŽJ ‚|||d |t d   t|ƒt| | ƒt
|| ƒ t|ƒ ks²J ‚||t |d |t d   t|ƒ|| t| | ƒ t
|| ƒ | t| | ƒt|| ƒ   t|ƒ 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sJJ ‚|t| t ƒt|  t|ƒ| | |d  t| | ƒ t	|ƒ kspJ ‚|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 ‚|dtd td d   t|ƒ|t|ƒ t
|ƒt|ƒ  ksÐJ ‚|td td d  t|ƒt
|ƒ t|ƒ | |ƒ ksîJ ‚|ddtd d   t|ƒt
|ƒ t|ƒ | |ƒ ksJ ‚|dtd  t|ƒ|t|ƒ ksJ ‚|dtd  t|ƒ|d t|ƒ d ks6J ‚|dt|  t|ƒ||d  t|ƒ t	|ƒ ksQJ ‚|td td td dt  d   t|dddt
d| ƒ dtd| ƒ  t|ƒ d td| ƒ t|ƒ d ksŽJ ‚|||td |d   t|ƒƒt|| ƒt|ƒ ks«J ‚||ttd |d   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sJ ‚t||| tttd |d  ƒ |   ttd |d  ƒ t|ƒ t¡ƒt|ƒt||| ƒ ksCJ ‚||| tttd |d  ƒ |   ttd |d  ƒ t|dd t¡t|ƒt||| ƒ kswJ ‚|dtttd ƒ  t|ƒt|ƒtt|ƒƒ ks’J ‚|dtd td d   t|dd|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 |d   t|dd|t|| ƒ t|ƒ | ksJ ‚||td |d   t|ƒ|t|| ƒt|ƒ d|  t| | ƒt|ƒ d|    ks3J ‚||t td |d   t|dd|t|| ƒ t|ƒ ksTJ ‚||tt|   t|ƒ|t|ƒ| t| | ƒt|ƒ |   ksxJ ‚||t t| d  t|ƒ|| | t| | ƒ t|ƒ t| | ƒt|ƒ   ks¦J ‚||t| t|   t|dd|t|| ƒt|| ƒ  t| ||  ƒ t|ƒ ||  ksØJ ‚||t t| t|   t|dd||t|| ƒ |t|| ƒ   t| ||  ƒ t|ƒ ||  ksJ ‚||t |d t| d  |d   t|dd|| | t
|| | ƒ |t|| | ƒ   t| | ƒ t|ƒ ||d   ksUJ ‚|td |
tƒ t|ƒtt|
tƒt|d ƒ|dfƒksrJ ‚||td  | |d td  d  t|dd|d | | t|| ƒ d|d  | t
|| ƒ  || | t|| ƒ  |t
|| ƒ  t|ƒ d|d   ksÇJ ‚|ttƒt|ƒtttƒt|d ƒksÛJ ‚|dt|ƒ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s5J ‚|tdt ƒt|ƒttdt ƒt|d ƒksMJ ‚tddd}|t|t ƒt|ƒtt|t ƒt|d ƒkskJ ‚|td |d td   t|ƒ| t
|| ƒ t|ƒ | |ƒ ksJ ‚|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 ]#}	|t|	t ƒt|dd}
t|
tƒr|
d n|
}
|
j| ksJ ‚qð|d|td  |t  |  t|ƒdt| | d|  ƒ t
|td| | |d  ƒ d|  ƒ t|ƒ td| | |d  ƒ ksZJ ‚d S )#Nr   )r$   za b c dTr8   zn, r)Úrealzt zr9   c                 S   s   t t| ƒƒ t¡S )N)r   r   Úrewriter!   )Úexprra   ra   rb   Úsimp_hyp–  s   z0test_inverse_laplace_transform.<locals>.simp_hypr<   r?   r   rF   r;   rC   rE   rQ   r>   rA   rM   rL   r@   rJ   iùÿÿÿrB   rH   ÚrFrT   Úz)Zextended_real)Z	imaginaryÚzero)Ú'sympy.functions.special.delta_functionsr$   r   r   r   r7   r%   r   r"   r#   r!   r    r   r0   r)   r   r   r   re   r/   r6   r5   Úinvr
   r   r   Ú
isinstanceÚtupleÚfunc)r$   ZILTrY   rZ   r[   ÚdÚnrh   r\   ri   r9   rg   ra   ra   rb   Útest_inverse_laplace_transform  s<  2<*ÿ0ÿÿ
ÿÿ$ÿD4<66H"H
ÿ.B" 
ÿ&
ÿ$
ÿ6
ÿ<<&.6.B
ÿ::"
ÿ$
ÿ<ÿÿ
þ:ÿ
þ6"
ÿ*.
ÿ
ÿ<
ÿ"
ÿ(
ÿ8
ÿ<
ÿ"D
ÿ.R
ÿ
ÿ.Nÿÿ
ÿ
ÿ(&ÿ
ÿ.00 
ÿ* ÿ
ÿ
ÿ
ÿ"@ÿÿÿrr   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 )
Nr^   rY   ÚuT)Zpolarr<   r;   r   r   )r   r   r,   r7   r   r-   r&   r   r   r   r   r   rW   r   Úexpandr%   re   r.   r+   )r^   rY   rs   ra   ra   rb   Útest_expint  s:   4&ÿÿÿ$ÿÿþþ
 ýru   N)MZsympy.integrals.laplacer   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.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.trigonometricr   r    r!   Z'sympy.functions.special.gamma_functionsr"   r#   rk   r$   r%   Z&sympy.functions.special.zeta_functionsr&   Z'sympy.functions.special.error_functionsr'   r(   r)   r*   r+   r,   r-   r.   Zsympy.functions.special.besselr/   r0   r1   r2   Zsympy.testing.pytestr3   r4   Zsympy.matricesr5   r6   Z	sympy.abcr7   rc   rr   ru   ra   ra   ra   rb   Ú<module>   s<    (
  s
 