C
        SUBROUTINE PEScvqz(V1,x,y,z)
        IMPLICIT DOUBLE PRECISION (A-H,O-Z)
        Degree=1.0d0
        E=2.718281828d0
C

        v1=367.4546631599754-427.2312909598234*(x-y)**2- 
     -  89.2867191912118*(x - y)**4 -
     -  0.6160025866071237*(x - y)**6 + 
     -  992.822414424169*(-1.3378d0 + 0.5d0*(x + y)) + 
     -  1366.032503523549*(x-y)**2*(-1.3378d0+0.5d0*(x+y))+ 
     -  210.5770888210323*(x-y)**4*(-1.3378d0+0.5d0*(x+y))- 
     -  1797.202403342277*(-1.3378d0 + 0.5d0*(x + y))**2 - 
     -  1883.307050371728*(x-y)**2*(-1.3378d0+0.5d0*(x+y))**2- 
     -  155.7503191002783*(x-y)**4*(-1.3378d0+0.5d0*(x+y))**2+ 
     -  1776.033583531499*(-1.3378d0 + 0.5d0*(x + y))**3 + 
     -  1138.676338966165*(x-y)**2*(-1.3378d0+0.5d0*(x+y))**3 - 
     -  1015.792295031244*(-1.3378d0+0.5d0*(x+y))**4+ 
     -  629.2208896781913*(-1.3378d0 + 0.5d0*(x + y))**5 - 
     -  775.9881693988549*(-1.3378d0 + 0.5d0*(x + y))**6 - 
     -  35.30930640108655*(0.04032362334446276+
     -  Cos(Degree*z)) + 
     -  16.38765482902978*(x - y)**2*(0.04032362334446276+
     -  Cos(Degree*z)) + 
     -  41.66794689019272*(-1.3378d0 + 0.5d0*(x + y))*
     -   (0.04032362334446276 + Cos(Degree*z)) + 
     -  15.18267969667999*(x-y)**2*(-1.3378d0+0.5d0*(x+y))*
     -   (0.04032362334446276 + Cos(Degree*z)) - 
     -  41.93598321052467*(x-y)**4*(-1.3378d0+0.5d0*(x+y))*
     -   (0.04032362334446276 + Cos(Degree*z)) + 
     -  17.67928209891903*(-1.3378d0 + 0.5d0*(x + y))**2*
     -   (0.04032362334446276 + Cos(Degree*z)) - 
     -  25.1933118396414*(x-y)**2*(-1.3378d0+0.5d0*(x+y))**2*
     -   (0.04032362334446276 + Cos(Degree*z)) + 
     -  3.376367187698178*(-1.3378d0 + 0.5d0*(x + y))**3*
     -   (0.04032362334446276 + Cos(Degree*z)) - 
     -  39.1551625418906*(-1.3378d0 + 0.5d0*(x + y))**4*
     -   (0.04032362334446276 + Cos(Degree*z)) - 
     -  39.52493130897267*(0.04032362334446276 + 
     -  Cos(Degree*z))**2 + 
     -  19.77557605502124*(x - y)**2*
     -   (0.04032362334446276 + Cos(Degree*z))**2 + 
     -  4.97750975514657*(x - y)**4*
     -   (0.04032362334446276 + Cos(Degree*z))**2 + 
     -  1.69156787407384*(-1.3378d0 + 0.5d0*(x + y))*
     -   (0.04032362334446276 + Cos(Degree*z))**2 + 
     -  20.07409377446032*(x-y)**2*(-1.3378d0+0.5d0*(x+y))*
     -   (0.04032362334446276 + Cos(Degree*z))**2 + 
     -  104.6464240770881*(-1.3378d0 + 0.5d0*(x + y))**2*
     -   (0.04032362334446276 + Cos(Degree*z))**2 + 
     -  71.18549298799598*(-1.3378d0 + 0.5d0*(x + y))**3*
     -   (0.04032362334446276 + Cos(Degree*z))**2 + 
     -  88.7982045948556*(-1.3378d0 + 0.5d0*(x + y))**4*
     -   (0.04032362334446276 + Cos(Degree*z))**2 - 
     -  45.94604668270892*(0.04032362334446276 + 
     -  Cos(Degree*z))**3 + 
     -  15.32192426225214*(x - y)**2*
     -   (0.04032362334446276 + Cos(Degree*z))**3 + 
     -  30.14322287839019*(-1.3378d0 + 0.5d0*(x + y))*
     -   (0.04032362334446276 + Cos(Degree*z))**3 + 
     -  148.0234744105225*(-1.3378d0 + 0.5d0*(x + y))**2*
     -   (0.04032362334446276 + Cos(Degree*z))**3 + 
     -  166.3946700986707*(-1.3378d0 + 0.5d0*(x + y))**3*
     -   (0.04032362334446276 + Cos(Degree*z))**3 - 
     -  414.3848562797138*(-1.3378d0 + 0.5d0*(x + y))**4*
     -   (0.04032362334446276 + Cos(Degree*z))**3 - 
     -  24.24860722565262*(0.04032362334446276 + 
     -  Cos(Degree*z))**4 - 
     -  65.16355631168778*(x - y)**2*
     -   (0.04032362334446276 + Cos(Degree*z))**4 + 
     -  38.57998737409728*(-1.3378d0 + 0.5d0*(x + y))*
     -   (0.04032362334446276 + Cos(Degree*z))**4 - 
     -  540.4081770075879*(-1.3378d0 + 0.5d0*(x + y))**2*
     -   (0.04032362334446276 + Cos(Degree*z))**4 - 
     -  422.0788151937211*(-1.3378d0 + 0.5d0*(x + y))**3*
     -   (0.04032362334446276 + Cos(Degree*z))**4 + 
     -  134.2117929813882*(0.04032362334446276 + 
     -  Cos(Degree*z))**5 - 
     -  197.0809475555731*(-1.3378d0 + 0.5d0*(x + y))*
     -   (0.04032362334446276 + Cos(Degree*z))**5 - 
     -  18.17098672436418*(-1.3378d0 + 0.5d0*(x + y))**2*
     -   (0.04032362334446276 + Cos(Degree*z))**5 + 
     -  125.0935290711315*(0.04032362334446276 + 
     -  Cos(Degree*z))**6 + 
     -  130.7504240131748*(x - y)**2*
     -   (0.04032362334446276 + Cos(Degree*z))**6 - 
     -  147.3176696896054*(-1.3378d0 + 0.5d0*(x + y))*
     -   (0.04032362334446276 + Cos(Degree*z))**6 + 
     -  959.238970740647*(-1.3378d0 + 0.5d0*(x + y))**2*
     -   (0.04032362334446276 + Cos(Degree*z))**6 - 
     -  117.8316665428289*(0.04032362334446276 + 
     -  Cos(Degree*z))**7 - 
     -  135.4388288253769*(0.04032362334446276 + 
     -  Cos(Degree*z))**8
        RETURN
         END