function test_OpHat()
  d0 = %pi
  d2 = [1 2 3 4 ; 5 6 7 8 ; 9 10 11 12 ; 13 14 15 16]
  z0 = d0 + %i * d0
  z2 = d2 + %i * d2
  
  disp(d0^d0)
  disp(z0^z0)
  disp(d2^d0)
  disp(z2^z0)

  disp(d0^d2)
  disp(z0^z2)

  //disp(d2^d2) // ! Not implemented in Scilab !
  //disp(z2^z2) // ! Not implemented in Scilab !
  
  disp(d2^z0)
  disp(z2^d0)
  
  disp(d0^z2)
  disp(z0^d2)
  
  disp(d0^z0)
  disp(z0^d0)

  //disp(d2^z2) // ! Not implemented in Scilab !
  //disp(z2^d2) // ! Not implemented in Scilab !
endfunction