functiony=cbrtx(x)%CBRTX The real cube root%% CBRTX(x) is the real cube root of x (assuming x is real). x% can be any shape.persistentoctavepifisempty(octavep)octavep=exist('OCTAVE_VERSION','builtin')~=0;endifoctavepy=cbrt(x);elsey=abs(x).^(1/3);y(x<0)=-y(x<0);y(x==0)=x(x==0);endend