function f=powerfun(x,p)
    f=x^p;