mpower

[Eric and Diana Luft]

Is there any Octave equivalent to the matlab "mpower" command, which would compute a matrix to an arbitrary, eg fractional, power?  I want to compute a matrix to the 2/3 power, but I get a complex-valued matrix as a result, where it should be real-valued.  Details follow below, where A2 comes out to be complex-valued when I think it should be real-valued:

 

M=[1  1  1; 1   2  2; 1  2  3]

K=[2 -1 .5;-1 1.2 .4;.5 .4 1.8]

A2=logm(inv(M)*K)^(2/3)

 

Thanks for any help or advice,

 

Eric Luft

 

address@hidden             http://www.lds.org