function y = pow(x,n); y = x.^n;