Anonymous user
Wilson primes of order n: Difference between revisions
m
→{{header|REXX}}: used an IF/THEN/ELSE statement, changed some comments.
m (→{{header|REXX}}: changed a comment.) |
m (→{{header|REXX}}: used an IF/THEN/ELSE statement, changed some comments.) |
||
Line 247:
<lang rexx>/*REXX program finds and displays Wilson primes: a prime P such that P**2 divides:*/
/*────────────────── (n-1)! * (P-n)! - (-1)**n where n is 1 ──◄ 11, and P < 18.*/
parse arg oLO oHI hip . /*obtain optional
if oLO=='' | oLO=="," then oLO= 1 /*Not specified? Then use the default.*/
if oHI=='' | oHI=="," then oHI= 11 /* " " " " " " */
Line 256:
parse value bignum 'E0' with ex 'E' ex . /*obtain possible exponent of factorial*/
numeric digits (max(9, ex+2) ) /*calculate max # of dec. digits needed*/
call facts hip /*
title= ' Wilson primes P of order ' oLO " ──► " oHI', where P < ' commas(hip)
w= length(title) + 1 /*width of columns of possible numbers.*/
say ' order │'center(title, w )
say '───────┼'center("" , w, '─')
do n=oLO to oHI; pom=
nmf= !(n-1)
if n==1 then lim= 103 /*limit to known primes for 1st order. */
else lim= #
$=
$= $ ' ' commas(p) /*add a commatized prime ──► $ list.*/
end /*p*/
if $=='' then $= ' (none found within the range specified)'
say center(n, 7)'│' substr($, 2)
end /*n*/
say '───────┴'center("" , w, '─')
exit 0 /*stick a fork in it, we're all done. */
|