Hex Artifact Content
Not logged in

Artifact 461a3190dc2984b63f1ba4ecb12ce84daa0e8482:


0000: 0d 0a 0d 0a 2f 2f 20 20 6d 61 74 72 69 78 20 61  ....//  matrix a
0010: 6e 64 20 76 65 63 74 6f 72 20 66 75 6e 63 74 69  nd vector functi
0020: 6f 6e 73 2e 0d 0a 2f 2f 0d 0a 2f 2f 0d 0a 0d 0a  ons...//..//....
0030: 66 75 6e 63 74 69 6f 6e 20 76 65 63 74 6f 72 5f  function vector_
0040: 61 64 64 28 61 2c 20 62 2c 20 72 65 73 75 6c 74  add(a, b, result
0050: 29 7b 0d 0a 20 20 20 20 69 66 28 21 72 65 73 75  ){..    if(!resu
0060: 6c 74 29 20 72 65 73 75 6c 74 20 3d 20 5b 5d 3b  lt) result = [];
0070: 0d 0a 09 69 66 28 21 61 29 20 74 68 72 6f 77 20  ...if(!a) throw 
0080: 22 66 69 72 73 74 20 61 72 67 75 6d 65 6e 74 20  "first argument 
0090: 74 6f 20 76 65 63 74 6f 72 5f 61 64 64 20 6d 69  to vector_add mi
00a0: 73 73 69 6e 67 22 3b 0d 0a 09 69 66 28 21 62 29  ssing";...if(!b)
00b0: 20 74 68 72 6f 77 20 22 73 65 63 6f 6e 64 20 61   throw "second a
00c0: 72 67 75 6d 65 6e 74 20 74 6f 20 76 65 63 74 6f  rgument to vecto
00d0: 72 5f 61 64 64 20 6d 69 73 73 69 6e 67 22 3b 0d  r_add missing";.
00e0: 0a 0d 0a 20 20 20 20 69 66 28 61 2e 6c 65 6e 67  ...    if(a.leng
00f0: 74 68 20 3e 20 62 2e 6c 65 6d 67 74 68 29 7b 0d  th > b.lemgth){.
0100: 0a 20 20 20 20 20 20 20 20 76 61 72 20 74 6d 70  .        var tmp
0110: 20 3d 20 61 3b 0d 0a 20 20 20 20 20 20 20 20 61   = a;..        a
0120: 20 3d 20 62 3b 20 62 20 3d 20 74 6d 70 3b 20 7d   = b; b = tmp; }
0130: 0d 0a 0d 0a 20 20 20 20 66 6f 72 28 76 61 72 20  ....    for(var 
0140: 69 20 3d 20 30 3b 20 69 20 3c 20 61 2e 6c 65 6e  i = 0; i < a.len
0150: 67 74 68 3b 20 2b 2b 69 29 7b 0d 0a 20 20 20 20  gth; ++i){..    
0160: 20 20 20 20 72 65 73 75 6c 74 5b 69 5d 20 3d 20      result[i] = 
0170: 61 5b 69 5d 20 2b 20 62 5b 69 5d 3b 20 7d 0d 0a  a[i] + b[i]; }..
0180: 0d 0a 20 20 20 20 66 6f 72 28 76 61 72 20 69 20  ..    for(var i 
0190: 3d 20 61 2e 6c 65 6e 67 74 68 3b 20 69 20 3c 20  = a.length; i < 
01a0: 62 2e 6c 65 6e 67 74 68 3b 20 2b 2b 69 29 7b 0d  b.length; ++i){.
01b0: 0a 20 20 20 20 20 20 20 20 72 65 73 75 6c 74 5b  .        result[
01c0: 69 5d 20 3d 20 62 5b 69 5d 3b 20 7d 0d 0a 0d 0a  i] = b[i]; }....
01d0: 20 20 20 20 72 65 74 75 72 6e 20 72 65 73 75 6c      return resul
01e0: 74 3b 20 7d 0d 0a 0d 0a 20 20 20 20 0d 0a 66 75  t; }....    ..fu
01f0: 6e 63 74 69 6f 6e 20 76 65 63 74 6f 72 5f 63 70  nction vector_cp
0200: 79 28 72 65 73 75 6c 74 2c 20 61 29 7b 0d 0a 09  y(result, a){...
0210: 69 66 28 21 61 29 20 74 68 72 6f 77 20 22 66 69  if(!a) throw "fi
0220: 72 73 74 20 61 72 67 75 6d 65 6e 74 20 74 6f 20  rst argument to 
0230: 76 65 63 74 6f 72 5f 63 70 79 20 6d 69 73 73 69  vector_cpy missi
0240: 6e 67 22 3b 0d 0a 20 20 20 20 66 6f 72 28 76 61  ng";..    for(va
0250: 72 20 69 20 3d 20 30 3b 20 69 20 3c 20 61 2e 6c  r i = 0; i < a.l
0260: 65 6e 67 74 68 3b 20 2b 2b 69 29 7b 0d 0a 20 20  ength; ++i){..  
0270: 20 20 20 20 20 20 72 65 73 75 6c 74 5b 69 5d 20        result[i] 
0280: 3d 20 61 5b 69 5d 3b 20 7d 0d 0a 7d 0d 0a 0d 0a  = a[i]; }..}....
0290: 0d 0a 0d 0a 66 75 6e 63 74 69 6f 6e 20 76 65 63  ....function vec
02a0: 74 6f 72 5f 6d 69 6e 75 73 28 61 2c 20 62 2c 20  tor_minus(a, b, 
02b0: 72 65 73 75 6c 74 29 7b 0d 0a 20 20 20 20 69 66  result){..    if
02c0: 28 21 72 65 73 75 6c 74 29 20 72 65 73 75 6c 74  (!result) result
02d0: 20 3d 20 5b 5d 3b 0d 0a 09 69 66 28 21 61 29 20   = [];...if(!a) 
02e0: 74 68 72 6f 77 20 22 66 69 72 73 74 20 61 72 67  throw "first arg
02f0: 75 6d 65 6e 74 20 74 6f 20 76 65 63 74 6f 72 5f  ument to vector_
0300: 6d 69 6e 75 73 20 6d 69 73 73 69 6e 67 22 3b 0d  minus missing";.
0310: 0a 09 69 66 28 21 62 29 20 74 68 72 6f 77 20 22  ..if(!b) throw "
0320: 73 65 63 6f 6e 64 20 61 72 67 75 6d 65 6e 74 20  second argument 
0330: 74 6f 20 76 65 63 74 6f 72 5f 6d 69 6e 75 73 20  to vector_minus 
0340: 6d 69 73 73 69 6e 67 22 3b 0d 0a 0d 0a 20 20 20  missing";....   
0350: 20 76 61 72 20 6c 69 6d 20 3d 20 4d 61 74 68 2e   var lim = Math.
0360: 6d 61 78 28 61 2e 6c 65 6e 67 74 68 2c 20 62 2e  max(a.length, b.
0370: 6c 65 6e 67 74 68 29 3b 0d 0a 20 20 20 20 0d 0a  length);..    ..
0380: 20 20 20 20 66 6f 72 28 76 61 72 20 69 20 3d 20      for(var i = 
0390: 30 3b 20 69 20 3c 20 6c 69 6d 3b 20 2b 2b 69 29  0; i < lim; ++i)
03a0: 7b 0d 0a 09 09 76 61 72 20 61 56 61 6c 75 65 2c  {....var aValue,
03b0: 20 62 56 61 6c 75 65 3b 0d 0a 09 09 69 66 28 69   bValue;....if(i
03c0: 20 3c 20 61 2e 6c 65 6e 67 74 68 29 0d 0a 09 09   < a.length)....
03d0: 09 61 56 61 6c 75 65 20 3d 20 61 5b 69 5d 3b 0d  .aValue = a[i];.
03e0: 0a 09 09 65 6c 73 65 0d 0a 09 09 09 61 56 61 6c  ...else.....aVal
03f0: 75 65 20 3d 20 30 3b 0d 0a 0d 0a 09 09 69 66 28  ue = 0;......if(
0400: 69 20 3c 20 62 2e 6c 65 6e 67 74 68 29 0d 0a 09  i < b.length)...
0410: 09 09 62 56 61 6c 75 65 20 3d 20 62 5b 69 5d 3b  ..bValue = b[i];
0420: 0d 0a 09 09 65 6c 73 65 0d 0a 09 09 09 62 56 61  ....else.....bVa
0430: 6c 75 65 20 3d 20 30 3b 0d 0a 09 09 09 0d 0a 09  lue = 0;........
0440: 09 72 65 73 75 6c 74 5b 69 5d 20 3d 20 61 56 61  .result[i] = aVa
0450: 6c 75 65 20 2d 20 62 56 61 6c 75 65 3b 20 7d 0d  lue - bValue; }.
0460: 0a 20 20 20 20 20 20 20 20 0d 0a 20 20 20 20 20  .        ..     
0470: 20 20 20 0d 0a 20 20 20 20 72 65 74 75 72 6e 20     ..    return 
0480: 72 65 73 75 6c 74 3b 20 7d 0d 0a 20 20 20 0d 0a  result; }..   ..
0490: 0d 0a 0d 0a 66 75 6e 63 74 69 6f 6e 20 76 65 63  ....function vec
04a0: 74 6f 72 5f 64 6f 74 28 61 2c 20 62 29 7b 0d 0a  tor_dot(a, b){..
04b0: 20 20 20 20 76 61 72 20 6c 69 6d 20 3d 20 4d 61      var lim = Ma
04c0: 74 68 2e 6d 69 6e 28 61 2e 6c 65 6e 67 74 68 2c  th.min(a.length,
04d0: 20 62 2e 6c 65 6e 67 74 68 29 3b 0d 0a 09 69 66   b.length);...if
04e0: 28 21 61 29 20 74 68 72 6f 77 20 22 66 69 72 73  (!a) throw "firs
04f0: 74 20 61 72 67 75 6d 65 6e 74 20 74 6f 20 76 65  t argument to ve
0500: 63 74 6f 72 5f 64 6f 74 20 6d 69 73 73 69 6e 67  ctor_dot missing
0510: 22 3b 0d 0a 09 69 66 28 21 62 29 20 74 68 72 6f  ";...if(!b) thro
0520: 77 20 22 73 65 63 6f 6e 64 20 61 72 67 75 6d 65  w "second argume
0530: 6e 74 20 74 6f 20 76 65 63 74 6f 72 5f 64 6f 74  nt to vector_dot
0540: 20 6d 69 73 73 69 6e 67 22 3b 0d 0a 09 0d 0a 20   missing";..... 
0550: 20 20 20 76 61 72 20 72 65 73 75 6c 74 20 3d 20     var result = 
0560: 30 3b 0d 0a 0d 0a 20 20 20 20 66 6f 72 28 76 61  0;....    for(va
0570: 72 20 69 20 3d 20 30 3b 20 69 20 3c 20 6c 69 6d  r i = 0; i < lim
0580: 3b 20 2b 2b 69 29 7b 0d 0a 20 20 20 20 20 20 20  ; ++i){..       
0590: 20 72 65 73 75 6c 74 20 2b 3d 20 61 5b 69 5d 20   result += a[i] 
05a0: 2a 20 62 5b 69 5d 3b 20 7d 0d 0a 0d 0a 20 20 20  * b[i]; }....   
05b0: 20 72 65 74 75 72 6e 20 72 65 73 75 6c 74 3b 20   return result; 
05c0: 7d 0d 0a 20 20 20 20 0d 0a 20 20 20 20 0d 0a 66  }..    ..    ..f
05d0: 75 6e 63 74 69 6f 6e 20 76 65 63 74 6f 72 5f 6e  unction vector_n
05e0: 6f 72 6d 28 61 29 7b 0d 0a 09 69 66 28 21 61 29  orm(a){...if(!a)
05f0: 20 74 68 72 6f 77 20 22 61 72 67 75 6d 65 6e 74   throw "argument
0600: 20 74 6f 20 76 65 63 74 6f 72 5f 6e 6f 72 6d 20   to vector_norm 
0610: 6d 69 73 73 69 6e 67 22 3b 0d 0a 20 20 20 20 72  missing";..    r
0620: 65 74 75 72 6e 20 4d 61 74 68 2e 73 71 72 74 28  eturn Math.sqrt(
0630: 76 65 63 74 6f 72 5f 64 6f 74 28 61 2c 61 29 29  vector_dot(a,a))
0640: 3b 20 7d 0d 0a 0d 0a 20 20 20 20 0d 0a 20 20 20  ; }....    ..   
0650: 20 0d 0a 20 20 20 20 0d 0a 66 75 6e 63 74 69 6f   ..    ..functio
0660: 6e 20 76 65 63 74 6f 72 5f 73 63 61 6c 65 28 61  n vector_scale(a
0670: 2c 20 73 2c 20 72 65 73 75 6c 74 29 7b 0d 0a 09  , s, result){...
0680: 69 66 28 21 61 29 20 74 68 72 6f 77 20 22 66 69  if(!a) throw "fi
0690: 72 73 74 20 61 72 67 75 6d 65 6e 74 20 74 6f 20  rst argument to 
06a0: 76 65 63 74 6f 72 5f 73 63 61 6c 65 20 6d 69 73  vector_scale mis
06b0: 73 69 6e 67 22 3b 0d 0a 09 69 66 28 73 20 3d 3d  sing";...if(s ==
06c0: 3d 20 75 6e 64 65 66 69 6e 65 64 20 7c 7c 20 73  = undefined || s
06d0: 20 3d 3d 3d 20 6e 75 6c 6c 29 20 74 68 72 6f 77   === null) throw
06e0: 20 22 73 63 61 6c 61 72 20 61 72 67 75 6d 65 6e   "scalar argumen
06f0: 74 20 74 6f 20 76 65 63 74 6f 72 5f 73 63 61 6c  t to vector_scal
0700: 65 20 6d 69 73 73 69 6e 67 22 3b 0d 0a 0d 0a 20  e missing";.... 
0710: 20 20 20 69 66 28 21 72 65 73 75 6c 74 29 20 72     if(!result) r
0720: 65 73 75 6c 74 20 3d 20 5b 5d 3b 0d 0a 20 20 20  esult = [];..   
0730: 20 0d 0a 20 20 20 20 66 6f 72 28 76 61 72 20 69   ..    for(var i
0740: 20 3d 20 30 3b 20 69 20 3c 20 61 2e 6c 65 6e 67   = 0; i < a.leng
0750: 74 68 3b 20 2b 2b 69 29 7b 0d 0a 20 20 20 20 20  th; ++i){..     
0760: 20 20 20 72 65 73 75 6c 74 5b 69 5d 20 3d 20 73     result[i] = s
0770: 20 2a 20 61 5b 69 5d 3b 20 7d 0d 0a 0d 0a 20 20   * a[i]; }....  
0780: 20 20 72 65 74 75 72 6e 20 72 65 73 75 6c 74 3b    return result;
0790: 20 7d 0d 0a 0d 0a 0d 0a 0d 0a 2f 2f 20 72 65 74   }........// ret
07a0: 75 72 6e 73 20 74 72 75 65 20 69 66 66 20 74 68  urns true iff th
07b0: 65 20 7a 65 72 6f 20 74 61 69 6c 20 70 61 64 64  e zero tail padd
07c0: 65 64 20 76 65 63 74 6f 72 73 20 6d 61 74 63 68  ed vectors match
07d0: 0d 0a 66 75 6e 63 74 69 6f 6e 20 76 65 63 74 6f  ..function vecto
07e0: 72 5f 63 6d 70 28 61 2c 20 62 29 7b 0d 0a 09 69  r_cmp(a, b){...i
07f0: 66 28 21 61 29 20 74 68 72 6f 77 20 22 66 69 72  f(!a) throw "fir
0800: 73 74 20 61 72 67 75 6d 65 6e 74 20 74 6f 20 76  st argument to v
0810: 65 63 74 6f 72 5f 63 6d 70 20 6d 69 73 73 69 6e  ector_cmp missin
0820: 67 22 3b 0d 0a 09 69 66 28 21 62 29 20 74 68 72  g";...if(!b) thr
0830: 6f 77 20 22 73 65 63 6f 6e 64 20 61 72 67 75 6d  ow "second argum
0840: 65 6e 74 20 74 6f 20 76 65 63 74 6f 72 5f 63 6d  ent to vector_cm
0850: 70 20 6d 69 73 73 69 6e 67 22 3b 0d 0a 0d 0a 20  p missing";.... 
0860: 20 20 20 69 66 28 61 2e 6c 65 6e 67 74 68 20 3e     if(a.length >
0870: 20 62 2e 6c 65 6e 67 74 68 29 7b 0d 0a 20 20 20   b.length){..   
0880: 20 20 20 20 20 20 76 61 72 20 74 6d 70 20 3d 20        var tmp = 
0890: 61 3b 0d 0a 20 20 20 20 20 20 20 20 20 61 20 3d  a;..         a =
08a0: 20 62 3b 20 62 20 3d 20 74 6d 70 3b 20 7d 0d 0a   b; b = tmp; }..
08b0: 0d 0a 20 20 20 20 66 6f 72 28 76 61 72 20 69 20  ..    for(var i 
08c0: 3d 20 30 3b 20 69 20 3c 20 61 2e 6c 65 6e 67 74  = 0; i < a.lengt
08d0: 68 3b 20 2b 2b 69 29 7b 0d 0a 20 20 20 20 20 20  h; ++i){..      
08e0: 20 20 69 66 28 61 5b 69 5d 20 21 3d 20 62 5b 69    if(a[i] != b[i
08f0: 5d 29 20 0d 0a 09 09 09 72 65 74 75 72 6e 20 66  ]) .....return f
0900: 61 6c 73 65 3b 20 7d 0d 0a 0d 0a 20 20 20 20 66  alse; }....    f
0910: 6f 72 28 76 61 72 20 69 20 3d 20 61 2e 6c 65 6e  or(var i = a.len
0920: 67 74 68 3b 20 69 20 3c 20 62 2e 6c 65 6e 67 74  gth; i < b.lengt
0930: 68 3b 20 2b 2b 69 29 7b 0d 0a 20 20 20 20 20 20  h; ++i){..      
0940: 20 20 69 66 28 62 5b 69 5d 20 21 3d 20 30 29 20    if(b[i] != 0) 
0950: 0d 0a 09 09 09 72 65 74 75 72 6e 20 66 61 6c 73  .....return fals
0960: 65 3b 20 7d 0d 0a 20 20 20 20 20 20 20 20 20 20  e; }..          
0970: 20 20 0d 0a 20 20 20 20 72 65 74 75 72 6e 20 74    ..    return t
0980: 72 75 65 3b 20 7d 0d 0a 20 20 20 20 0d 0a 20 20  rue; }..    ..  
0990: 20 20 0d 0a 20 20 20 20 0d 0a 66 75 6e 63 74 69    ..    ..functi
09a0: 6f 6e 20 76 65 63 74 6f 72 5f 6d 69 64 70 6f 69  on vector_midpoi
09b0: 6e 74 28 61 2c 20 62 2c 20 72 65 73 75 6c 74 29  nt(a, b, result)
09c0: 7b 0d 0a 09 69 66 28 21 61 29 20 74 68 72 6f 77  {...if(!a) throw
09d0: 20 22 66 69 72 73 74 20 61 72 67 75 6d 65 6e 74   "first argument
09e0: 20 74 6f 20 76 65 63 74 6f 72 5f 6d 69 64 70 6f   to vector_midpo
09f0: 69 6e 74 20 6d 69 73 73 69 6e 67 22 3b 0d 0a 09  int missing";...
0a00: 69 66 28 21 62 29 20 74 68 72 6f 77 20 22 73 65  if(!b) throw "se
0a10: 63 6f 6e 64 20 61 72 67 75 6d 65 6e 74 20 74 6f  cond argument to
0a20: 20 76 65 63 74 6f 72 5f 6d 69 64 70 6f 69 6e 74   vector_midpoint
0a30: 20 6d 69 73 73 69 6e 67 22 3b 0d 0a 0d 0a 20 20   missing";....  
0a40: 20 20 69 66 28 62 2e 6c 65 6e 67 74 68 20 3e 20    if(b.length > 
0a50: 61 2e 6c 65 6e 67 74 68 29 7b 0d 0a 20 20 20 20  a.length){..    
0a60: 20 20 20 20 76 61 72 20 74 6d 70 20 3d 20 61 3b      var tmp = a;
0a70: 0d 0a 20 20 20 20 20 20 20 20 61 20 3d 20 62 3b  ..        a = b;
0a80: 20 62 20 3d 20 74 6d 70 3b 20 7d 0d 0a 20 20 20   b = tmp; }..   
0a90: 20 20 20 20 20 0d 0a 20 20 20 20 69 66 28 21 72       ..    if(!r
0aa0: 65 73 75 6c 74 29 20 72 65 73 75 6c 74 20 3d 20  esult) result = 
0ab0: 6e 65 77 20 41 72 72 61 79 28 62 2e 6c 65 6e 67  new Array(b.leng
0ac0: 74 68 29 3b 0d 0a 20 20 20 20 20 20 20 20 0d 0a  th);..        ..
0ad0: 20 20 20 20 66 6f 72 28 76 61 72 20 69 20 3d 20      for(var i = 
0ae0: 30 3b 20 69 20 3c 20 61 2e 6c 65 6e 67 74 68 3b  0; i < a.length;
0af0: 20 2b 2b 69 29 7b 0d 0a 20 20 20 20 20 20 20 20   ++i){..        
0b00: 72 65 73 75 6c 74 5b 69 5d 20 3d 20 28 61 5b 69  result[i] = (a[i
0b10: 5d 20 2b 20 62 5b 69 5d 29 2f 32 3b 20 7d 0d 0a  ] + b[i])/2; }..
0b20: 20 20 20 20 20 20 20 20 0d 0a 20 20 20 20 66 6f          ..    fo
0b30: 72 28 76 61 72 20 69 20 3d 20 61 2e 6c 65 6e 67  r(var i = a.leng
0b40: 74 68 3b 20 69 20 3c 20 62 2e 6c 65 6e 67 74 68  th; i < b.length
0b50: 3b 20 2b 2b 69 29 7b 0d 0a 20 20 20 20 20 20 20  ; ++i){..       
0b60: 20 72 65 73 75 6c 74 5b 69 5d 20 3d 20 62 5b 69   result[i] = b[i
0b70: 5d 3b 20 7d 0d 0a 20 20 20 20 0d 0a 20 20 20 20  ]; }..    ..    
0b80: 72 65 74 75 72 6e 20 72 65 73 75 6c 74 3b 0d 0a  return result;..
0b90: 7d 0d 0a 0d 0a 0d 0a 0d 0a 0d 0a 2f 2f 74 68 69  }..........//thi
0ba0: 73 20 66 75 6e 63 74 69 6f 6e 20 6d 61 79 20 6e  s function may n
0bb0: 6f 74 20 63 61 72 65 2c 20 62 75 74 20 74 68 65  ot care, but the
0bc0: 20 65 78 70 65 63 74 65 64 20 6d 61 74 72 69 78   expected matrix
0bd0: 20 66 6f 72 6d 61 74 20 69 6e 20 74 68 69 73 20   format in this 
0be0: 70 72 6f 6a 65 63 74 20 69 73 20 61 20 6c 69 73  project is a lis
0bf0: 74 20 6f 66 20 63 6f 6c 75 6d 6e 20 76 65 63 74  t of column vect
0c00: 6f 72 73 2e 0d 0a 66 75 6e 63 74 69 6f 6e 20 6d  ors...function m
0c10: 61 74 72 69 78 5f 74 72 61 6e 73 70 6f 73 65 28  atrix_transpose(
0c20: 6d 61 74 72 69 78 29 7b 0d 0a 09 69 66 28 21 6d  matrix){...if(!m
0c30: 61 74 72 69 78 29 20 74 68 72 6f 77 20 22 61 72  atrix) throw "ar
0c40: 67 75 6d 65 6e 74 20 74 6f 20 6d 61 74 72 69 78  gument to matrix
0c50: 5f 74 72 61 6e 73 70 6f 73 65 20 6d 69 73 73 69  _transpose missi
0c60: 6e 67 22 3b 0d 0a 0d 0a 20 20 20 20 76 61 72 20  ng";....    var 
0c70: 72 65 73 75 6c 74 20 3d 20 5b 5d 3b 0d 0a 20 20  result = [];..  
0c80: 20 20 76 61 72 20 64 69 6d 30 20 3d 20 6d 61 74    var dim0 = mat
0c90: 72 69 78 2e 6c 65 6e 67 74 68 3b 0d 0a 20 20 20  rix.length;..   
0ca0: 20 69 66 28 21 64 69 6d 30 29 20 0d 0a 09 09 72   if(!dim0) ....r
0cb0: 65 74 75 72 6e 20 72 65 73 75 6c 74 3b 0d 0a 0d  eturn result;...
0cc0: 0a 20 20 20 20 76 61 72 20 64 69 6d 31 20 3d 20  .    var dim1 = 
0cd0: 6d 61 74 72 69 78 5b 30 5d 2e 6c 65 6e 67 74 68  matrix[0].length
0ce0: 3b 0d 0a 0d 0a 20 20 20 20 66 6f 72 28 76 61 72  ;....    for(var
0cf0: 20 69 20 3d 20 30 3b 20 69 20 3c 20 64 69 6d 31   i = 0; i < dim1
0d00: 3b 20 2b 2b 69 29 7b 0d 0a 20 20 20 20 20 20 20  ; ++i){..       
0d10: 20 76 61 72 20 63 6f 6c 76 65 63 74 6f 72 20 3d   var colvector =
0d20: 20 5b 5d 3b 0d 0a 20 20 20 20 20 20 20 20 0d 0a   [];..        ..
0d30: 20 20 20 20 20 20 20 20 66 6f 72 28 76 61 72 20          for(var 
0d40: 6a 20 3d 20 30 3b 20 6a 20 3c 20 64 69 6d 30 3b  j = 0; j < dim0;
0d50: 20 2b 2b 6a 29 7b 0d 0a 20 20 20 20 20 20 20 20   ++j){..        
0d60: 20 20 20 20 63 6f 6c 76 65 63 74 6f 72 2e 70 75      colvector.pu
0d70: 73 68 28 6d 61 74 72 69 78 5b 6a 5d 5b 69 5d 29  sh(matrix[j][i])
0d80: 3b 20 7d 0d 0a 0d 0a 20 20 20 20 20 20 20 20 72  ; }....        r
0d90: 65 73 75 6c 74 2e 70 75 73 68 28 63 6f 6c 76 65  esult.push(colve
0da0: 63 74 6f 72 29 3b 20 7d 0d 0a 20 20 20 20 20 20  ctor); }..      
0db0: 20 20 0d 0a 20 20 20 20 72 65 74 75 72 6e 20 72    ..    return r
0dc0: 65 73 75 6c 74 3b 20 7d 0d 0a 20 20 20 20 0d 0a  esult; }..    ..
0dd0: 0d 0a 0d 0a 20 20 20 20 0d 0a 20 20 20 20 0d 0a  ....    ..    ..
0de0: 0d 0a 66 75 6e 63 74 69 6f 6e 20 6d 61 74 72 69  ..function matri
0df0: 78 5f 6d 75 6c 74 28 61 2c 20 62 29 7b 0d 0a 09  x_mult(a, b){...
0e00: 69 66 28 21 61 29 20 74 68 72 6f 77 20 22 66 69  if(!a) throw "fi
0e10: 72 73 74 20 61 72 67 75 6d 65 6e 74 20 74 6f 20  rst argument to 
0e20: 6d 61 74 72 69 78 5f 6d 75 6c 74 20 6d 69 73 73  matrix_mult miss
0e30: 69 6e 67 22 3b 0d 0a 09 69 66 28 21 62 29 20 74  ing";...if(!b) t
0e40: 68 72 6f 77 20 22 73 65 63 6f 6e 64 20 61 72 67  hrow "second arg
0e50: 75 6d 65 6e 74 20 74 6f 20 6d 61 74 72 69 78 5f  ument to matrix_
0e60: 6d 75 6c 74 20 6d 69 73 73 69 6e 67 22 3b 0d 0a  mult missing";..
0e70: 0d 0a 20 20 20 20 76 61 72 20 72 65 73 75 6c 74  ..    var result
0e80: 20 3d 20 5b 5d 3b 0d 0a 20 20 20 20 76 61 72 20   = [];..    var 
0e90: 61 20 3d 20 6d 61 74 72 69 78 5f 74 72 61 6e 73  a = matrix_trans
0ea0: 70 6f 73 65 28 61 29 3b 0d 0a 20 20 20 20 76 61  pose(a);..    va
0eb0: 72 20 64 69 6d 30 20 3d 20 61 2e 6c 65 6e 67 74  r dim0 = a.lengt
0ec0: 68 3b 0d 0a 20 20 20 20 76 61 72 20 64 69 6d 31  h;..    var dim1
0ed0: 20 3d 20 62 2e 6c 65 6e 67 74 68 3b 0d 0a 0d 0a   = b.length;....
0ee0: 20 20 20 20 66 6f 72 28 76 61 72 20 69 20 3d 20      for(var i = 
0ef0: 30 3b 20 69 20 3c 20 64 69 6d 31 3b 20 2b 2b 69  0; i < dim1; ++i
0f00: 29 7b 0d 0a 20 20 20 20 20 20 20 20 76 61 72 20  ){..        var 
0f10: 63 6f 6c 76 65 63 74 6f 72 20 3d 20 5b 5d 3b 0d  colvector = [];.
0f20: 0a 20 20 20 20 20 20 20 20 0d 0a 20 20 20 20 20  .        ..     
0f30: 20 20 20 66 6f 72 28 76 61 72 20 6a 20 3d 20 30     for(var j = 0
0f40: 3b 20 6a 20 3c 20 64 69 6d 30 3b 20 2b 2b 6a 29  ; j < dim0; ++j)
0f50: 7b 0d 0a 20 20 20 20 20 20 20 20 20 20 20 20 63  {..            c
0f60: 6f 6c 76 65 63 74 6f 72 2e 70 75 73 68 28 20 76  olvector.push( v
0f70: 65 63 74 6f 72 5f 64 6f 74 28 61 5b 6a 5d 2c 20  ector_dot(a[j], 
0f80: 62 5b 69 5d 29 20 29 3b 20 20 7d 0d 0a 20 20 20  b[i]) );  }..   
0f90: 20 20 20 20 20 72 65 73 75 6c 74 2e 70 75 73 68       result.push
0fa0: 28 63 6f 6c 76 65 63 74 6f 72 29 3b 20 7d 0d 0a  (colvector); }..
0fb0: 0d 0a 20 20 20 20 72 65 74 75 72 6e 20 72 65 73  ..    return res
0fc0: 75 6c 74 3b 20 7d 0d 0a 0d 0a 0d 0a 20 20 20 20  ult; }......    
0fd0: 0d 0a 2f 2f 20 66 69 6e 64 73 20 74 68 65 20 69  ..// finds the i
0fe0: 6e 76 65 72 73 65 20 6f 66 20 61 20 33 20 62 79  nverse of a 3 by
0ff0: 20 33 20 6d 61 74 72 69 78 2c 20 72 65 74 75 72   3 matrix, retur
1000: 6e 73 20 66 61 6c 73 65 20 69 66 20 6d 61 74 72  ns false if matr
1010: 69 78 20 68 61 73 20 64 65 74 65 72 6d 69 6e 61  ix has determina
1020: 6e 74 20 7a 65 72 6f 2e 0d 0a 2f 2f 0d 0a 0d 0a  nt zero...//....
1030: 66 75 6e 63 74 69 6f 6e 20 6d 61 74 72 69 78 33  function matrix3
1040: 33 69 6e 76 28 6d 29 7b 0d 0a 20 20 20 20 69 66  3inv(m){..    if
1050: 28 21 6d 29 20 74 68 72 6f 77 20 22 61 72 67 75  (!m) throw "argu
1060: 6d 65 6e 74 20 74 6f 20 6d 61 74 72 69 78 33 33  ment to matrix33
1070: 69 6e 76 20 6d 69 73 73 69 6e 67 22 3b 0d 0a 0d  inv missing";...
1080: 0a 20 20 20 20 2f 2f 20 74 68 65 20 6d 61 74 72  .    // the matr
1090: 69 78 20 66 6f 72 6d 61 74 20 69 73 20 61 20 6c  ix format is a l
10a0: 69 73 74 20 6f 66 20 63 6f 6c 75 6d 6e 20 76 65  ist of column ve
10b0: 63 74 6f 72 73 2e 0d 0a 20 20 20 20 76 61 72 20  ctors...    var 
10c0: 61 20 3d 20 6d 5b 30 5d 5b 30 5d 3b 0d 0a 09 76  a = m[0][0];...v
10d0: 61 72 20 62 20 3d 20 6d 5b 31 5d 5b 30 5d 3b 0d  ar b = m[1][0];.
10e0: 0a 09 76 61 72 20 63 20 3d 20 6d 5b 32 5d 5b 30  ..var c = m[2][0
10f0: 5d 3b 0d 0a 09 76 61 72 20 64 20 3d 20 6d 5b 30  ];...var d = m[0
1100: 5d 5b 31 5d 3b 0d 0a 09 76 61 72 20 65 20 3d 20  ][1];...var e = 
1110: 6d 5b 31 5d 5b 31 5d 3b 0d 0a 09 76 61 72 20 66  m[1][1];...var f
1120: 20 3d 20 6d 5b 32 5d 5b 31 5d 3b 0d 0a 09 76 61   = m[2][1];...va
1130: 72 20 67 20 3d 20 6d 5b 30 5d 5b 32 5d 3b 0d 0a  r g = m[0][2];..
1140: 09 76 61 72 20 68 20 3d 20 6d 5b 31 5d 5b 32 5d  .var h = m[1][2]
1150: 3b 0d 0a 09 76 61 72 20 6b 20 3d 20 6d 5b 32 5d  ;...var k = m[2]
1160: 5b 32 5d 3b 0d 0a 09 0d 0a 09 2f 2f 20 5b 20 61  [2];......// [ a
1170: 20 62 20 63 0d 0a 20 20 20 20 2f 2f 20 20 20 64   b c..    //   d
1180: 20 65 20 66 0d 0a 20 20 20 20 2f 2f 20 20 20 67   e f..    //   g
1190: 20 68 20 6b 20 5d 0d 0a 09 0d 0a 09 0d 0a 20 20   h k ]........  
11a0: 20 20 76 61 72 20 41 20 3d 20 65 2a 6b 20 2d 20    var A = e*k - 
11b0: 66 2a 68 3b 0d 0a 09 76 61 72 20 42 20 3d 20 66  f*h;...var B = f
11c0: 2a 67 20 2d 20 64 2a 6b 3b 0d 0a 09 76 61 72 20  *g - d*k;...var 
11d0: 43 20 3d 20 64 2a 68 20 2d 20 65 2a 67 3b 09 0d  C = d*h - e*g;..
11e0: 0a 09 76 61 72 20 44 20 3d 20 68 2a 63 20 2d 20  ..var D = h*c - 
11f0: 62 2a 6b 3b 0d 0a 09 76 61 72 20 45 20 3d 20 61  b*k;...var E = a
1200: 2a 6b 20 2d 20 63 2a 67 3b 0d 0a 09 76 61 72 20  *k - c*g;...var 
1210: 46 20 3d 20 67 2a 62 20 2d 20 61 2a 68 3b 0d 0a  F = g*b - a*h;..
1220: 09 76 61 72 20 47 20 3d 20 62 2a 66 20 2d 20 65  .var G = b*f - e
1230: 2a 63 3b 0d 0a 09 76 61 72 20 48 20 3d 20 64 2a  *c;...var H = d*
1240: 63 20 2d 20 61 2a 66 3b 0d 0a 09 76 61 72 20 4b  c - a*f;...var K
1250: 20 3d 20 61 2a 65 20 2d 20 64 2a 62 3b 0d 0a 09   = a*e - d*b;...
1260: 0d 0a 09 0d 0a 20 20 20 20 76 61 72 20 64 65 74  .....    var det
1270: 20 3d 20 28 61 2a 41 20 2b 20 62 2a 42 20 2b 20   = (a*A + b*B + 
1280: 63 2a 43 29 3b 0d 0a 09 0d 0a 09 0d 0a 09 69 66  c*C);.........if
1290: 28 64 65 74 20 3d 3d 20 30 29 20 72 65 74 75 72  (det == 0) retur
12a0: 6e 20 66 61 6c 73 65 3b 0d 0a 09 0d 0a 09 2f 2f  n false;......//
12b0: 20 74 68 65 20 6d 75 6c 74 69 70 6c 69 63 61 74   the multiplicat
12c0: 69 76 65 20 69 6e 76 65 72 73 65 20 6f 66 20 74  ive inverse of t
12d0: 68 65 20 64 65 74 65 72 6d 69 6e 61 6e 74 2e 0d  he determinant..
12e0: 0a 09 76 61 72 20 64 69 20 3d 20 31 2f 64 65 74  ..var di = 1/det
12f0: 3b 0d 0a 09 0d 0a 09 2f 2f 20 74 68 65 20 6d 61  ;......// the ma
1300: 74 72 69 78 20 66 6f 72 6d 61 74 20 69 73 20 61  trix format is a
1310: 20 6c 69 73 74 20 6f 66 20 63 6f 6c 75 6d 6e 20   list of column 
1320: 76 65 63 74 6f 72 73 2e 0d 0a 09 76 61 72 20 72  vectors....var r
1330: 65 73 75 6c 74 20 3d 20 5b 5b 64 69 2a 41 2c 20  esult = [[di*A, 
1340: 64 69 2a 42 2c 20 64 69 2a 43 5d 2c 0d 0a 09 09  di*B, di*C],....
1350: 20 20 20 20 20 20 20 20 20 20 5b 64 69 2a 44 2c            [di*D,
1360: 20 64 69 2a 45 2c 20 64 69 2a 46 5d 2c 0d 0a 09   di*E, di*F],...
1370: 09 20 20 20 20 20 20 20 20 20 20 5b 64 69 2a 47  .          [di*G
1380: 2c 20 64 69 2a 48 2c 20 64 69 2a 4b 5d 5d 3b 0d  , di*H, di*K]];.
1390: 0a 0d 0a 09 72 65 74 75 72 6e 20 72 65 73 75 6c  ....return resul
13a0: 74 3b 20 20 7d 0d 0a 0d 0a 09 0d 0a 2f 2f 20 64  t;  }.......// d
13b0: 6f 65 73 20 6e 6f 74 20 63 68 65 63 6b 20 69 66  oes not check if
13c0: 20 66 75 6e 63 74 69 6f 6e 20 69 73 20 64 69 66   function is dif
13d0: 66 65 72 65 6e 74 69 61 62 6c 65 2e 0d 0a 0d 0a  ferentiable.....
13e0: 66 75 6e 63 74 69 6f 6e 20 63 6f 6d 70 75 74 65  function compute
13f0: 5f 67 72 61 64 69 65 6e 74 28 66 75 6e 63 2c 20  _gradient(func, 
1400: 70 6f 69 6e 74 2c 20 65 72 72 6f 72 29 7b 0d 0a  point, error){..
1410: 09 69 66 28 21 66 75 6e 63 29 20 74 68 72 6f 77  .if(!func) throw
1420: 20 22 66 75 6e 63 74 69 6f 6e 20 61 72 67 75 6d   "function argum
1430: 65 6e 74 20 74 6f 20 63 6f 6d 70 75 74 65 5f 67  ent to compute_g
1440: 72 61 64 69 65 6e 74 20 6d 69 73 73 69 6e 67 22  radient missing"
1450: 3b 0d 0a 09 69 66 28 21 70 6f 69 6e 74 29 20 74  ;...if(!point) t
1460: 68 72 6f 77 20 22 70 6f 69 6e 74 20 61 72 67 75  hrow "point argu
1470: 6d 65 6e 74 20 74 6f 20 63 6f 6d 70 75 74 65 5f  ment to compute_
1480: 67 72 61 64 69 65 6e 74 20 6d 69 73 73 69 6e 67  gradient missing
1490: 22 3b 0d 0a 0d 0a 20 20 20 20 69 66 28 21 65 72  ";....    if(!er
14a0: 72 6f 72 29 20 65 72 72 6f 72 20 3d 20 30 2e 30  ror) error = 0.0
14b0: 30 30 30 31 3b 0d 0a 0d 0a 20 20 20 20 2f 2f 20  0001;....    // 
14c0: 66 69 6e 64 20 67 72 61 64 69 65 6e 74 0d 0a 20  find gradient.. 
14d0: 20 20 20 76 61 72 20 6f 74 68 65 72 70 6f 69 6e     var otherpoin
14e0: 74 20 3d 20 70 6f 69 6e 74 2e 73 6c 69 63 65 28  t = point.slice(
14f0: 30 29 3b 0d 0a 20 20 20 20 76 61 72 20 67 72 61  0);..    var gra
1500: 64 69 65 6e 74 20 3d 20 6e 65 77 20 41 72 72 61  dient = new Arra
1510: 79 28 70 6f 69 6e 74 2e 6c 65 6e 67 74 68 29 3b  y(point.length);
1520: 0d 0a 20 20 20 20 76 61 72 20 70 6f 69 6e 74 76  ..    var pointv
1530: 61 6c 75 65 20 3d 20 66 75 6e 63 28 70 6f 69 6e  alue = func(poin
1540: 74 29 3b 0d 0a 20 20 20 20 0d 0a 20 20 20 20 66  t);..    ..    f
1550: 6f 72 28 76 61 72 20 69 20 3d 20 30 3b 20 69 20  or(var i = 0; i 
1560: 3c 20 70 6f 69 6e 74 2e 6c 65 6e 67 74 68 3b 20  < point.length; 
1570: 2b 2b 69 29 7b 0d 0a 20 20 20 20 20 20 20 20 6f  ++i){..        o
1580: 74 68 65 72 70 6f 69 6e 74 5b 69 5d 20 2b 3d 20  therpoint[i] += 
1590: 65 72 72 6f 72 3b 0d 0a 20 20 20 20 20 20 20 20  error;..        
15a0: 76 61 72 20 6f 74 68 65 72 76 61 6c 75 65 20 3d  var othervalue =
15b0: 20 66 75 6e 63 28 6f 74 68 65 72 70 6f 69 6e 74   func(otherpoint
15c0: 29 0d 0a 20 20 20 20 20 20 20 20 67 72 61 64 69  )..        gradi
15d0: 65 6e 74 5b 69 5d 20 3d 20 28 6f 74 68 65 72 76  ent[i] = (otherv
15e0: 61 6c 75 65 20 2d 20 70 6f 69 6e 74 76 61 6c 75  alue - pointvalu
15f0: 65 29 20 2f 20 65 72 72 6f 72 3b 0d 0a 20 20 20  e) / error;..   
1600: 20 20 20 20 20 0d 0a 20 20 20 20 20 20 20 20 6f       ..        o
1610: 74 68 65 72 70 6f 69 6e 74 5b 69 5d 20 2d 3d 20  therpoint[i] -= 
1620: 65 72 72 6f 72 3b 20 7d 0d 0a 20 20 20 20 0d 0a  error; }..    ..
1630: 20 20 20 20 72 65 74 75 72 6e 20 67 72 61 64 69      return gradi
1640: 65 6e 74 3b 0d 0a 7d 0d 0a 20 20 20 20 0d 0a 20  ent;..}..    .. 
1650: 20 20 20 0d 0a 20 20 20 20 0d 0a 0d 0a 2f 2f 0d     ..    ....//.
1660: 0a 2f 2f 20 65 71 75 61 74 69 6f 6e 20 61 74 3a  .// equation at:
1670: 0d 0a 2f 2f 20 68 74 74 70 3a 2f 2f 69 6e 73 69  ..// http://insi
1680: 64 65 2e 6d 69 6e 65 73 2e 65 64 75 2f 7e 67 6d  de.mines.edu/~gm
1690: 75 72 72 61 79 2f 41 72 62 69 74 72 61 72 79 41  urray/ArbitraryA
16a0: 78 69 73 52 6f 74 61 74 69 6f 6e 2f 41 72 62 69  xisRotation/Arbi
16b0: 74 72 61 72 79 41 78 69 73 52 6f 74 61 74 69 6f  traryAxisRotatio
16c0: 6e 2e 68 74 6d 6c 0d 0a 2f 2f 20 63 72 65 61 74  n.html..// creat
16d0: 65 73 20 61 20 72 6f 74 61 74 69 6f 6e 20 6d 61  es a rotation ma
16e0: 74 72 69 78 20 61 62 6f 75 74 20 74 68 65 20 67  trix about the g
16f0: 69 76 65 6e 20 76 65 63 74 6f 72 2e 0d 0a 2f 2f  iven vector...//
1700: 0d 0a 66 75 6e 63 74 69 6f 6e 20 76 65 63 74 6f  ..function vecto
1710: 72 5f 72 6f 74 61 74 69 6f 6e 28 76 65 63 74 6f  r_rotation(vecto
1720: 72 2c 20 73 29 7b 0d 0a 09 69 66 28 21 76 65 63  r, s){...if(!vec
1730: 74 6f 72 29 20 74 68 72 6f 77 20 22 76 65 63 74  tor) throw "vect
1740: 6f 72 20 61 72 67 75 6d 65 6e 74 20 74 6f 20 76  or argument to v
1750: 65 63 74 6f 72 5f 72 6f 74 61 74 69 6f 6e 20 6d  ector_rotation m
1760: 69 73 73 69 6e 67 22 3b 0d 0a 09 69 66 28 73 20  issing";...if(s 
1770: 3d 3d 3d 20 6e 75 6c 6c 20 7c 7c 20 73 20 3d 3d  === null || s ==
1780: 3d 20 75 6e 64 65 66 69 6e 65 64 29 20 74 68 72  = undefined) thr
1790: 6f 77 20 22 61 6e 67 6c 65 20 61 72 67 75 6d 65  ow "angle argume
17a0: 6e 74 20 74 6f 20 76 65 63 74 6f 72 5f 72 6f 74  nt to vector_rot
17b0: 61 74 69 6f 6e 20 6d 69 73 73 69 6e 67 22 3b 0d  ation missing";.
17c0: 0a 09 0d 0a 20 20 20 20 76 61 72 20 6c 65 6e 20  ....    var len 
17d0: 3d 20 4d 61 74 68 2e 73 71 72 74 28 76 65 63 74  = Math.sqrt(vect
17e0: 6f 72 5b 30 5d 2a 76 65 63 74 6f 72 5b 30 5d 0d  or[0]*vector[0].
17f0: 0a 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20  .               
1800: 20 20 20 20 20 20 20 20 20 2b 20 76 65 63 74 6f           + vecto
1810: 72 5b 31 5d 2a 76 65 63 74 6f 72 5b 31 5d 0d 0a  r[1]*vector[1]..
1820: 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20                  
1830: 20 20 20 20 20 20 20 20 2b 20 76 65 63 74 6f 72          + vector
1840: 5b 32 5d 2a 76 65 63 74 6f 72 5b 32 5d 29 3b 0d  [2]*vector[2]);.
1850: 0a 0d 0a 20 20 20 20 76 61 72 20 75 20 3d 20 76  ...    var u = v
1860: 65 63 74 6f 72 5b 30 5d 2f 6c 65 6e 3b 0d 0a 20  ector[0]/len;.. 
1870: 20 20 20 76 61 72 20 76 20 3d 20 76 65 63 74 6f     var v = vecto
1880: 72 5b 31 5d 2f 6c 65 6e 3b 0d 0a 20 20 20 20 76  r[1]/len;..    v
1890: 61 72 20 77 20 3d 20 76 65 63 74 6f 72 5b 32 5d  ar w = vector[2]
18a0: 2f 6c 65 6e 3b 0d 0a 20 20 20 20 0d 0a 20 20 20  /len;..    ..   
18b0: 20 76 61 72 20 5f 63 6f 73 20 3d 20 4d 61 74 68   var _cos = Math
18c0: 2e 63 6f 73 28 73 29 3b 0d 0a 20 20 20 20 76 61  .cos(s);..    va
18d0: 72 20 5f 73 69 6e 20 3d 20 4d 61 74 68 2e 73 69  r _sin = Math.si
18e0: 6e 28 73 29 3b 0d 0a 20 20 20 20 76 61 72 20 5f  n(s);..    var _
18f0: 31 6d 63 6f 73 20 3d 20 31 20 2d 20 5f 63 6f 73  1mcos = 1 - _cos
1900: 3b 0d 0a 0d 0a 20 20 20 20 72 65 74 75 72 6e 20  ;....    return 
1910: 5b 5b 75 2a 75 2a 5f 31 6d 63 6f 73 20 2b 20 5f  [[u*u*_1mcos + _
1920: 63 6f 73 2c 20 75 2a 76 2a 5f 31 6d 63 6f 73 20  cos, u*v*_1mcos 
1930: 2d 20 77 2a 5f 73 69 6e 2c 20 75 2a 77 2a 5f 31  - w*_sin, u*w*_1
1940: 6d 63 6f 73 20 2b 20 76 2a 5f 73 69 6e 5d 2c 0d  mcos + v*_sin],.
1950: 0a 20 20 20 20 20 20 20 20 20 20 20 20 5b 75 2a  .            [u*
1960: 76 2a 5f 31 6d 63 6f 73 20 2b 20 77 2a 5f 73 69  v*_1mcos + w*_si
1970: 6e 2c 20 76 2a 76 2a 5f 31 6d 63 6f 73 20 2b 20  n, v*v*_1mcos + 
1980: 5f 63 6f 73 2c 20 76 2a 77 2a 5f 31 6d 63 6f 73  _cos, v*w*_1mcos
1990: 20 2d 20 75 2a 5f 73 69 6e 5d 2c 0d 0a 20 20 20   - u*_sin],..   
19a0: 20 20 20 20 20 20 20 20 20 5b 75 2a 77 2a 5f 31           [u*w*_1
19b0: 6d 63 6f 73 20 2d 20 76 2a 5f 73 69 6e 2c 20 76  mcos - v*_sin, v
19c0: 2a 77 2a 5f 31 6d 63 6f 73 20 2b 20 75 2a 5f 73  *w*_1mcos + u*_s
19d0: 69 6e 2c 20 77 2a 77 2a 5f 31 6d 63 6f 73 20 2b  in, w*w*_1mcos +
19e0: 20 5f 63 6f 73 5d 5d 3b 20 7d 0d 0a               _cos]]; }..