Correction: After fixing a bad filename.
- Demonstrative program about the reliability of the comment made by R. J. Cano at the entry A198960 on Jan 2012.
- 1) Donwload b198960.txt
- 2) Delete the offset column. For example in Unix/Linux systems with GNU file utils available try: user@host$ cut -d' ' -f2 b198960.txt > w.txt
- 3) Run this PARI/GP script. You never should be in the output the banner "The statement about a198960 is wrong!!". (Please first ensure to have allocated enough memory)
- v=vector(#w-1,k,w[k+1]); \\ Delete the first term a(1)=2;
- u=vector(#v,k,(v[k]-26)/6); \\ Convert into "y" parameter of (6*y+26 == 3*9^n-1) each case;
- s=vector(#u-1,k,u[k+1]-u[k]); \\ Take the first differences;
- print("Thirty-six is: ",s); \\ Consistency checking;
- r=vector(#s-1,k,s[k+1]/s[k]); \\ Compute the ratios between each pair of terms.
- r/=9; \\ PARI/GP does it like a human.
- if(#r==sum(j=1,#r,r[j]),print("Test approved for the given b-file."),print("The statement about a198960 is wrong!!"));
RAW Paste Data Copied