10 f = rotvec([1;2;3],pi/2);
 
   12 q_ = quaternion([0 0 1], pi/4);
 
   14 q_ = quaternion([pi/4, pi/6, pi/8, 1]);
 
   16 h  = qu_.r;    % Rotation Matrix from quaternion
 
   17 i  = qu_.t;    % Transformation 
"            " 
   23 % -------------- R O B O T -----------------
 
   26 q_   = [pi/4 pi/4 pi/4 pi/4 pi/4 pi/4];
 
   33 o = rne_dh(p560, q_, qd_, qdd_);
 
   34 p = inertia(p560, q_);
 
   37 q = rne_dh(p560, q_, qd_, qdd_);
 
   39 run stanford_no_motor;
 
   41 s = jacob0(stanf, q_);
 
   42 t = jacobn(stanf, q_);
 
   43 u = rne_dh(stanf, q_, qd_, qdd_);
 
   45 f_ = [10;5;7;11;22;33];
 
   46 v = rne_dh(stanf, q_, qd_, qdd_, g_, f_);
 
   47 w = inertia(stanf, q_);
 
   49 run puma560akb_no_motor;
 
   51 y = jacob0(p560m, q_);
 
   52 z = jacobn(p560m, q_);
 
   53 aa = rne_mdh(p560m, q_, qd_, qdd_);
 
   54 bb= inertia(p560m, q_);
 
   57 cc = rne_mdh(p560m, q_, qd_, qdd_);
 
   58 dd = rne_mdh(p560m, q_, qd_, qdd_, g_, f_);
 
   60 save ../source/test.txt -ascii -
double a b c d e f g h i j k l m n o p q r s t u v w x y z aa bb cc dd