real function f1(xarg)
  implicit none
  real ::xarg
  
  real ::p1,p2,p3
  
  p1=xarg
  p2=xarg**2
  p3=xarg**3
  f1=p1+p2+p3
end function f1