(*
The examples given in this notebook were generated by
running Mathematica Version 4.1 on a personal computer
(CPU: Pentium III, 800 MHz; OS: Windows ME).
*)

Needs["MatrixExp1`"]

(*
Matrix A is defined:
*)

Clear[t];

A={{   -2214499,      6961302,    -15057756,      23912112,     -27724410,
	         22991724,    -13243692,       5031936,    -1134747,      115150},
   {  -12376112,     98738712,   -361394496,     771497664,   -1059285024,
	        971579952,   -596081472,     236045376,   -54756144,     5668664},
   {  -99827109,    797360436,  -2842484292,    5931517032,   -7995875958,
	       7226909424,  -4381413876,    1718119512,  -395309997,    40641948},
   { -374179072,   2941573536, -10322399232,   21243450144,  -28288092672,
	      25293223872, -15189436416,    5906747232, -1349040384,   137790112},
   { -950140535,   7338687570, -25308482220,   51222262560,  -67135298706,
	      59143279860, -35033305020,   13453049520, -3037422735,   307006826},
   {-1857263184,  13986839592, -46979994816,   92548455552, -118044591840,
	     101234571984, -58420861632,   21879559296, -4823794512,   476716680},
   {-2903360377,  20906981712, -66753632916,  124253885784, -148871914254,
	     119272241592, -63982659972,   22173897096, -4504155489,   408353944},
   {-3489723104,  22568624496, -62326035072,   94573806432,  -83167495488,
	      38558638944,  -3332305536,   -5734135008,  2764938528,  -416677072},
   {-2377347147,   7903893006,   9821594628,  -95552651376,  223174371798,
	    -280042519428, 212802787764,  -98587002528, 25792860861, -2936350458},
   { 2594437520, -42075226440, 223387439040, -612148125120, 1009646243424,
	   -1066174315920, 730918587840, -316381688640, 78901848720, -8669563304}};
A=A t/362880;

(*
The matrix exponential of A is computed by using
M a t r i x E x p 1:
*)

Print["\nMatrixExp1[A][[1,1]]="]
MatrixExp1[A][[1,1]]
Print["\nCPU time used to evaluate MatrixExp1[A]:"]
Timing[MatrixExp1[A];][[1]]

(*
The matrix exponential of A is computed by using
M a t r i x E x p:
*)

Print["\nMatrixExp[A][[1,1]]="]
MatrixExp[A][[1,1]]
Print["\nCPU time used to evaluate MatrixExp[A]:"]
Timing[MatrixExp[A];][[1]]