function y = testfun (x)

x = x/2;

sx = localsqrt(x);

y = log(sx);

end

function b = localsqrt (a)

  b = sqrt (a);
  
endfunction