function y=cbrt(x);
y=x.^(1/3);