function [fx]=rad(x)

fx=sign(x).*sqrt(abs(x));