function y=paramUpdate(x,j,k) if x(1)==0 error('Initial parameter must be greater than 0'); end y=x(1)*(x(2)/x(1))^(j/k);