       Re: Assumptions question (1/m^x,x>1,m=Infinity)

• To: mathgroup at smc.vnet.net
• Subject: [mg31077] Re: Assumptions question (1/m^x,x>1,m=Infinity)
• From: hugo at doemaarwat.nl (BlackShift)
• Date: Sun, 7 Oct 2001 03:11:42 -0400 (EDT)
• Organization: Rijksuniversiteit Groningen
• References: <9pmdb3\$62v\$1@smc.vnet.net>
• Sender: owner-wri-mathgroup at wolfram.com

```On Sat, 6 Oct 2001 07:54:11 +0000 (UTC), Andrzej Kozlowski
<andrzej at bekkoame.ne.jp> wrotf:

>Still, there are several possible ways to deal with some problems of
>this kind. However, the example you give is just too trivial to serve as
>a good illustration and actually ought to be done by hand. Still, here

Can do it by hand of course, but since the formula's got quite large
at the moment I wanted to make the substitution it wasn't that easy to
do by hand...

>is one way one could try to use Mathematica to solve this sort of
>problem. Unfortunately it is not guaranteed to work in other cases,
>since Integrate with Assumptions on which it depends on is very erratic.
>
>Observe that:

>In:=
>Integrate[(-m^(-1 - x))*x, {m, 1, Infinity},
>     Assumptions -> {Re[x] > 1}] + Integrate[(-m^(-1 - x))*x,
>     m] /. m -> 1
>
>Out=
>0
>
>must be the result you wanted.

Hmm, that's actually quite true, I get that formula by integration, so
I can use the assumption there already (why didn't I think of that).
Nonetheless, isn't it possible to do it afterwards? That would be a
'nicer' thing to do, since it is part of a model, where m doesn't have
to be Infity in all cases

>There are other ways, but you would have to present your real problem
>first.

I think I can make mathematica do what I want now, but if you know a
method to make the assumption that m is infinity later on in the
calculations that would be very great, so here is the real problem

(x=x, MU=m)
Background: It is a model of starforming in galaxies,
Phi[M_]dM is the ratio of stars formed with a mass between M and M+dM
x is just a parameter for the model

In:= Phi[M_]=Cp*M^(-1-x)

To normalize Cp (which is just a normalization factor) I Integrate
over al possible M, from ML (lower mass, about .1 solar masses) to MU
(upper mass, about 32 solar masses)

In:= subC=Solve[Integrate[Phi[M],{M,ML,MU}]==1,Cp][]
1
Out= {Cp -> -------------}
1       1
----- - -----
x       x
ML  x   MU  x

with which I calculate further, with Cp in the expresions, in the
final result /.subC them. In some cases it is better to choose for ML
and MU the numbers 0.1 and 32, but sometimes it is algebraically
easyer to assume ML->0 or MU->Infinity, so sometimes I want to do
/.subC/.{ML->0.1,MU->32}, and sometimes
/.subC/.{ML->0.1,MU->Infinity}, depending on the equations, but the
latter isn't possible:

In:= subC/.{ML->0.1,MU->Infinity}

Out= {Cp -> Indeterminate}

Because it is indeterminate when x=0, which is not the case (x~1.3).

Is there anyway I can do this? without explicitly giving a value for
x, since that is the last valuable to enter (so I can test what value
for x is likely)

BTW, why doesn't mathematica not result with an If statement which
says that it is only indeterminate when x=0 and give Cp->0 otherwise.
That would be more logical behavoir I think.

groetjes,
hugo

```

• Prev by Date: Re: Bug in FourierTransform?
• Next by Date: AVI-Video with Mathematica V4.2
• Previous by thread: Re: Assumptions question (1/m^x,x>1,m=Infinity)
• Next by thread: Re: Assumptions question (1/m^x,x>1,m=Infinity)