| 1 | pagelength = 1000; |
|---|
| 2 | proc sum |
|---|
| 3 | { |
|---|
| 4 | if (#[1] <= 1) { return(#[1]); } |
|---|
| 5 | return (#[1] + sum(#[1] -1)); |
|---|
| 6 | } |
|---|
| 7 | proc fak |
|---|
| 8 | { |
|---|
| 9 | if (#[1] <= 1) { return(1); } |
|---|
| 10 | return (#[1] * fak(#[1] -1)); |
|---|
| 11 | } |
|---|
| 12 | proc binom |
|---|
| 13 | { |
|---|
| 14 | return (fak(#[1]) / (fak(#[1]-#[2]) * fak(#[2]))); |
|---|
| 15 | } |
|---|
| 16 | "nesting too deep for fak(30)!!!!";fak(30); |
|---|
| 17 | int i; |
|---|
| 18 | int j; |
|---|
| 19 | //for (i = 2; i<=10; i=i+1) |
|---|
| 20 | //{ |
|---|
| 21 | // i, "sum", sum(i), "fak", fak(i); |
|---|
| 22 | //} |
|---|
| 23 | //"binom(5,1)", binom(5,1); |
|---|
| 24 | //$; |
|---|
| 25 | string b; |
|---|
| 26 | for (i=1; i<=16; i=i+1) |
|---|
| 27 | { |
|---|
| 28 | b = "binoms " + string(i) + " "; |
|---|
| 29 | for (j=0; j<=i; j=j+1) |
|---|
| 30 | { |
|---|
| 31 | b = b + " " + string(binom(i,j)); |
|---|
| 32 | } |
|---|
| 33 | b; |
|---|
| 34 | } |
|---|
| 35 | LIB "tst.lib";tst_status(1);$; |
|---|
| 36 | "binom(5,0)", binom(5,0); |
|---|
| 37 | "binom(5,1)", binom(5,1); |
|---|
| 38 | "binom(5,2)", binom(5,2); |
|---|
| 39 | "binom(5,3)", binom(5,3); |
|---|
| 40 | "binom(5,4)", binom(5,4); |
|---|
| 41 | "binom(5,5)", binom(5,5); |
|---|
| 42 | LIB "tst.lib";tst_status(1);$; |
|---|
