function fun(x)
      common/pawpar/pa(1)
      fun=sin(pa(1)*x)
      end