subroutine sayhello(x,r,v)
        double precision x,r,v
        r =        (sin(2.7646015351590181E+3*((v-1.0000000000000001E-1)*v**2*sin(1.6
     1   587609210954106E+4*x)/3.6E+3+x))+v**(5.0E+0/3.0E+0)*sin(1.38230
     2   0767579509E+4*(v**(4.0E+0/3.0E+0)*sin(1.1058406140636072E+4*x)/
     3   2.7182818284590451E+0**(8*x)/4.0E+2+x))/5.0E+0+v**(1.0E+1/3.0E+
     4   0)*sin(2.764601535159018E+4*x)/1.0E+1+v**3*sin(2.48814138164311
     5   63E+4*x)/9.0E+0+v**(8.0E+0/3.0E+0)*sin(2.2116812281272145E+4*x)
     6   /8.0E+0+v**(7.0E+0/3.0E+0)*sin(1.9352210746113127E+4*x)/7.0E+0+
     7   v**2*sin(1.6587609210954106E+4*x)/6.0E+0+v**(4.0E+0/3.0E+0)*sin
     8   (1.1058406140636072E+4*x)/4.0E+0+sin(8.293804605477053E+3*x)/2.
     9   7182818284590451E+0**(8*x)+v**(2.0E+0/3.0E+0)*sin(5.52920307031
     :   80362E+3*x)/2.0E+0)/2.7182818284590451E+0**(3*x)/2.0E+0
        return
      end