dyn.load("vecsqr.so")    # load the compiled code

vecsqr = function(x)     # our R interface to the compiled code
{
  n = length(x)
  result  = .C("vecsqr",
        y = double(n),
        as.double(x),
        as.integer(n) )
  result[["y"]]
}