home *** CD-ROM | disk | FTP | other *** search
/ Celestin Apprentice 2 / Apprentice-Release2.iso / Tools / Languages / RLaB 1.15c / examples / matrix_pow.r < prev    next >
Encoding:
Text File  |  1994-02-21  |  358 b   |  24 lines  |  [TEXT/RLAB]
  1. //
  2. // A simple examples of computing
  3. // A^p, where A is a matrix, and p
  4. // is an integer power.
  5. //
  6.  
  7. pow = function (a, p)
  8. {
  9.   local (mtmp1, mtmp2);
  10.  
  11.   if (p > 4)
  12.   {
  13.     mtmp1 = a*a;
  14.     for (i in 1:(p/2 - 2))
  15.     {
  16.       mtmp2 = mtmp1*a;
  17.       mtmp1 = mtmp2;
  18.     }
  19.     mtmp1 = mtmp1*mtmp1;
  20.     if (mod(p,2) != 0) { mtmp1 = mtmp1*a; }
  21.   }
  22.   return mtmp1;
  23. };
  24.