This is a couple of classes written in C++ that handles
Bignums (integers whose size is limited only by available memory)
and matrices.
The matrix directory contains several versions of the code, use whatever suits
you best. (There is even a C version includd that was my first attempt this)
The PMatrix class is a permutation matrix and the LUMatrix is a matrix that
has been LU-decomposed. (if you didn't understand the last sentence so don't worry,
you probably won't need them then :-) )
The documentation is sorely lacking but using the classes is fairly 
straight-forward. (The Bignum class is fairly well commented - in Swedish only. :-( )

This code is hereby placed into the public domain.

 PERMSSN    UID  GID    PACKED    SIZE  RATIO     CRC       STAMP          NAME
---------- ----------- ------- ------- ------ ---------- ------------ -------------
drwx--s---  3533/39          0       0 ****** -lhd- 0000 Jan 11  1997 math_classes/
-rw-------  3533/39        472     851  55.5% -lh5- ebf4 Jan 11  1997 math_classes/ReadMe
drwx--s---  3533/39          0       0 ****** -lhd- 0000 Jan 11  1997 math_classes/bignum/
-rw-------  3533/39       4309   18536  23.2% -lh5- a655 Apr 14  1996 math_classes/bignum/bignum.cc
-rw-------  3533/39        288    1113  25.9% -lh5- 6b92 Mar 19  1995 math_classes/bignum/bignum.h
-rw-------  3533/39        383    1919  20.0% -lh5- 5a34 Mar 19  1995 math_classes/bignum/bignumtest.cc
-rw-------  3533/39        121     268  45.1% -lh5- 9dae Mar 24  1996 math_classes/bignum/makefile
drwx--s---  3533/39          0       0 ****** -lhd- 0000 Jan 11  1997 math_classes/matrix/
-rw-------  3533/39       2870   14206  20.2% -lh5- 2d49 Sep 17  1996 math_classes/matrix/Matrix.cc
-rw-------  3533/39        826    3304  25.0% -lh5- 3e18 Jun 10  1996 math_classes/matrix/Vector.cc
-rw-------  3533/39       5350   10612  50.4% -lh5- 4737 Oct 27  1996 math_classes/matrix/crank
-rw-------  3533/39        333     647  51.5% -lh5- 7fee Oct 27  1996 math_classes/matrix/crank.cc
-rw-------  3533/39        594    1267  46.9% -lh5- 94ac Sep 19  1996 math_classes/matrix/jacobi.cc
drwx--s---  3533/39          0       0 ****** -lhd- 0000 Jan 11  1997 math_classes/matrix/libmatrix/
-rw-------  3533/39        412    1082  38.1% -lh5- 0b57 Oct  5  1996 math_classes/matrix/libmatrix/Makefile
-rw-------  3533/39        695    2698  25.8% -lh5- 347e Oct  5  1996 math_classes/matrix/libmatrix/Matrix.h
-rw-------  3533/39        350     694  50.4% -lh5- ac5e Sep 18  1996 math_classes/matrix/libmatrix/determinant.cc
-rw-------  3533/39      36181  189904  19.1% -lh5- f8b5 Oct  5  1996 math_classes/matrix/libmatrix/libmatrix.a.020881
-rw-------  3533/39        153     244  62.7% -lh5- 920a Sep 18  1996 math_classes/matrix/libmatrix/lumat_getl.cc
-rw-------  3533/39         78      86  90.7% -lh5- 9b36 Sep 18  1996 math_classes/matrix/libmatrix/lumat_getp.cc
-rw-------  3533/39        139     186  74.7% -lh5- d898 Sep 18  1996 math_classes/matrix/libmatrix/lumat_getu.cc
-rw-------  3533/39        414    1038  39.9% -lh5- 124c Sep 18  1996 math_classes/matrix/libmatrix/lumat_ludecomp.cc
-rw-------  3533/39         92     110  83.6% -lh5- ec29 Sep 18  1996 math_classes/matrix/libmatrix/lumat_op_assign.cc
-rw-------  3533/39        348     774  45.0% -lh5- a1bc Oct  5  1996 math_classes/matrix/libmatrix/lusolve.cc
-rw-------  3533/39         62      66  93.9% -lh5- d498 Sep 18  1996 math_classes/matrix/libmatrix/mat_create.cc
-rw-------  3533/39        166     248  66.9% -lh5- 34fc Sep 18  1996 math_classes/matrix/libmatrix/mat_create_int_int.cc
-rw-------  3533/39        207     331  62.5% -lh5- 018c Sep 18  1996 math_classes/matrix/libmatrix/mat_create_int_int_ptr.cc
-rw-------  3533/39        174     261  66.7% -lh5- d423 Oct  5  1996 math_classes/matrix/libmatrix/mat_create_mat.cc
-rw-------  3533/39        184     272  67.6% -lh5- 40ec Sep 18  1996 math_classes/matrix/libmatrix/mat_create_pmat.cc
-rw-------  3533/39         51      51 100.0% -lh0- ce5d Sep 19  1996 math_classes/matrix/libmatrix/mat_dabs.cc
-rw-------  3533/39         73      74  98.6% -lh5- f806 Sep 18  1996 math_classes/matrix/libmatrix/mat_destr.cc
-rw-------  3533/39        127     175  72.6% -lh5- 523b Oct  5  1996 math_classes/matrix/libmatrix/mat_getcol.cc
-rw-------  3533/39        127     175  72.6% -lh5- cadd Oct  5  1996 math_classes/matrix/libmatrix/mat_getrow.cc
-rw-------  3533/39        159     234  67.9% -lh5- c1fc Oct  5  1996 math_classes/matrix/libmatrix/mat_inverse.cc
-rw-------  3533/39        198     343  57.7% -lh5- ef69 Sep 18  1996 math_classes/matrix/libmatrix/mat_op_add.cc
-rw-------  3533/39        178     280  63.6% -lh5- 0fec Sep 18  1996 math_classes/matrix/libmatrix/mat_op_addassign.cc
-rw-------  3533/39        244     565  43.2% -lh5- 135f Sep 18  1996 math_classes/matrix/libmatrix/mat_op_assign.cc
-rw-------  3533/39        171     302  56.6% -lh5- 0c06 Sep 18  1996 math_classes/matrix/libmatrix/mat_op_eq.cc
-rw-------  3533/39        147     205  71.7% -lh5- 0d15 Sep 18  1996 math_classes/matrix/libmatrix/mat_op_mult_dbl_mat.cc
-rw-------  3533/39        146     204  71.6% -lh5- bfae Sep 18  1996 math_classes/matrix/libmatrix/mat_op_mult_mat_dbl.cc
-rw-------  3533/39        234     391  59.8% -lh5- 2200 Sep 18  1996 math_classes/matrix/libmatrix/mat_op_mult_mat_mat.cc
-rw-------  3533/39        185     283  65.4% -lh5- e976 Sep 18  1996 math_classes/matrix/libmatrix/mat_op_mult_pmat_mat.cc
-rw-------  3533/39        134     194  69.1% -lh5- 5a05 Sep 18  1996 math_classes/matrix/libmatrix/mat_op_multassign.cc
-rw-------  3533/39        168     300  56.0% -lh5- 5701 Sep 18  1996 math_classes/matrix/libmatrix/mat_op_neq.cc
-rw-------  3533/39        199     343  58.0% -lh5- d279 Sep 18  1996 math_classes/matrix/libmatrix/mat_op_sub.cc
-rw-------  3533/39        178     279  63.8% -lh5- 7ac5 Sep 18  1996 math_classes/matrix/libmatrix/mat_op_subassign.cc
-rw-------  3533/39        294     516  57.0% -lh5- 86d2 Oct  5  1996 math_classes/matrix/libmatrix/mat_print.cc
-rw-------  3533/39        181     293  61.8% -lh5- 1716 Sep 18  1996 math_classes/matrix/libmatrix/mat_rxchg.cc
-rw-------  3533/39        147     193  76.2% -lh5- 76e1 Oct  5  1996 math_classes/matrix/libmatrix/mat_setcol.cc
-rw-------  3533/39        145     192  75.5% -lh5- b2f5 Oct  5  1996 math_classes/matrix/libmatrix/mat_setrow.cc
-rw-------  3533/39        176     293  60.1% -lh5- f617 Sep 18  1996 math_classes/matrix/libmatrix/mat_setunity.cc
-rw-------  3533/39        175     263  66.5% -lh5- 408b Sep 18  1996 math_classes/matrix/libmatrix/mat_transpose.cc
-rw-------  3533/39        125     165  75.8% -lh5- 4531 Sep 18  1996 math_classes/matrix/libmatrix/mat_unity.cc
-rw-------  3533/39         61      61 100.0% -lh0- f573 Sep 18  1996 math_classes/matrix/libmatrix/pmat_create.cc
-rw-------  3533/39        145     206  70.4% -lh5- 8bb9 Sep 18  1996 math_classes/matrix/libmatrix/pmat_create_int.cc
-rw-------  3533/39        166     257  64.6% -lh5- b5bd Sep 18  1996 math_classes/matrix/libmatrix/pmat_create_pmat.cc
-rw-------  3533/39         81      86  94.2% -lh5- 64e9 Sep 18  1996 math_classes/matrix/libmatrix/pmat_destr.cc
-rw-------  3533/39        122     156  78.2% -lh5- 0c9e Sep 18  1996 math_classes/matrix/libmatrix/pmat_inverse.cc
-rw-------  3533/39        142     167  85.0% -lh5- da50 Sep 18  1996 math_classes/matrix/libmatrix/pmat_op_arr.cc
-rw-------  3533/39        154     219  70.3% -lh5- 091b Sep 18  1996 math_classes/matrix/libmatrix/pmat_op_assign.cc
-rw-------  3533/39        181     238  76.1% -lh5- c7c0 Sep 18  1996 math_classes/matrix/libmatrix/pmat_op_index.cc
-rw-------  3533/39        165     240  68.8% -lh5- 0f06 Sep 18  1996 math_classes/matrix/libmatrix/pmat_op_mult_pmat_pmat.cc
-rw-------  3533/39        124     158  78.5% -lh5- 3812 Sep 18  1996 math_classes/matrix/libmatrix/pmat_transpose.cc
-rw-------  3533/39         84     100  84.0% -lh5- 9c95 Oct  5  1996 math_classes/matrix/libmatrix/solve_lumat.cc
-rw-------  3533/39         85     108  78.7% -lh5- 5da7 Oct  5  1996 math_classes/matrix/libmatrix/solve_mat.cc
-rw-------  3533/39        128     198  64.6% -lh5- b527 Sep 18  1996 math_classes/matrix/libmatrix/swap_mat_mat.cc
-rw-------  3533/39        248     526  47.1% -lh5- 4dd3 Sep 22  1996 math_classes/matrix/matrixtest.cc
drwx--s---  3533/39          0       0 ****** -lhd- 0000 Jan 11  1997 math_classes/matrix/old/
-rw-------  3533/39       1780    7105  25.1% -lh5- 15ff Jun  4  1996 math_classes/matrix/old/cMatrix.c
-rw-------  3533/39        337     977  34.5% -lh5- 6bb3 Jun  4  1996 math_classes/matrix/old/cMatrix.h
-rw-------  3533/39        453     970  46.7% -lh5- a643 Jun  4  1996 math_classes/matrix/old/cmtest.c
-rw-------  3533/39       2694   12580  21.4% -lh5- 9bd7 Feb 19  1995 math_classes/matrix/old/oldmatrix.cc
-rw-------  3533/39        825    3316  24.9% -lh5- 912e May 21  1995 math_classes/matrix/old/oldmatrix.h
-rw-------  3533/39        857    3448  24.9% -lh5- 8f74 Feb 12  1996 math_classes/matrix/old/oldvector.cc
-rw-------  3533/39        359    1339  26.8% -lh5- d16f May 21  1995 math_classes/matrix/old/oldvector.h
-rw-------  3533/39        229     523  43.8% -lh5- 64b5 Oct  5  1996 math_classes/matrix/test.cc
-rw-------  3533/39       2656    5928  44.8% -lh5- 68ba Oct  5  1996 math_classes/matrix/test.o
-rw-------  3533/39        176     280  62.9% -lh5- 2f3d Sep 19  1996 math_classes/matrix/vectest.cc
---------- ----------- ------- ------- ------ ---------- ------------ -------------
 Total        78 files   71415  296216  24.1%            Jan 16  1997