cancel
Showing results for 
Search instead for 
Did you mean: 
cancel
Showing results for 
Search instead for 
Did you mean: 

How to correct the Pell(n) function ? (for n = 181)

SOLVED
Alexandrite

How to correct the Pell(n) function ? (for n = 181)

Hello Everyone,

Pell%28181%29.PNG

How to correct the Pell(n) function ? ( for n = 181)

Thanks in advance for your time and help.

Regards.

1 ACCEPTED SOLUTION

Accepted Solutions

Re: How to correct the Pell(n) function ? (for n = 181)

For some N (like 2,5,10) m is just 1 (second element of C is already greater than the root of N) and so your calculation of the rearranged vector c fails. Either include another if-statement to cope with that case (then c<--C) or do it completely different like in the attached.

1.png

0.png

View solution in original post

14 REPLIES 14

Re: How to correct the Pell(n) function ? (for n = 181)

Looking at Numerical Recipes, I think you need to add a correction factor in for when w[0,1 is zero.

Stuart

Highlighted

Re: How to correct the Pell(n) function ? (for n = 181)

I greatly appreciate your help, Stuart.

With :

Pell

http://www.jakebakermaths.org.uk/maths/jshtmlpellsolverbigintegerv10.html

(2).PNG

Best Regards.

Re: How to correct the Pell(n) function ? (for n = 181)

Unfortunately, I'm Mathcadless for a week or so, but you can define your function as follows.

Rename my function to PellA.

Define

Pell(n):=p <-PellA(n, infinity)

stack (p[1, p[2)

thenPell(P181).

Stuart

Re: How to correct the Pell(n) function ? (for n = 181)

Hello Everyone again,

1.PNG

2.PNG

Now I need help with correcting the Pell_Test(N) function, the above. ("This array index is invalid for this array")

Thanks in advance for your time and help.

Best Regards.

Re: How to correct the Pell(n) function ? (for n = 181)

For some N (like 2,5,10) m is just 1 (second element of C is already greater than the root of N) and so your calculation of the rearranged vector c fails. Either include another if-statement to cope with that case (then c<--C) or do it completely different like in the attached.

1.png

0.png

View solution in original post

Re: How to correct the Pell(n) function ? (for n = 181)

I greatly appreciate your help. The answer is very neat . Many, many thanks, Werner.

3.PNG

By the way, I need help with how to type "the thin equal sign" (the green circle) in Mathcad 14.

Best Regards.

Message was edited by: Loi Le

Re: How to correct the Pell(n) function ? (for n = 181)

By the way, I need help with how to type "the thin equal sign" (the green circle) in Mathcad 14.

Same way as the * and the +. If you create a math region and you type Ctrl-Shift-K (cursor gets red as in a text region) you can use literally any character in a variable name. Pressing Ctrl-Shift-K again brings you back to normal mode. That way you can use characters in variable names which normally would be taken as operators. The sam way you can also insert some spaces instead of using the underline character as you had done.

Re: How to correct the Pell(n) function ? (for n = 181)

I got it. . Many, many, many thanks, Werner. .

Best Regards.

Re: How to correct the Pell(n) function ? (for n = 181)

I guess the Pell_Test(N) function need to adjust a little bit :

4.PNG

after adjusting a little bit :

5.PNG

And while I am learning about match( ) , submatrix( ) function, I adjust the function a little bit :

3.PNG

Thanks again.

Best Regards.

Message was edited by: Loi Le