Anonymous user
Pascal matrix generation: Difference between revisions
m
→{{header|REXX}}: added/changed comments and whitespace, used a template for the output section.
No edit summary |
m (→{{header|REXX}}: added/changed comments and whitespace, used a template for the output section.) |
||
Line 3,097:
numeric digits 50 /*be able to calculate huge factorials.*/
parse arg N . /*obtain the optional matrix size (N).*/
if N==''
call show N, upp(N), 'Pascal upper triangular matrix'
call show N, low(N), 'Pascal lower triangular matrix'
call show N, sym(N), 'Pascal symmetric matrix'
exit /*stick a fork in it, we're all done. */
/*──────────────────────────────────────────────────────────────────────────────────────*/
Line 3,112:
do i=0 for N; do j=0 for N; $=$ comb(i+j, i); end; end; return $
/*──────────────────────────────────────────────────────────────────────────────────────*/
!: procedure; parse arg x; !=1; do j=2 to x; != !*j; end;
/*──────────────────────────────────────────────────────────────────────────────────────*/
comb: procedure; parse arg x,y; if x=y then return 1 /* {=} case.*/
if y>x then return 0 /* {>} case.*/
if x-y<y then y= x-y; _= 1; do j=x-y+1 to x; _= _*j; end;
/*──────────────────────────────────────────────────────────────────────────────────────*/
show: procedure; parse arg s,@; w=0; #=0 /*get args. */
<pre>
──────────Pascal upper triangular matrix──────────
Line 3,158:
numeric digits 50 /*be able to calculate huge factorials.*/
parse arg N . /*obtain the optional matrix size (N).*/
if N==''
call show N, Pmat(N, 'upper'), 'Pascal upper triangular matrix'
call show N, Pmat(N, 'lower'), 'Pascal lower triangular matrix'
call show N, Pmat(N, 'sym') , 'Pascal symmetric matrix'
exit /*stick a fork in it, we're all done. */
/*──────────────────────────────────────────────────────────────────────────────────────*/
Pmat: procedure; parse arg N; $= /*generate a format of a Pascal matrix.*/
arg , ? /*get uppercase version of the 2nd arg.*/
return $
/*──────────────────────────────────────────────────────────────────────────────────────*/
!: procedure; parse arg x; !=1; do j=2 to x; != ! * j; end; return !
/*──────────────────────────────────────────────────────────────────────────────────────*/
comb: procedure; parse arg x,y; if x=y then return 1 /* {=} case.*/
if y>x then return 0 /* {>} case.*/
if x-y<y then y=x-y; _= 1; do j=x-y+1 to x; _= _ * j; end;
/*──────────────────────────────────────────────────────────────────────────────────────*/
show: procedure; parse arg s,@; w=0; #=0 /*get args. */
=={{header|Ring}}==
|