       Re: real world solutions for a fractional permutation

• To: mathgroup at smc.vnet.net
• Subject: [mg77033] Re: real world solutions for a fractional permutation
• From: dimitris <dimmechan at yahoo.com>
• Date: Fri, 1 Jun 2007 02:41:44 -0400 (EDT)
• References: <f3lud3\$6nq\$1@smc.vnet.net>

```Solve is not intended for this kind of solution.

For equations like these one should consult FindRoot.

Here is a nice way to use FindRoot.

f[x_] := Gamma[5 + x] - 40
g[x_] := (4 + x)! - 40

In:=
plot = Plot[f[x], {x, 0, 1}]
Print[StyleForm["the points used by the Plot function", FontColor ->
Blue]];
Short[points = Cases[plot, {(x_)?NumberQ, (y_)?NumberQ}, Infinity], 2]
Print[StyleForm["find where the function changes sign", FontColor ->
Blue]];
seeds = Position[Apply[Times, Partition[points[[All,2]], 2, 1], {1}],
x_ /; x <= 0]
Print[StyleForm["between this points in x axis there is a change in
sign of f[x]", FontColor -> Blue]];
samples = Extract[Partition[points[[All,1]], 2, 1], seeds]
Print[StyleForm["the root(s), at last!", FontColor -> Blue]];
(FindRoot[Gamma[5 + x] - 40, {x, #1[], #1[]}, WorkingPrecision -
> 40, PrecisionGoal -> 30] & ) /@ samples

plot = Plot[g[x], {x, 0, 1}]
Print[StyleForm["the points used by the Plot function", FontColor ->
Blue]];
Short[points = Cases[plot, {(x_)?NumberQ, (y_)?NumberQ}, Infinity], 2]
Print[StyleForm["find where the function changes sign", FontColor ->
Blue]];
seeds = Position[Apply[Times, Partition[points[[All, 2]], 2, 1],
{1}],  x_ /; x   0]
Print[StyleForm["between this points in x axis there is a change in
sign of
f[x]", FontColor -> Blue]];
samples = Extract[Partition[points[[All, 1]], 2, 1], seeds]
Print[StyleForm["the root(s), at last!", FontColor -> Blue]];
(FindRoot[Gamma[5 + x] - 40, {x, #1[], #1[]}, WorkingPrecision -
> 40, PrecisionGoal -> 30] &) /@ samples

BTW,

In:=
FunctionExpand[Gamma[5 + q]-40 == (4 + q)!-40]
Out=
True

/  Roger Bagula       :
> This type of solution comes up in algebra theory for symmetric type groups:
> Solve[(4+q)!-40==0,q]
> Solve[Gamma[5+q]-40==0,q]
> The Mathematica output refuses to give a number.
> I worked at it a little.
> Factorial:
> (4+q)!=40
> gives:
> 0.331291631797621
> Gamma[5+q]=40
> 0.3312924244499
> They stop agreeing at the 6th place in Mathematica.
> I don't have a good program for either of them and essentually did them
> by hand
> one digit at a time.
>
> Something is wrong with how I'm calculating them for sure.

```

• Prev by Date: Re: Iterating List
• Next by Date: Re: Re: String to List
• Previous by thread: Re: Iterating List
• Next by thread: Re: Re: String to List