       Re: Using FindRoot in complex funtion

• To: mathgroup at smc.vnet.net
• Subject: [mg71945] Re: Using FindRoot in complex funtion
• From: Peter Pein <petsie at dordos.net>
• Date: Tue, 5 Dec 2006 06:04:43 -0500 (EST)
• Organization: 1&1 Internet AG
• References: <ekucvg\$830\$1@smc.vnet.net>

tatsec schrieb:
> The program is :
>
> dotnum=40;
> sigma=Sqrt;
> For[L=2,L<=dotnum,L=L+2,
> FindRoot[Integrate[2^(1/2-L/2)*sigma^(-1-L)*(Abs[t])^(1/2*(-1-L))*BesselK[(1-L)/2,Abs[t]/sigma]/(Sqrt[pi]*Gamma[L/2]),{t,0,y}],{y,1}];]
>
> when L<30, FindRoot can find the root of function.but L=>30,FindRoot
> can't calculation.So How to solve this problem?
> Thanks in advance for any help you can give me.
>
> Raymond
>
Hi Raymond,

have a closer look at your integrand:

in the case L==2 you get

FullSimplify[#1, t >= 0]&/@ Collect[ExpandAll[Normal[Series[
(2^(1/2 - L/2)*sigma^(-1 - L)*Abs[t]^((1/2)*(-1 - L))* BesselK[(1 - L)/2,
Abs[t]/sigma])/(Sqrt[Pi]*Gamma[L/2]) /. {sigma -> Sqrt, L -> 2},
{t, 0, 10}]]], t]

--> 1/(E^(t/Sqrt)*(4*2^(1/4)*t^2))

This can't give a value < infinity, when integrated from 0 to somewhere (David
Park mentioned this).

For other L you will get similar "results".

P²

• Prev by Date: Re: Making plots using transformation rules
• Next by Date: Re: Making plots using transformation rules
• Previous by thread: Re: Using FindRoot in complex funtion
• Next by thread: cyclic graphs