Wrote

Thanks, will explore.

Wrote

Try the other one for gamma calculation.

Wrote

Wrote

Try the other one for gamma calculation.

round(■ ,■ ,■ ) ... undefined
H.T.Davis [Abramowitz & Stegun] -> very objective.
In the mean time I'm finishing two superb applications examples
... examples based on your first Gamma(a,x) version,
Thanks for that one, gorgeous ... Jean.

Wrote

As offered ... two applications.
Now, the old Rooster is going in the bed marmite ... Jean.

gamma(a,x) Applications.sm (23 КиБ) скачан 116 раз(а).

Wrote

I suggest you to change gamma function with gamma_r4.
It is 3 times faster and much more precise.

Wrote

Thanks for gamma_r4
1. NO gain timing both applications
2. gamma_r4 Does NOT solve ... first version solves.

Wrote

What do you mean by 'does not solve'?

Wrote

Wrote

What do you mean by 'does not solve'?

Wrote

1. NO gain timing both applications

Wrote

Wrote

What do you mean by 'does not solve'?

Wrote

It solves, not with conventional methods.
They tend to give errors, not gamma_r4's fault.
alglib is the answer, it usually gives result.

Wrote

Good rescue ... OK.
As it looks, sr4 resides at the kernel level.
Thus, it plots but not solve(,,,,)
because not scalar wrt 'x' for the solve bloc.

Wrote

Check the sample below

Wrote

I know I have flooded this topic too much

Wrote

Wrote

I know I have flooded this topic too much

You may not like my verdicts:
1. on-line faster has no or little interest, neither up ^307.
2. Your first version runs fine the two applications.
3. The champion is the long time ago published Alvaro Γ(x)
4. for these two applications [2] & [3] are same
but => [3] drops timing [2] from 24 s down 18 s
The drop in timing results from Alvaro Γ(x)
running at the kernel scalar level.
BTW, my original H.T. Davis sanity Mathcad & Alvaro Γ(x)
By same token, thanks Alvaro for your Γ(x).
Cheers ... Jean

Maths Special Gamma(a,x) Incomplete APPLICATIONS [Alvaro G(x)].sm (26 КиБ) скачан 114 раз(а).

Показать спойлер