C
C
        SUBROUTINE PESrel(V1,x,y,z)
        IMPLICIT DOUBLE PRECISION (A-H,O-Z)
        Degree=1.0d0
        E=2.718281828d0
C
        v1=80.0701060280091 - 152.3561227776055*(x-y)**2 - 
     -  71.04287940624135*(x-y)**4-4.386716923011888*(x-y)**6+
     -  6.262482097058905*(-1.3378d0+0.5d0*(x+y))+
     -  839.377713504004*(x-y)**2*(-1.3378d0+0.5d0*(x+y))+
     -  220.3350879282827*(x-y)**4*(-1.3378d0+0.5d0*(x+y)) - 
     -  697.6983929689319*(-1.3378d0+0.5d0*(x+y))**2 - 
     -  1528.068314013101*(x-y)**2*(-1.3378d0+0.5d0*(x+y))**2 - 
     -  204.2051289515273*(x-y)**4*(-1.3378d0+0.5d0*(x+y))**2+
     -  1149.37293339893*(-1.3378d0+0.5d0*(x+y))**3+
     -  1117.112519040496*(x-y)**2*(-1.3378d0+0.5d0*(x+y))**3 - 
     -  967.272374980957*(-1.3378d0+0.5d0*(x+y))**4+
     -  929.268634882678*(-1.3378d0+0.5d0*(x+y))**5 - 
     -  860.77528444627*(-1.3378d0+0.5d0*(x+y))**6 - 
     -  63.690468391533*(0.04032362334446276+Cos(Degree*z))+
     -  21.28016866950892*(x-y)**2*(0.04032362334446276+
     -  Cos(Degree*z))+
     -  52.67468265463152*(-1.3378d0+0.5d0*(x+y))*
     -   (0.04032362334446276+Cos(Degree*z))+
     -  49.56942014290817*(x-y)**2*(-1.3378d0+0.5d0*(x+y))*
     -   (0.04032362334446276+Cos(Degree*z)) - 
     -  14.63537131341098*(x-y)**4*(-1.3378d0+0.5d0*(x+y))*
     -   (0.04032362334446276+Cos(Degree*z))+
     -  205.6822442876155*(-1.3378d0+0.5d0*(x+y))**2*
     -   (0.04032362334446276+Cos(Degree*z)) - 
     -  206.3184988394598*(x-y)**2*(-1.3378d0+0.5d0*(x+y))**2*
     -   (0.04032362334446276+Cos(Degree*z))+
     -  23.17264696870488*(-1.3378d0+0.5d0*(x+y))**3*
     -   (0.04032362334446276+Cos(Degree*z)) - 
     -  466.2606077231819*(-1.3378d0+0.5d0*(x+y))**4*
     -   (0.04032362334446276+Cos(Degree*z)) - 
     -  6.252752004329815*(0.04032362334446276+
     -  Cos(Degree*z))**2 - 
     -  31.77034937141531*(x-y)**2*
     -   (0.04032362334446276+Cos(Degree*z))**2 - 
     -  9.30153456744338*(x-y)**4*
     -   (0.04032362334446276+Cos(Degree*z))**2+
     -  53.99390285719871*(-1.3378d0+0.5d0*(x+y))*
     -   (0.04032362334446276+Cos(Degree*z))**2+
     -  52.62078757957153*(x-y)**2*(-1.3378d0+0.5d0*(x+y))*
     -   (0.04032362334446276+Cos(Degree*z))**2 - 
     -  638.0597162619703*(-1.3378d0+0.5d0*(x+y))**2*
     -   (0.04032362334446276+Cos(Degree*z))**2+
     -  528.0975480795673*(-1.3378d0+0.5d0*(x+y))**3*
     -   (0.04032362334446276+Cos(Degree*z))**2 - 
     -  162.2323259917546*(-1.3378d0+0.5d0*(x+y))**4*
     -   (0.04032362334446276+Cos(Degree*z))**2+
     -  112.7284595265815*(0.04032362334446276+
     -  Cos(Degree*z))**3+
     -  43.83199768092012*(x-y)**2*
     -   (0.04032362334446276+Cos(Degree*z))**3+
     -  126.2709420581237*(-1.3378d0+0.5d0*(x+y))*
     -   (0.04032362334446276+Cos(Degree*z))**3 - 
     -  85.9330626944742*(-1.3378d0+0.5d0*(x+y))**2*
     -   (0.04032362334446276+Cos(Degree*z))**3 - 
     -  1303.344922178826*(-1.3378d0+0.5d0*(x+y))**3*
     -   (0.04032362334446276+Cos(Degree*z))**3+
     -  1643.14063062605*(-1.3378d0+0.5d0*(x+y))**4*
     -   (0.04032362334446276+Cos(Degree*z))**3+
     -  357.1392075672478*(0.04032362334446276+
     -  Cos(Degree*z))**4+
     -  110.6437280453673*(x-y)**2*
     -   (0.04032362334446276+Cos(Degree*z))**4 - 
     -  811.593559361352*(-1.3378d0+0.5d0*(x+y))*
     -   (0.04032362334446276+Cos(Degree*z))**4+
     -  2971.778942721558*(-1.3378d0+0.5d0*(x+y))**2*
     -   (0.04032362334446276+Cos(Degree*z))**4 - 
     -  546.7376682764434*(-1.3378d0+0.5d0*(x+y))**3*
     -   (0.04032362334446276+Cos(Degree*z))**4 - 
     -  700.9692326561216*(0.04032362334446276+
     -  Cos(Degree*z))**5 - 
     -  159.009347716965*(-1.3378d0+0.5d0*(x+y))*
     -   (0.04032362334446276+Cos(Degree*z))**5+
     -  874.793680151463*(-1.3378d0+0.5d0*(x+y))**2*
     -   (0.04032362334446276+Cos(Degree*z))**5 - 
     -  1296.827028982462*(0.04032362334446276+
     -  Cos(Degree*z))**6 - 
     -  182.8232076518137*(x-y)**2*
     -   (0.04032362334446276+Cos(Degree*z))**6+
     -  852.000488852307*(-1.3378d0+0.5d0*(x+y))*
     -   (0.04032362334446276+Cos(Degree*z))**6 - 
     -  3526.305585133939*(-1.3378d0+0.5d0*(x+y))**2*
     -   (0.04032362334446276+Cos(Degree*z))**6+
     -  669.5480220734656*(0.04032362334446276+
     -  Cos(Degree*z))**7+
     -  1278.119766773477*(0.04032362334446276+
     -  Cos(Degree*z))**8
        RETURN
        END