login  home  contents  what's new  discussion  bug reports help  links  subscribe  changes  refresh  edit

Below there is a function which generates Pascal triangle by iteratively filling rows of a matrix using recurence formula. Spad array indexing is one based, so we cheat and add 1 to indices.

fricas
fricas
)abbrev package PASCAL Pascal
Pascal : with ( make_pascal : Integer -> Matrix(Integer) )
make_pascal(N : Integer) : Matrix(Integer) ==
-- coerce here to NonNegativeInteger because this is what
-- 'zero' requires
nn := N::NonNegativeInteger
tab := zero(nn, nn)\$Matrix(Integer)
n : Integer
tab(1, 1) := 1
for n in 2..nn repeat
-- Fill row numer n
k : Integer
tab(n, 1) := 1
tab(n, n) := 1
for k in 2..(n - 1) repeat
tab(n, k) := tab(n-1, k) + tab(n-1, k-1)
fricas
Compiling FriCAS source code from file
using old system compiler.
PASCAL abbreviates package Pascal
------------------------------------------------------------------------
initializing NRLIB PASCAL for Pascal
compiling into NRLIB PASCAL
compiling exported make_pascal : Integer -> Matrix Integer
Time: 0.07 SEC.
(time taken in buildFunctor: 0)
;;; *** |Pascal| REDEFINED
;;; *** |Pascal| REDEFINED Time: 0.01 SEC.
Cumulative Statistics for Constructor Pascal Time: 0.08 seconds
finalizing NRLIB PASCAL Processing Pascal for Browser database: --->-->Pascal(constructor): Not documented!!!! --->-->Pascal((make_pascal ((Matrix (Integer)) (Integer)))): Not documented!!!! --->-->Pascal(): Missing Description ; compiling file "/var/aw/var/LatexWiki/PASCAL.NRLIB/PASCAL.lsp" (written 21 FEB 2023 12:45:59 PM):
; /var/aw/var/LatexWiki/PASCAL.NRLIB/PASCAL.fasl written ; compilation finished in 0:00:00.033 ------------------------------------------------------------------------ Pascal is now explicitly exposed in frame initial Pascal will be automatically loaded when needed from /var/aw/var/LatexWiki/PASCAL.NRLIB/PASCAL

Try it out:

fricas
make_pascal(5)
 (1)
Type: Matrix(Integer)