I want to calculate the inverse Laplace transform of ,
where G(r)= , s is the variable.
My code can’t compute a Hypergeometric function function with symbolic variable, and it can’t get a concrete inverse Laplace transformation.
here is my code
phi=4;alpha=2.2;beta=4.7e-05;rho=25;r1=500;r2=400;lambda=4.0e-05;
syms s r t
G_r=0.5*r^2*(hypergeom([phi,-2/alpha],(alpha-2)/alpha,-r^(-alpha)*beta*rho*s )-1);
L_IBd00_s=exp(2*pi*lambda*( vpa(subs(G_r,r,r1),5) -vpa(subs(G_r,r,r2),5) ));
L_IBd00_inv=ilaplace(1/s*L_IBd00_s);
the result is
>> vpa(subs(G_r,r,radius),5)

ans =

125000.0*hypergeom([-0.90909, 4.0], 0.090909, -1.342e-9*s) – 125000.0

>> vpa(L_IBd00_inv,5)

ans =

0.000012253*heaviside(t)*ilaplace(exp(10.0*pi*hypergeom([-0.90909, 4.0], 0.090909, -1.342e-9*s) – 6.4*pi*hypergeom([-0.90909, 4.0], 0.090909, -2.1926e-9*s))/s, s, t)

So how can i get an expression of L_IBd00_inv with respect to the variable t?I want to calculate the inverse Laplace transform of ,
where G(r)= , s is the variable.
My code can’t compute a Hypergeometric function function with symbolic variable, and it can’t get a concrete inverse Laplace transformation.
here is my code
phi=4;alpha=2.2;beta=4.7e-05;rho=25;r1=500;r2=400;lambda=4.0e-05;
syms s r t
G_r=0.5*r^2*(hypergeom([phi,-2/alpha],(alpha-2)/alpha,-r^(-alpha)*beta*rho*s )-1);
L_IBd00_s=exp(2*pi*lambda*( vpa(subs(G_r,r,r1),5) -vpa(subs(G_r,r,r2),5) ));
L_IBd00_inv=ilaplace(1/s*L_IBd00_s);
the result is
>> vpa(subs(G_r,r,radius),5)

ans =

125000.0*hypergeom([-0.90909, 4.0], 0.090909, -1.342e-9*s) – 125000.0

>> vpa(L_IBd00_inv,5)

ans =

0.000012253*heaviside(t)*ilaplace(exp(10.0*pi*hypergeom([-0.90909, 4.0], 0.090909, -1.342e-9*s) – 6.4*pi*hypergeom([-0.90909, 4.0], 0.090909, -2.1926e-9*s))/s, s, t)

So how can i get an expression of L_IBd00_inv with respect to the variable t? I want to calculate the inverse Laplace transform of ,
where G(r)= , s is the variable.
My code can’t compute a Hypergeometric function function with symbolic variable, and it can’t get a concrete inverse Laplace transformation.
here is my code
phi=4;alpha=2.2;beta=4.7e-05;rho=25;r1=500;r2=400;lambda=4.0e-05;
syms s r t
G_r=0.5*r^2*(hypergeom([phi,-2/alpha],(alpha-2)/alpha,-r^(-alpha)*beta*rho*s )-1);
L_IBd00_s=exp(2*pi*lambda*( vpa(subs(G_r,r,r1),5) -vpa(subs(G_r,r,r2),5) ));
L_IBd00_inv=ilaplace(1/s*L_IBd00_s);
the result is
>> vpa(subs(G_r,r,radius),5)

ans =

125000.0*hypergeom([-0.90909, 4.0], 0.090909, -1.342e-9*s) – 125000.0

>> vpa(L_IBd00_inv,5)

ans =

0.000012253*heaviside(t)*ilaplace(exp(10.0*pi*hypergeom([-0.90909, 4.0], 0.090909, -1.342e-9*s) – 6.4*pi*hypergeom([-0.90909, 4.0], 0.090909, -2.1926e-9*s))/s, s, t)

So how can i get an expression of L_IBd00_inv with respect to the variable t? hypergeometric, matlab, laplace MATLAB Answers — New Questions

​