      function [rotate] = rot_i
%     ---------------------------------------
%     Isometric view rotation matrix 
%     (45 deg y-rot then 35.26 deg x-rot)
%     ---------------------------------------
%     rotate = rotation matrix (sub-set of finite)
      rotate   = eye(3)   ; % define sizes
      rotate_y = eye(3)   ;
      rotate_x = eye(3)   ;
%       y-rotation
      ang_y      = pi*45./180.   ;
      [rotate_y] = rot_y (ang_y)  ;
%c    disp(rotate_y)  ;
%       x-rotation
      ang_x      = pi*35.26/180. ;
      [rotate_x] = rot_x (ang_x)  ;
%c    disp(rotate_x)  ;
%      net rotation matrix, right to left
      rotate = rotate_x*rotate_y  ;
%     end of rot_i
