22 lines
377 B
Mathematica
22 lines
377 B
Mathematica
|
function y = mpower(x,d)
|
||
|
|
%MPOWER (overloaded)
|
||
|
|
|
||
|
|
%Sanity check
|
||
|
|
if prod(size(d))>1 | ~((fix(d) == d)) | ~(d>=0)
|
||
|
|
error('The power must be scalar.');
|
||
|
|
end
|
||
|
|
if x.dim(1)~=x.dim(2)
|
||
|
|
error('Matrix must be square.')
|
||
|
|
end
|
||
|
|
|
||
|
|
switch d
|
||
|
|
case 0
|
||
|
|
y = eye(x.dim(1),x.dim(2))^0;
|
||
|
|
case 1
|
||
|
|
y = x;
|
||
|
|
case 2
|
||
|
|
y = x * x;
|
||
|
|
otherwise
|
||
|
|
y = x*mpower(x,d-1);
|
||
|
|
end
|