!**************************************************************
!* AceGen    5.101 Windows (28 May 13)                        *
!*           Co. J. Korelc  2007            19 Mar 14 10:40:42*
!**************************************************************
! User     : Full professional version
! Notebook : BP_yield_function.nb
! Evaluation time                 : 40 s    Mode  : Optimal
! Number of formulae              : 874     Method: Automatic
! Subroutine                      : BPYieldFunction size :15011
! Total size of Mathematica  code : 15011 subexpressions
! Total size of Fortran code      : 49072 bytes

!******************* S U B R O U T I N E **********************
      SUBROUTINE BPYieldFunction(v,sigma,data,F,DF,DDF,dorder,iter)
      IMPLICIT NONE
      include 'sms.h'
      INTEGER dorder,iter,i01,i02,i15,i34
      LOGICAL b65,b417,b433,b466
      DOUBLE PRECISION v(5005),sigma(6),data(8),F,DF(13),DDF
     &(13,13)
      v(1121)=0.7071067811865476d0*sigma(6)
      v(1039)=(-0.5235987755982988d0)*data(5)
      v(1038)=(-sigma(1)-sigma(2)-sigma(3))/3d0
      v(1037)=data(8)**2
      v(1035)=data(7)
      v(1036)=2d0*v(1035)
      v(1034)=data(6)
      v(1086)=v(1034)**3
      v(1033)=data(4)
      v(1032)=data(3)
      v(1031)=data(2)
      v(1030)=data(1)
      v(1029)=sigma(6)**2
      v(1104)=-v(1029)/2d0
      v(1028)=sigma(5)*v(1121)
      v(1027)=sigma(5)**2
      v(1103)=-v(1027)/2d0
      v(1026)=sigma(4)**2
      v(96)=2d0*v(1030)
      v(39)=(v(1030)*v(1030))
      v(1077)=v(1035)*v(39)
      v(37)=v(1030)+v(1031)
      v(279)=(-2d0)+v(1032)
      v(50)=(-1d0)+v(1032)
      v(45)=1d0-v(1033)
      v(1040)=2d0*v(45)
      v(234)=(v(1034)*v(1034))
      v(264)=v(1036)*v(39)
      v(46)=(v(1035)*v(1035))
      v(1071)=v(1030)*v(46)
      v(1041)=-(v(39)*v(46))
      v(280)=v(1040)*v(1041)
      v(100)=v(1037)*v(96)
      i15=int(dorder)
      v(17)=sigma(1)+v(1038)
      v(18)=sigma(2)+v(1038)
      v(1105)=-v(1026)/2d0+v(17)*v(18)
      v(19)=sigma(3)+v(1038)
      v(20)=(v(1026)+v(1027)+v(1029)+(v(17)*v(17))+(v(18)*v(18))+(v
     &(19)*v(19)))/2d0+v(1037)*v(39)
      v(1106)=sqrt(v(20))
      v(449)=1d0/v(1106)**3
      v(439)=0.25980762113533157d1*v(449)
      v(437)=3d0*v(20)
      v(24)=sigma(4)*v(1028)+v(1103)*v(17)+v(1104)*v(18)+v(1105)*v(19
     &)
      v(25)=sqrt(v(437))
      v(1081)=(v(25)*v(25))
      v(26)=v(37)/2d0
      v(27)=-v(1031)+v(26)
      v(31)=v(1038)-v(27)
      v(133)=2d0*v(31)
      v(74)=(v(31)*v(31))
      v(185)=v(1081)+v(74)
      v(28)=sqrt(v(437)+v(74))
      v(29)=v(24)*v(439)
      v(1102)=v(29)**3
      v(1072)=v(1034)*v(29)
      v(242)=(v(29)*v(29))
      v(1091)=-(v(234)*v(242))
      v(232)=1d0+v(1091)
      v(168)=1d0/Sqrt(v(232))
      v(1090)=(2d0/3d0)*v(168)
      v(82)=v(1039)+dacos(v(1072))/3d0
      v(178)=dtan(v(82))
      v(30)=1d0/dcos(v(82))
      v(236)=1d0/v(30)**2
      v(1044)=2d0*v(236)
      v(32)=datan2(v(25),v(31))
      v(43)=dsin(v(32))
      v(1045)=v(1044)*(v(43)*v(43))
      v(41)=dcos(v(32))
      v(42)=v(41)/v(37)
      v(57)=v(1040)*v(42)
      v(33)=v(26)
      DO i34=1,15
       v(1054)=v(33)*v(41)
       v(71)=v(33)
       v(76)=v(1031)+v(1054)+v(27)
       v(36)=v(33)*v(43)
       v(307)=0.20943951023931953d1*v(36)
       v(295)=(4d0/3d0)*v(168)*v(36)
       v(181)=v(178)*v(295)
       v(83)=(v(36)*v(36))
       v(69)=v(1044)*v(36)
       v(38)=v(76)/v(37)
       v(1042)=2d0*v(38)
       v(151)=1d0-v(1042)
       v(152)=v(1041)*v(151)
       v(149)=v(1042)*v(45)
       v(53)=dsign(1.d0,v(38))
       v(1063)=v(38)*v(53)
       v(67)=-(v(1063)*v(50))
       v(54)=v(42)*v(53)
       v(49)=dabs(v(38))
       v(278)=v(49)**((-3d0)+v(1032))
       v(147)=v(49)**v(279)
       v(1048)=v(147)*v(54)
       v(146)=dlog(v(49))
       v(136)=v(278)*v(279)*v(67)
       v(56)=v(1048)*v(50)
       v(48)=v(49)**v(50)
       v(1043)=(-1d0)+v(48)
       v(281)=-(v(1043)*v(42))-v(38)*v(56)
       v(85)=-(v(1043)*v(38))
       v(58)=v(1041)*v(281)
       v(47)=v(1033)+v(149)
       v(66)=v(1041)*v(47)
       v(1049)=v(42)*v(66)
       v(44)=v(1041)*v(85)
       v(40)=v(44)*v(47)+v(236)*v(83)
       v(52)=v(44)*v(57)+v(47)*v(58)+v(43)*v(69)
       v(63)=v(1045)+v(1049)*(v(136)*v(54)-2d0*v(56))+2d0*v(57)*v(58)
       v(64)=-(v(40)/v(52))
       v(33)=v(33)+v(64)
       IF(dabs(v(64)).lt.0.1d-9.and.dabs(v(40)).lt.0.1d-7) THEN
        v(1097)=v(38)*v(66)
        v(1079)=-(v(43)/v(37))
        v(1078)=(-2d0)*(v(281)*v(47)+v(57)*v(85))
        v(1075)=(-2d0)*v(44)
        v(1074)=v(31)/v(185)
        v(1073)=2d0*v(47)
        v(1068)=v(41)*v(69)
        v(1064)=1d0/v(37)**3
        v(1062)=v(146)*v(48)
        v(1061)=v(147)*v(50)
        v(1067)=((-2d0)*v(1061)+v(136))*v(66)
        v(1059)=v(1074)*v(41)
        v(1058)=v(25)/v(185)
        v(1057)=v(1043)-v(147)*v(67)
        v(1055)=v(1054)*v(236)
        v(1056)=v(1055)*v(181)
        v(1053)=1d0/v(185)**2
        v(1052)=1d0/v(232)**0.15d1
        v(1050)=dsin(v(82))
        v(1047)=1d0/v(37)**2
        v(1046)=dcos(v(82))
        v(1051)=v(1046)*v(1050)
        v(275)=-(v(1057)*v(66))+v(280)*v(85)
        v(287)=v(275)/v(37)
        v(272)=-(v(1041)*v(149))
        v(273)=v(42)*(v(272)-v(66))
        v(284)=v(280)*v(281)+(v(1049)*v(136)+v(1061)*v(273))*v(53)-v
     &  (1041)*v(1043)*v(57)-v(56)*v(66)
        v(285)=v(284)/v(37)
        v(288)=2d0*v(1068)-v(43)*(v(287)+v(285)*v(71))
        v(293)=v(1058)*v(288)
        v(1080)=2d0*v(293)
        v(290)=v(285)+v(293)
        v(292)=(v(290)-2d0*v(1047)*(v(275)*v(41)+v(284)*v(76)))/2d0
        v(268)=0.3490658503988659d0*v(168)*v(83)
        v(253)=-(v(1071)*v(85))
        v(249)=4d0*v(253)*v(45)
        v(243)=((2d0/9d0)*v(83))/v(232)
        v(306)=v(1051)*v(43)
        v(1076)=v(295)*v(306)
        v(238)=(v(1050)*v(1050))
        v(1088)=-v(236)+v(238)
        v(231)=v(1051)*v(83)
        v(241)=(2d0/3d0)*v(1052)*v(231)
        v(183)=v(1056)*v(29)
        v(180)=v(1055)*v(178)*v(307)
        v(170)=v(1034)*v(1056)
        v(166)=v(1044)*v(1054)
        v(156)=v(264)*(v(1057)*v(47)-v(1040)*v(85))
        v(164)=v(156)/v(37)
        v(184)=-(v(164)*v(36))
        v(153)=v(1075)-v(1057)*v(152)
        v(163)=v(153)/v(37)
        v(177)=-(v(163)*v(36))
        v(138)=-(v(1057)*v(280))
        v(1066)=2d0*v(138)
        v(251)=v(1062)*v(38)
        v(121)=-(v(1047)*v(76))
        v(104)=v(1058)/2d0
        v(105)=-(v(1058)*v(41))
        v(107)=v(104)*v(41)
        v(1065)=v(107)*v(71)
        v(108)=v(1058)*v(43)
        v(109)=-(v(104)*v(43))
        v(1060)=v(109)*v(71)
        v(173)=v(1060)*v(69)
        v(111)=v(105)*v(71)
        v(112)=0.5d0+v(1060)
        v(114)=0.5d0-v(1060)
        v(174)=v(1065)*v(166)
        v(117)=(v(108)*v(71))/v(37)
        v(118)=v(111)/v(37)
        v(120)=v(121)+v(112)/v(37)
        v(122)=v(121)+v(114)/v(37)
        v(123)=v(117)*v(53)
        v(124)=v(118)*v(53)
        v(125)=v(120)*v(53)
        v(126)=v(122)*v(53)
        v(129)=v(1061)*v(125)
        v(246)=-(v(1043)*v(120))-v(129)*v(38)
        v(142)=((-2d0)*v(129)+v(125)*v(136))*v(66)-v(1057)*(v(120)*v
     &  (280)-v(46)*v(47)*v(96))-v(1040)*v(46)*(v(246)*v(39)+v(85)*v
     &  (96))
        v(144)=v(1066)*v(122)+v(1067)*v(126)
        v(256)=v(144)/v(37)
        v(150)=v(1062)*v(272)-v(66)*(v(1062)+v(147)*(v(1063)-v(146)*v
     &  (67)))
        v(162)=v(150)/v(37)
        v(176)=-(v(162)*v(36))
        v(94)=v(1073)*v(85)
        v(90)=-(v(231)/v(52))
        v(92)=v(1090)*v(90)
        v(206)=v(1058)*v(184)
        v(207)=v(164)+v(206)
        v(218)=v(121)*v(156)+v(207)/2d0
        v(205)=v(1058)*v(183)
        v(217)=v(205)/2d0
        v(204)=v(1058)*v(180)
        v(216)=v(204)/2d0
        v(202)=v(1058)*v(177)
        v(203)=v(163)+v(202)
        v(215)=v(121)*v(153)+v(203)/2d0
        v(200)=v(1058)*v(176)
        v(201)=v(162)+v(200)
        v(214)=v(121)*v(150)+v(201)/2d0
        v(194)=v(1058)*v(170)
        v(210)=v(194)/2d0
        v(212)=2d0*v(1064)*v(275)*v(76)
        v(160)=-(v(1047)*v(275))
        v(161)=v(160)+v(256)
        v(159)=v(160)+v(142)/v(37)
        v(172)=-(v(1065)*v(287))
        v(175)=-v(172)-v(173)-v(174)-v(161)*v(36)
        v(171)=v(172)+v(173)+v(174)-v(159)*v(36)
        v(167)=(v(1079)*(v(1066)*v(118)+v(1067)*v(124))+v(1059)*(v
     &  (166)-v(287))+v(105)*v(69))*v(71)
        v(72)=(v(1068)-v(287)*v(43))*v(71)
        v(1069)=v(1053)*v(72)
        v(1070)=v(1069)*v(25)
        v(224)=v(1069)*v(74)
        v(223)=-v(72)/(2d0*v(185))
        v(1085)=v(223)+v(224)
        v(219)=-(v(1070)*v(133))
        v(197)=v(1070)*v(31)
        v(198)=v(1058)*v(175)-v(197)
        v(199)=v(161)+v(198)
        v(213)=v(199)/2d0+v(212)-v(1047)*(v(114)*v(275)+v(144)*v(76))
        v(195)=v(1058)*v(171)+v(197)
        v(78)=(-2d0)*v(223)*v(25)
        v(75)=v(287)+v(78)
        v(77)=v(75)/2d0+v(160)*v(76)
        v(79)=v(78)/v(52)
        v(1082)=-(v(63)*v(79))
        v(1083)=v(1082)+v(293)
        v(81)=(v(133)*v(223))/v(52)
        v(1084)=v(63)*v(81)
        v(84)=v(1034)*v(92)
        v(1087)=v(63)*v(84)
        v(86)=(-v(77)+v(1071)*v(94))/v(52)
        v(1093)=v(63)*v(86)
        v(87)=-((v(287)-v(75)+v(77))/v(52))
        v(1095)=v(63)*v(87)
        v(88)=(v(251)*v(66))/v(52)
        v(1098)=v(63)*v(88)
        v(89)=-((v(151)*v(44))/v(52))
        v(1099)=v(63)*v(89)
        v(91)=0.10471975511965977d1*v(90)
        v(1101)=v(63)*v(91)
        v(93)=v(29)*v(92)
        v(95)=(v(1077)*v(94))/v(52)
        v(235)=v(234)*v(243)
        v(239)=v(1034)*v(268)
        v(244)=v(1072)*v(243)
        v(247)=-(v(1073)*v(46))
        v(1094)=-(v(1030)*v(247))
        v(266)=0.5483113556160755d0*v(83)
        v(294)=v(1074)*v(288)
        v(297)=v(1034)*v(1076)
        v(1089)=2d0*v(297)
        v(1092)=-v(1087)-v(1089)
        v(300)=v(1071)*v(1078)+v(292)
        v(301)=v(285)-v(290)+v(292)
        v(1096)=-v(1095)-v(301)
        v(304)=v(1062)*v(273)+v(1097)*(-v(1048)-v(146)*v(56))
        v(305)=v(1075)*v(42)+v(151)*v(58)
        v(308)=v(306)*v(307)
        v(1100)=(-2d0)*v(308)
        v(309)=v(1076)*v(29)
        v(310)=v(1077)*v(1078)
        v(311)=(v(219)+v(1058)*(v(111)*v(166)+(v(1079)*(v(1066)*v(117
     &  )+v(1067)*v(123))-v(105)*v(287)+v(108)*v(69))*v(71))+v(1080
     &  )*v(79)-v(63)*(v(79)*v(79)))/v(52)
        v(313)=(v(1058)*v(167)+((-2d0)*v(1053)*v(1081)+1d0/v(185))*v
     &  (72)+(v(1080)+v(1082))*v(81))/v(52)
        v(314)=(v(194)-v(297)*v(79)+v(1083)*v(84))/v(52)
        v(315)=(v(195)-v(300)*v(79)+v(1083)*v(86))/v(52)
        v(316)=(v(198)-v(301)*v(79)+v(1083)*v(87))/v(52)
        v(317)=(v(200)-v(304)*v(79)+v(1083)*v(88))/v(52)
        v(318)=(v(202)-v(305)*v(79)+v(1083)*v(89))/v(52)
        v(319)=(v(204)-v(308)*v(79)+v(1083)*v(91))/v(52)
        v(320)=(v(205)-v(309)*v(79)+v(1083)*v(93))/v(52)
        v(321)=(v(206)-v(310)*v(79)+v(1083)*v(95))/v(52)
        v(322)=v(313)
        v(323)=-((v(1074)*v(167)+v(219)+2d0*v(294)*v(81)+v(63)*(v(81
     &  )*v(81)))/v(52))
        v(324)=-((v(1074)*v(170)+v(297)*v(81)+v(1084)*v(84)+v(294)*v
     &  (84))/v(52))
        v(325)=-((v(171)*v(31)+v(185)*(v(1085)+v(300)*v(81)+v(1084)*v
     &  (86)+v(294)*v(86)))/(v(185)*v(52)))
        v(326)=(v(1085)-v(1074)*v(175)-v(301)*v(81)+(-v(1084)-v(294)
     &  )*v(87))/v(52)
        v(327)=-((v(1074)*v(176)+v(304)*v(81)+v(1084)*v(88)+v(294)*v
     &  (88))/v(52))
        v(328)=-((v(1074)*v(177)+v(305)*v(81)+v(1084)*v(89)+v(294)*v
     &  (89))/v(52))
        v(329)=-((v(1074)*v(180)+v(308)*v(81)+v(1084)*v(91)+v(294)*v
     &  (91))/v(52))
        v(330)=-((v(1074)*v(183)+v(309)*v(81)+v(1084)*v(93)+v(294)*v
     &  (93))/v(52))
        v(331)=-((v(1074)*v(184)+v(310)*v(81)+v(1084)*v(95)+v(294)*v
     &  (95))/v(52))
        v(332)=v(314)
        v(333)=v(324)
        v(334)=-((v(1088)*v(235)+v(1086)*v(241)*v(29)+v(1089)*v(84)+v
     &  (63)*(v(84)*v(84)))/v(52))
        v(335)=-((v(210)+v(300)*v(84)+v(1087)*v(86)+v(297)*v(86))/v
     &  (52))
        v(336)=(v(210)-v(301)*v(84)+(-v(1087)-v(297))*v(87))/v(52)
        v(337)=-((v(304)*v(84)+v(1087)*v(88)+v(297)*v(88))/v(52))
        v(338)=-((v(305)*v(84)+v(1087)*v(89)+v(297)*v(89))/v(52))
        v(340)=(-(v(1088)*v(239))+v(1092)*v(91))/v(52)
        v(342)=(-(v(1090)*v(231))+v(1091)*v(241)-v(1088)*v(244)+v
     &  (1092)*v(93))/v(52)
        v(343)=-((v(310)*v(84)+v(1087)*v(95)+v(297)*v(95))/v(52))
        v(344)=v(315)
        v(345)=v(325)
        v(346)=v(335)
        v(347)=-(2d0*v(121)*v(142)+v(159)+2d0*v(112)*v(160)+v(195)
     &  +2d0*v(212)+2d0*v(120)*v(249)+2d0*v(247)*v(85)+4d0*v(300)*v
     &  (86)+2d0*v(63)*(v(86)*v(86))+v(246)*v(247)*v(96))/(2d0*v(52)
     &  )
        v(348)=-((v(1057)*v(1094)*v(122)+v(213)+v(122)*v(249)+v(301
     &  )*v(86)+v(1093)*v(87)+v(300)*v(87))/v(52))
        v(349)=-((v(214)+v(1094)*v(251)+v(304)*v(86)+v(1093)*v(88)+v
     &  (300)*v(88))/v(52))
        v(350)=-((v(215)+2d0*v(151)*v(253)+v(305)*v(86)+v(1093)*v(89)
     &  +v(300)*v(89))/v(52))
        v(351)=-((v(216)+v(308)*v(86)+v(1093)*v(91)+v(300)*v(91))/v
     &  (52))
        v(352)=-((v(217)+v(309)*v(86)+v(1093)*v(93)+v(300)*v(93))/v
     &  (52))
        v(353)=-((v(218)-4d0*v(1030)*v(1035)*v(47)*v(85)+v(310)*v(86)
     &  +v(1093)*v(95)+v(300)*v(95))/v(52))
        v(354)=v(316)
        v(355)=v(326)
        v(356)=v(336)
        v(357)=v(348)
        v(358)=-((v(161)-v(199)+v(213)+2d0*v(301)*v(87)+v(63)*(v(87
     &  )*v(87)))/v(52))
        v(359)=-((v(162)-v(201)+v(214)+v(304)*v(87)+v(1095)*v(88)+v
     &  (301)*v(88))/v(52))
        v(360)=-((v(163)-v(203)+v(215)+v(305)*v(87)+v(1095)*v(89)+v
     &  (301)*v(89))/v(52))
        v(361)=(v(216)-v(308)*v(87)+v(1096)*v(91))/v(52)
        v(362)=(v(217)-v(309)*v(87)+v(1096)*v(93))/v(52)
        v(363)=-((v(164)-v(207)+v(218)+v(310)*v(87)+v(1095)*v(95)+v
     &  (301)*v(95))/v(52))
        v(364)=v(317)
        v(365)=v(327)
        v(366)=v(337)
        v(367)=v(349)
        v(368)=v(359)
        v(369)=(v(1097)*(v(146)*v(146))*v(48)-2d0*v(304)*v(88)-v(63)*
     &  (v(88)*v(88)))/v(52)
        v(370)=(v(152)*v(251)-v(305)*v(88)+(-v(1098)-v(304))*v(89))/v
     &  (52)
        v(371)=-((v(308)*v(88)+v(1098)*v(91)+v(304)*v(91))/v(52))
        v(372)=-((v(309)*v(88)+v(1098)*v(93)+v(304)*v(93))/v(52))
        v(373)=-((v(251)*v(264)*v(47)+v(310)*v(88)+v(1098)*v(95)+v
     &  (304)*v(95))/v(52))
        v(374)=v(318)
        v(375)=v(328)
        v(376)=v(338)
        v(377)=v(350)
        v(378)=v(360)
        v(379)=v(370)
        v(380)=-(((v(1099)+2d0*v(305))*v(89))/v(52))
        v(381)=-((v(308)*v(89)+v(1099)*v(91)+v(305)*v(91))/v(52))
        v(382)=-((v(309)*v(89)+v(1099)*v(93)+v(305)*v(93))/v(52))
        v(383)=(v(151)*v(264)*v(85)-v(310)*v(89)+(-v(1099)-v(305))*v
     &  (95))/v(52)
        v(384)=v(319)
        v(385)=v(329)
        v(386)=v(340)
        v(387)=v(351)
        v(388)=v(361)
        v(389)=v(371)
        v(390)=v(381)
        v(391)=(-(v(1088)*v(266))+v(1100)*v(91)-v(63)*(v(91)*v(91)))
     &  /v(52)
        v(393)=(-(v(1088)*v(268)*v(29))+(v(1100)-v(1101))*v(93))/v(52
     &  )
        v(394)=-((v(310)*v(91)+v(1101)*v(95)+v(308)*v(95))/v(52))
        v(395)=v(320)
        v(396)=v(330)
        v(397)=v(342)
        v(398)=v(352)
        v(399)=v(362)
        v(400)=v(372)
        v(401)=v(382)
        v(402)=v(393)
        v(403)=(-(v(1034)*v(1102)*v(241))-v(1088)*v(242)*v(243)-2d0*v
     &  (309)*v(93)-v(63)*(v(93)*v(93)))/v(52)
        v(404)=-((v(310)*v(93)+v(309)*v(95)+v(63)*v(93)*v(95))/v(52))
        v(405)=v(321)
        v(406)=v(331)
        v(407)=v(343)
        v(408)=v(353)
        v(409)=v(363)
        v(410)=v(373)
        v(411)=v(383)
        v(412)=v(394)
        v(413)=v(404)
        v(414)=(v(39)*v(94)-2d0*v(310)*v(95)-v(63)*(v(95)*v(95)))/v
     &  (52)
        iter=i34
        exit
       ELSE
       ENDIF
       IF(i34.eq.15) THEN
        iter=(-1d0)
        exit
       ELSE
       ENDIF
      ENDDO
      v(420)=v(79)
      v(421)=v(81)
      v(422)=v(84)
      v(423)=v(86)
      v(424)=v(87)
      v(425)=v(88)
      v(426)=v(89)
      v(427)=v(91)
      v(428)=v(93)
      v(429)=v(95)
      v(430)=v(33)
      F=(-1d0)+v(28)/v(430)
      IF(i15.ge.1) THEN
       v(1123)=0.15d1/v(25)
       v(1122)=sigma(4)*sigma(5)
       v(1114)=2d0*sigma(6)
       v(1113)=2d0*sigma(5)
       v(1112)=2d0*sigma(4)
       v(1111)=v(17)*v(19)
       v(1110)=1d0/(2d0*v(28))
       v(1108)=1d0/(v(28)*v(430))
       v(802)=v(1103)+v(18)*v(19)
       v(785)=v(1104)+v(1111)
       v(736)=-(v(1123)*v(28))
       v(735)=v(421)*v(736)
       v(477)=1d0/v(20)**2
       v(734)=((-0.3897114317029974d1)*v(477))/v(1106)
       v(733)=v(24)*v(734)
       v(529)=v(1110)*v(31)
       v(458)=-(v(529)/v(430))
       v(436)=1d0/v(430)**2
       v(1109)=v(28)*v(436)
       v(434)=-(v(1109)*v(422))
       v(1107)=v(434)*v(439)
       v(705)=(-0.12990381056766578d1)*v(434)
       v(690)=0.18371173070873836d1*v(434)
       v(455)=v(449)*v(690)
       v(452)=v(449)*v(705)
       v(451)=sigma(6)*v(455)
       v(438)=0.15d1*v(1108)+v(434)*v(733)+v(436)*v(735)
       v(453)=v(438)/2d0
       v(906)=v(18)*v(452)+v(453)
       v(891)=v(17)*v(452)+v(453)
       v(876)=v(19)*v(452)+v(453)
       v(440)=v(1105)*v(1107)+v(19)*v(438)
       v(441)=v(18)*v(438)+v(1107)*v(785)
       v(442)=v(17)*v(438)+v(1107)*v(802)
       v(446)=(-(v(1108)*v(31))+v(1109)*v(420)-v(440)-v(441)-v(442))
     & /3d0
       DF(1)=v(442)+v(446)
       DF(2)=v(441)+v(446)
       DF(3)=v(440)+v(446)
       DF(4)=sigma(5)*v(451)+v(1112)*v(876)
       DF(5)=sigma(4)*v(451)+v(1113)*v(891)
       DF(6)=v(1122)*v(455)+v(1114)*v(906)
       DF(7)=-(v(1109)*v(423))+v(100)*v(438)+v(458)
       DF(8)=-(v(1109)*v(424))-v(458)
       DF(9)=-(v(1109)*v(425))
       DF(10)=-(v(1109)*v(426))
       DF(11)=-(v(1109)*v(427))
       DF(12)=-(v(1109)*v(428))
       DF(13)=-(v(1109)*v(429))
       IF(i15.eq.2) THEN
        v(1160)=-(v(28)*v(429))
        v(1151)=-(v(436)*v(529))
        v(1150)=-(v(31)/v(28))
        v(1142)=v(439)*v(802)
        v(1141)=v(439)*v(785)
        v(1133)=v(1105)*v(439)
        v(1131)=v(434)*v(734)
        v(1130)=0.15d1/v(430)
        v(1129)=0.15d1/v(28)
        v(1117)=3d0*sigma(6)
        v(1116)=3d0*sigma(5)
        v(1115)=3d0*sigma(4)
        v(851)=(2d0/3d0)*v(452)
        v(836)=-v(452)/3d0
        v(824)=v(1110)/v(430)
        v(817)=v(1108)/3d0
        v(789)=-v(19)/3d0
        v(787)=(2d0/3d0)*v(438)
        v(773)=-v(18)/3d0
        v(772)=-v(17)/3d0
        v(770)=-v(438)/3d0
        v(741)=(0.75d0*v(28))/v(437)
        v(740)=(0.5625d1*v(434))/v(20)**3
        v(739)=(0.17320508075688773d1*v(24))/v(1106)
        v(1132)=v(739)*v(740)
        v(707)=(-0.12990381056766578d1)*v(449)
        v(692)=0.18371173070873836d1*v(449)
        v(664)=1d0/v(430)**3
        v(1124)=(-2d0)*v(664)
        v(656)=v(1108)/4d0
        v(649)=-v(31)/(2d0*v(430))
        v(648)=v(1108)/6d0
        v(639)=1d0/v(28)**2
        v(618)=-v(420)/3d0
        v(610)=-v(405)/3d0
        v(602)=-v(395)/3d0
        v(594)=-v(384)/3d0
        v(586)=-v(374)/3d0
        v(578)=-v(364)/3d0
        v(570)=-v(354)/3d0
        v(562)=-v(344)/3d0
        v(554)=-v(332)/3d0
        v(546)=-v(322)/3d0
        v(538)=-v(311)/3d0
        v(522)=(-2d0/3d0)*v(31)
        v(507)=v(1029)/6d0
        v(506)=v(1027)/6d0
        v(505)=v(19)*v(773)
        v(1119)=v(505)+v(506)
        v(504)=-v(1111)/3d0
        v(1120)=v(504)+v(507)
        v(502)=v(1026)/6d0
        v(501)=v(18)*v(772)
        v(1118)=v(501)+v(502)
        v(479)=0.3849001794597505d0*v(734)
        v(474)=(-2d0/3d0)*v(18)
        v(473)=(-2d0/3d0)*v(17)
        v(471)=(-2d0/3d0)*v(19)
        v(478)=v(17)*v(479)
        v(480)=v(18)*v(479)
        v(481)=v(19)*v(479)
        v(482)=sigma(4)*v(479)
        v(483)=sigma(5)*v(479)
        v(484)=sigma(6)*v(479)
        v(485)=v(100)*v(479)
        v(486)=0.25980762113533157d1*v(478)
        v(487)=0.25980762113533157d1*v(480)
        v(488)=0.25980762113533157d1*v(481)
        v(489)=0.25980762113533157d1*v(482)
        v(490)=0.25980762113533157d1*v(483)
        v(491)=0.25980762113533157d1*v(484)
        v(492)=0.25980762113533157d1*v(485)
        v(493)=3d0*v(17)
        v(494)=3d0*v(18)
        v(495)=3d0*v(19)
        v(499)=3d0*v(100)
        v(500)=-v(1027)/3d0+v(1118)+v(1120)-v(18)*v(471)
        v(503)=-v(1029)/3d0+v(1118)+v(1119)-v(17)*v(471)
        v(508)=-v(1026)/3d0+v(1119)+v(1120)-v(18)*v(473)
        v(510)=v(1028)-(v(1112)*v(19))/2d0
        v(511)=sigma(4)*v(1121)-(v(1113)*v(17))/2d0
        v(512)=0.7071067811865476d0*v(1122)-(v(1114)*v(18))/2d0
        v(514)=v(1123)*v(17)
        v(515)=v(1123)*v(18)
        v(516)=v(1123)*v(19)
        v(517)=sigma(4)*v(1123)
        v(518)=sigma(5)*v(1123)
        v(519)=sigma(6)*v(1123)
        v(520)=v(100)*v(1123)
        v(521)=v(1110)*(v(493)+v(522))
        v(523)=v(1110)*(v(494)+v(522))
        v(524)=v(1110)*(v(495)+v(522))
        v(525)=v(1110)*v(1115)
        v(526)=v(1110)*v(1116)
        v(527)=v(1110)*v(1117)
        v(528)=v(1110)*(-v(31)+v(499))
        v(530)=v(24)*v(486)+v(439)*v(500)
        v(531)=v(24)*v(487)+v(439)*v(503)
        v(532)=v(24)*v(488)+v(439)*v(508)
        v(533)=v(24)*v(489)+v(439)*v(510)
        v(534)=v(24)*v(490)+v(439)*v(511)
        v(535)=v(24)*v(491)+v(439)*v(512)
        v(536)=v(24)*v(492)
        v(537)=v(313)*v(514)+v(314)*v(530)+v(538)
        v(539)=v(313)*v(515)+v(314)*v(531)+v(538)
        v(540)=v(313)*v(516)+v(314)*v(532)+v(538)
        v(541)=v(313)*v(517)+v(314)*v(533)
        v(542)=v(313)*v(518)+v(314)*v(534)
        v(543)=v(313)*v(519)+v(314)*v(535)
        v(544)=v(315)+v(313)*v(520)+v(314)*v(536)
        v(545)=v(323)*v(514)+v(324)*v(530)+v(546)
        v(547)=v(323)*v(515)+v(324)*v(531)+v(546)
        v(548)=v(323)*v(516)+v(324)*v(532)+v(546)
        v(549)=v(323)*v(517)+v(324)*v(533)
        v(550)=v(323)*v(518)+v(324)*v(534)
        v(551)=v(323)*v(519)+v(324)*v(535)
        v(552)=v(325)+v(323)*v(520)+v(324)*v(536)
        v(553)=v(333)*v(514)+v(334)*v(530)+v(554)
        v(555)=v(333)*v(515)+v(334)*v(531)+v(554)
        v(556)=v(333)*v(516)+v(334)*v(532)+v(554)
        v(557)=v(333)*v(517)+v(334)*v(533)
        v(558)=v(333)*v(518)+v(334)*v(534)
        v(559)=v(333)*v(519)+v(334)*v(535)
        v(560)=v(335)+v(333)*v(520)+v(334)*v(536)
        v(561)=v(345)*v(514)+v(346)*v(530)+v(562)
        v(563)=v(345)*v(515)+v(346)*v(531)+v(562)
        v(564)=v(345)*v(516)+v(346)*v(532)+v(562)
        v(565)=v(345)*v(517)+v(346)*v(533)
        v(566)=v(345)*v(518)+v(346)*v(534)
        v(567)=v(345)*v(519)+v(346)*v(535)
        v(568)=v(347)+v(345)*v(520)+v(346)*v(536)
        v(569)=v(355)*v(514)+v(356)*v(530)+v(570)
        v(571)=v(355)*v(515)+v(356)*v(531)+v(570)
        v(572)=v(355)*v(516)+v(356)*v(532)+v(570)
        v(573)=v(355)*v(517)+v(356)*v(533)
        v(574)=v(355)*v(518)+v(356)*v(534)
        v(575)=v(355)*v(519)+v(356)*v(535)
        v(576)=v(357)+v(355)*v(520)+v(356)*v(536)
        v(577)=v(365)*v(514)+v(366)*v(530)+v(578)
        v(579)=v(365)*v(515)+v(366)*v(531)+v(578)
        v(580)=v(365)*v(516)+v(366)*v(532)+v(578)
        v(581)=v(365)*v(517)+v(366)*v(533)
        v(582)=v(365)*v(518)+v(366)*v(534)
        v(583)=v(365)*v(519)+v(366)*v(535)
        v(584)=v(367)+v(365)*v(520)+v(366)*v(536)
        v(585)=v(375)*v(514)+v(376)*v(530)+v(586)
        v(587)=v(375)*v(515)+v(376)*v(531)+v(586)
        v(588)=v(375)*v(516)+v(376)*v(532)+v(586)
        v(589)=v(375)*v(517)+v(376)*v(533)
        v(590)=v(375)*v(518)+v(376)*v(534)
        v(591)=v(375)*v(519)+v(376)*v(535)
        v(592)=v(377)+v(375)*v(520)+v(376)*v(536)
        v(593)=v(385)*v(514)+v(386)*v(530)+v(594)
        v(595)=v(385)*v(515)+v(386)*v(531)+v(594)
        v(596)=v(385)*v(516)+v(386)*v(532)+v(594)
        v(597)=v(385)*v(517)+v(386)*v(533)
        v(598)=v(385)*v(518)+v(386)*v(534)
        v(599)=v(385)*v(519)+v(386)*v(535)
        v(600)=v(387)+v(385)*v(520)+v(386)*v(536)
        v(601)=v(396)*v(514)+v(397)*v(530)+v(602)
        v(603)=v(396)*v(515)+v(397)*v(531)+v(602)
        v(604)=v(396)*v(516)+v(397)*v(532)+v(602)
        v(605)=v(396)*v(517)+v(397)*v(533)
        v(606)=v(396)*v(518)+v(397)*v(534)
        v(607)=v(396)*v(519)+v(397)*v(535)
        v(608)=v(398)+v(396)*v(520)+v(397)*v(536)
        v(609)=v(406)*v(514)+v(407)*v(530)+v(610)
        v(611)=v(406)*v(515)+v(407)*v(531)+v(610)
        v(612)=v(406)*v(516)+v(407)*v(532)+v(610)
        v(613)=v(406)*v(517)+v(407)*v(533)
        v(614)=v(406)*v(518)+v(407)*v(534)
        v(615)=v(406)*v(519)+v(407)*v(535)
        v(616)=v(408)+v(406)*v(520)+v(407)*v(536)
        v(617)=v(421)*v(514)+v(422)*v(530)+v(618)
        v(619)=v(421)*v(515)+v(422)*v(531)+v(618)
        v(620)=v(421)*v(516)+v(422)*v(532)+v(618)
        v(621)=v(421)*v(517)+v(422)*v(533)
        v(622)=v(421)*v(518)+v(422)*v(534)
        v(623)=v(421)*v(519)+v(422)*v(535)
        v(624)=v(423)+v(421)*v(520)+v(422)*v(536)
        v(625)=-(v(436)*v(617))
        v(626)=-(v(436)*v(619))
        v(627)=-(v(436)*v(620))
        v(628)=-(v(436)*v(621))
        v(629)=-(v(436)*v(622))
        v(630)=-(v(436)*v(623))
        v(631)=-(v(436)*v(624))
        v(632)=-(v(424)*v(436))
        v(935)=-(v(529)*v(632))
        v(633)=-(v(425)*v(436))
        v(634)=-(v(426)*v(436))
        v(635)=-(v(427)*v(436))
        v(636)=-(v(428)*v(436))
        v(637)=-(v(429)*v(436))
        v(638)=-(v(521)*v(639))
        v(640)=-(v(523)*v(639))
        v(641)=-(v(524)*v(639))
        v(642)=-(v(525)*v(639))
        v(643)=-(v(526)*v(639))
        v(644)=-(v(527)*v(639))
        v(645)=-(v(528)*v(639))
        v(646)=-(v(529)*v(639))
        v(647)=-(v(529)*v(625))+v(648)+v(638)*v(649)
        v(650)=-(v(529)*v(626))+v(648)+v(640)*v(649)
        v(651)=-(v(529)*v(627))+v(648)+v(641)*v(649)
        v(652)=-(v(529)*v(628))+v(642)*v(649)
        v(653)=-(v(529)*v(629))+v(643)*v(649)
        v(654)=-(v(529)*v(630))+v(644)*v(649)
        v(655)=-(v(529)*v(631))+v(645)*v(649)+v(656)
        v(657)=v(646)*v(649)-v(656)+v(935)
        v(658)=-(v(529)*v(633))
        v(659)=-(v(529)*v(634))
        v(660)=-(v(529)*v(635))
        v(661)=-(v(529)*v(636))
        v(662)=-(v(529)*v(637))
        v(663)=v(1124)*v(617)
        v(1143)=v(28)*v(663)
        v(665)=v(1124)*v(619)
        v(1144)=v(28)*v(665)
        v(666)=v(1124)*v(620)
        v(1145)=v(28)*v(666)
        v(667)=v(1124)*v(621)
        v(1146)=v(28)*v(667)
        v(668)=v(1124)*v(622)
        v(1147)=v(28)*v(668)
        v(669)=v(1124)*v(623)
        v(1148)=v(28)*v(669)
        v(670)=v(1124)*v(624)
        v(1149)=v(28)*v(670)
        v(671)=v(1124)*v(424)
        v(1125)=-(v(28)*v(671))
        v(1011)=v(1125)*v(429)
        v(1156)=v(1011)-v(662)
        v(994)=v(1125)*v(428)
        v(1155)=-v(661)+v(994)
        v(978)=v(1125)*v(427)
        v(1154)=-v(660)+v(978)
        v(963)=v(1125)*v(426)
        v(1153)=-v(659)+v(963)
        v(949)=v(1125)*v(425)
        v(1152)=-v(658)+v(949)
        v(672)=v(1124)*v(425)
        v(1126)=-(v(28)*v(672))
        v(1013)=v(1126)*v(429)
        v(996)=v(1126)*v(428)
        v(980)=v(1126)*v(427)
        v(965)=v(1126)*v(426)
        v(673)=v(1124)*v(426)
        v(1127)=-(v(28)*v(673))
        v(1015)=v(1127)*v(429)
        v(998)=v(1127)*v(428)
        v(982)=v(1127)*v(427)
        v(674)=v(1124)*v(427)
        v(1128)=-(v(28)*v(674))
        v(1017)=v(1128)*v(429)
        v(1000)=v(1128)*v(428)
        v(675)=v(1124)*v(428)
        v(1019)=v(1160)*v(675)
        v(676)=v(1124)*v(429)
        v(677)=-(v(422)*(v(1143)+v(436)*v(521)))-v(1109)*v(553)
        v(1134)=v(434)*v(486)+v(439)*v(677)
        v(678)=-(v(422)*(v(1144)+v(436)*v(523)))-v(1109)*v(555)
        v(1135)=v(434)*v(487)+v(439)*v(678)
        v(679)=-(v(422)*(v(1145)+v(436)*v(524)))-v(1109)*v(556)
        v(1136)=v(434)*v(488)+v(439)*v(679)
        v(680)=-(v(422)*(v(1146)+v(436)*v(525)))-v(1109)*v(557)
        v(1137)=v(434)*v(489)+v(439)*v(680)
        v(681)=-(v(422)*(v(1147)+v(436)*v(526)))-v(1109)*v(558)
        v(1138)=v(434)*v(490)+v(439)*v(681)
        v(682)=-(v(422)*(v(1148)+v(436)*v(527)))-v(1109)*v(559)
        v(1139)=v(434)*v(491)+v(439)*v(682)
        v(683)=-(v(422)*(v(1149)+v(436)*v(528)))-v(1109)*v(560)
        v(1140)=v(434)*v(492)+v(439)*v(683)
        v(684)=v(1151)*v(422)-v(28)*(v(336)*v(436)+v(422)*v(671))
        v(685)=-(v(28)*(v(337)*v(436)+v(422)*v(672)))
        v(686)=-(v(28)*(v(338)*v(436)+v(422)*v(673)))
        v(687)=-(v(28)*(v(340)*v(436)+v(422)*v(674)))
        v(688)=-(v(28)*(v(342)*v(436)+v(422)*v(675)))
        v(689)=-(v(28)*(v(343)*v(436)+v(422)*v(676)))
        v(691)=v(478)*v(690)+v(677)*v(692)
        v(693)=v(480)*v(690)+v(678)*v(692)
        v(694)=v(481)*v(690)+v(679)*v(692)
        v(695)=v(482)*v(690)+v(680)*v(692)
        v(696)=v(483)*v(690)+v(681)*v(692)
        v(697)=v(484)*v(690)+v(682)*v(692)
        v(698)=v(485)*v(690)+v(683)*v(692)
        v(699)=v(684)*v(692)
        v(700)=v(685)*v(692)
        v(701)=v(686)*v(692)
        v(702)=v(687)*v(692)
        v(703)=v(688)*v(692)
        v(704)=v(689)*v(692)
        v(706)=v(478)*v(705)+v(677)*v(707)
        v(708)=v(480)*v(705)+v(678)*v(707)
        v(709)=v(481)*v(705)+v(679)*v(707)
        v(710)=v(482)*v(705)+v(680)*v(707)
        v(711)=v(483)*v(705)+v(681)*v(707)
        v(712)=v(484)*v(705)+v(682)*v(707)
        v(713)=v(485)*v(705)+v(683)*v(707)
        v(714)=v(684)*v(707)
        v(715)=v(685)*v(707)
        v(716)=v(686)*v(707)
        v(717)=v(687)*v(707)
        v(718)=v(688)*v(707)
        v(719)=v(689)*v(707)
        v(720)=sigma(6)*v(691)
        v(721)=sigma(6)*v(693)
        v(722)=sigma(6)*v(694)
        v(723)=sigma(6)*v(695)
        v(724)=sigma(6)*v(696)
        v(725)=v(455)+sigma(6)*v(697)
        v(726)=sigma(6)*v(698)
        v(727)=sigma(6)*v(699)
        v(728)=sigma(6)*v(700)
        v(729)=sigma(6)*v(701)
        v(730)=sigma(6)*v(702)
        v(731)=sigma(6)*v(703)
        v(732)=sigma(6)*v(704)
        v(737)=(v(421)*v(436))/v(25)
        v(738)=v(1132)*v(17)+v(1131)*v(500)+v(1129)*v(625)+v(1130)*v
     &  (638)+v(677)*v(733)+v(663)*v(735)+v(436)*v(545)*v(736)+v(737
     &  )*((-0.15d1)*v(521)+v(493)*v(741))
        v(742)=v(1132)*v(18)+v(1131)*v(503)+v(1129)*v(626)+v(1130)*v
     &  (640)+v(678)*v(733)+(v(436)*v(547)+v(421)*v(665))*v(736)+v
     &  (737)*((-0.15d1)*v(523)+v(494)*v(741))
        v(744)=v(1132)*v(19)+v(1131)*v(508)+v(1129)*v(627)+v(1130)*v
     &  (641)+v(679)*v(733)+(v(436)*v(548)+v(421)*v(666))*v(736)+v
     &  (737)*((-0.15d1)*v(524)+v(495)*v(741))
        v(745)=sigma(4)*v(1132)+v(1131)*v(510)+v(1129)*v(628)+v(1130
     &  )*v(642)+v(680)*v(733)+(v(436)*v(549)+v(421)*v(667))*v(736)+v
     &  (737)*((-0.15d1)*v(525)+v(1115)*v(741))
        v(746)=sigma(5)*v(1132)+v(1131)*v(511)+v(1129)*v(629)+v(1130
     &  )*v(643)+v(681)*v(733)+(v(436)*v(550)+v(421)*v(668))*v(736)+v
     &  (737)*((-0.15d1)*v(526)+v(1116)*v(741))
        v(747)=sigma(6)*v(1132)+v(1131)*v(512)+v(1129)*v(630)+v(1130
     &  )*v(644)+v(682)*v(733)+(v(436)*v(551)+v(421)*v(669))*v(736)+v
     &  (737)*((-0.15d1)*v(527)+v(1117)*v(741))
        v(748)=v(100)*v(1132)+v(1129)*v(631)+v(1130)*v(645)+v(683)*v
     &  (733)+(v(436)*v(552)+v(421)*v(670))*v(736)+v(737)*((-0.15d1
     &  )*v(528)+v(499)*v(741))
        v(749)=v(1129)*v(632)+v(1130)*v(646)+v(684)*v(733)+(v(326)*v
     &  (436)+v(421)*v(671))*v(736)+(-0.15d1)*v(529)*v(737)
        v(750)=v(1129)*v(633)+v(685)*v(733)+(v(327)*v(436)+v(421)*v
     &  (672))*v(736)
        v(751)=v(1129)*v(634)+v(686)*v(733)+(v(328)*v(436)+v(421)*v
     &  (673))*v(736)
        v(752)=v(1129)*v(635)+v(687)*v(733)+(v(329)*v(436)+v(421)*v
     &  (674))*v(736)
        v(753)=v(1129)*v(636)+v(688)*v(733)+(v(330)*v(436)+v(421)*v
     &  (675))*v(736)
        v(754)=v(1129)*v(637)+v(689)*v(733)+(v(331)*v(436)+v(421)*v
     &  (676))*v(736)
        v(755)=v(738)/2d0
        v(1157)=v(755)+v(836)
        v(756)=v(742)/2d0
        v(1158)=v(756)+v(836)
        v(757)=v(744)/2d0
        v(1159)=v(757)+v(836)
        v(758)=v(745)/2d0
        v(759)=v(746)/2d0
        v(760)=v(747)/2d0
        v(761)=v(748)/2d0
        v(762)=v(749)/2d0
        v(763)=v(750)/2d0
        v(764)=v(751)/2d0
        v(765)=v(752)/2d0
        v(766)=v(753)/2d0
        v(767)=v(754)/2d0
        v(769)=v(1105)*v(1134)+v(19)*v(738)+v(770)+v(1107)*(-v(474)+v
     &  (772))
        v(771)=v(1105)*v(1135)+v(19)*v(742)+v(770)+v(1107)*(-v(473)+v
     &  (773))
        v(774)=v(1105)*v(1136)+v(19)*v(744)+v(787)-v(1107)*v(789)
        v(775)=v(1105)*v(1137)+v(1112)*v(452)+v(19)*v(745)
        v(776)=v(1105)*v(1138)+v(19)*v(746)
        v(777)=v(1105)*v(1139)+v(19)*v(747)
        v(778)=v(1105)*v(1140)+v(19)*v(748)
        v(779)=v(1133)*v(684)+v(19)*v(749)
        v(780)=v(1133)*v(685)+v(19)*v(750)
        v(781)=v(1133)*v(686)+v(19)*v(751)
        v(782)=v(1133)*v(687)+v(19)*v(752)
        v(783)=v(1133)*v(688)+v(19)*v(753)
        v(784)=v(1133)*v(689)+v(19)*v(754)
        v(786)=v(18)*v(738)+v(770)+v(1107)*(-v(471)+v(772))+v(1134)*v
     &  (785)
        v(788)=v(18)*v(742)-v(1107)*v(773)+v(1135)*v(785)+v(787)
        v(790)=v(18)*v(744)+v(770)+v(1136)*v(785)+v(1107)*(-v(473)+v
     &  (789))
        v(791)=v(18)*v(745)+v(1137)*v(785)
        v(792)=v(18)*v(746)+v(1138)*v(785)
        v(794)=v(1114)*v(452)+v(18)*v(747)+v(1139)*v(785)
        v(795)=v(18)*v(748)+v(1140)*v(785)
        v(796)=v(1141)*v(684)+v(18)*v(749)
        v(797)=v(1141)*v(685)+v(18)*v(750)
        v(798)=v(1141)*v(686)+v(18)*v(751)
        v(799)=v(1141)*v(687)+v(18)*v(752)
        v(800)=v(1141)*v(688)+v(18)*v(753)
        v(801)=v(1141)*v(689)+v(18)*v(754)
        v(803)=v(17)*v(738)-v(1107)*v(772)+v(787)+v(1134)*v(802)
        v(804)=v(17)*v(742)+v(770)+v(1107)*(-v(471)+v(773))+v(1135)*v
     &  (802)
        v(805)=v(17)*v(744)+v(770)+v(1107)*(-v(474)+v(789))+v(1136)*v
     &  (802)
        v(806)=v(17)*v(745)+v(1137)*v(802)
        v(807)=v(1113)*v(452)+v(17)*v(746)+v(1138)*v(802)
        v(808)=v(17)*v(747)+v(1139)*v(802)
        v(809)=v(17)*v(748)+v(1140)*v(802)
        v(810)=v(1142)*v(684)+v(17)*v(749)
        v(811)=v(1142)*v(685)+v(17)*v(750)
        v(812)=v(1142)*v(686)+v(17)*v(751)
        v(813)=v(1142)*v(687)+v(17)*v(752)
        v(814)=v(1142)*v(688)+v(17)*v(753)
        v(815)=v(1142)*v(689)+v(17)*v(754)
        v(816)=(v(1143)*v(420)+v(436)*(v(420)*v(521)+v(28)*v(537))+v
     &  (31)*(-(v(625)/v(28))-v(638)/v(430))-v(769)-v(786)-v(803)+v
     &  (817))/3d0
        v(818)=(v(1144)*v(420)+v(436)*(v(420)*v(523)+v(28)*v(539))+v
     &  (31)*(-(v(626)/v(28))-v(640)/v(430))-v(771)-v(788)-v(804)+v
     &  (817))/3d0
        v(819)=(v(1145)*v(420)+v(436)*(v(420)*v(524)+v(28)*v(540))+v
     &  (31)*(-(v(627)/v(28))-v(641)/v(430))-v(774)-v(790)-v(805)+v
     &  (817))/3d0
        v(820)=(v(1146)*v(420)+v(436)*(v(420)*v(525)+v(28)*v(541))+v
     &  (31)*(-(v(628)/v(28))-v(642)/v(430))-v(775)-v(791)-v(806))
     &  /3d0
        v(821)=(v(1147)*v(420)+v(436)*(v(420)*v(526)+v(28)*v(542))+v
     &  (31)*(-(v(629)/v(28))-v(643)/v(430))-v(776)-v(792)-v(807))
     &  /3d0
        v(822)=(v(1148)*v(420)+v(436)*(v(420)*v(527)+v(28)*v(543))+v
     &  (31)*(-(v(630)/v(28))-v(644)/v(430))-v(777)-v(794)-v(808))
     &  /3d0
        v(823)=(v(1149)*v(420)+v(436)*(v(420)*v(528)+v(28)*v(544))+v
     &  (31)*(-(v(631)/v(28))-v(645)/v(430))-v(778)-v(795)-v(809)+v
     &  (824))/3d0
        v(825)=(-(v(1125)*v(420))+v(436)*(v(28)*v(316)+v(420)*v(529))
     &  +v(31)*(-(v(632)/v(28))-v(646)/v(430))-v(779)-v(796)-v(810)-v
     &  (824))/3d0
        v(826)=(v(1150)*v(633)+v(28)*(v(317)*v(436)+v(420)*v(672))-v
     &  (780)-v(797)-v(811))/3d0
        v(827)=(v(1150)*v(634)+v(28)*(v(318)*v(436)+v(420)*v(673))-v
     &  (781)-v(798)-v(812))/3d0
        v(828)=(v(1150)*v(635)+v(28)*(v(319)*v(436)+v(420)*v(674))-v
     &  (782)-v(799)-v(813))/3d0
        v(829)=(v(1150)*v(636)+v(28)*(v(320)*v(436)+v(420)*v(675))-v
     &  (783)-v(800)-v(814))/3d0
        v(830)=(v(1150)*v(637)+v(28)*(v(321)*v(436)+v(420)*v(676))-v
     &  (784)-v(801)-v(815))/3d0
        v(934)=v(1151)*v(423)+v(657)-v(28)*(v(348)*v(436)+v(423)*v
     &  (671))+v(100)*v(749)
        v(936)=-v(657)-v(28)*(v(358)*v(436)+v(424)*v(671))-v(935)
        v(937)=v(1152)-v(1109)*v(368)
        v(938)=v(1153)-v(1109)*v(378)
        v(939)=v(1154)-v(1109)*v(388)
        v(940)=v(1155)-v(1109)*v(399)
        v(941)=v(1156)-v(1109)*v(409)
        v(948)=v(658)-v(28)*(v(349)*v(436)+v(423)*v(672))+v(100)*v
     &  (750)
        v(950)=v(1152)-v(1109)*v(359)
        v(951)=-(v(28)*(v(369)*v(436)+v(425)*v(672)))
        v(952)=-(v(1109)*v(379))+v(965)
        v(953)=-(v(1109)*v(389))+v(980)
        v(954)=-(v(1109)*v(400))+v(996)
        v(955)=v(1013)-v(1109)*v(410)
        v(962)=v(659)-v(28)*(v(350)*v(436)+v(423)*v(673))+v(100)*v
     &  (751)
        v(964)=v(1153)-v(1109)*v(360)
        v(966)=-(v(1109)*v(370))+v(965)
        v(967)=-(v(28)*(v(380)*v(436)+v(426)*v(673)))
        v(968)=-(v(1109)*v(390))+v(982)
        v(969)=-(v(1109)*v(401))+v(998)
        v(970)=v(1015)-v(1109)*v(411)
        v(977)=v(660)-v(28)*(v(351)*v(436)+v(423)*v(674))+v(100)*v
     &  (752)
        v(979)=v(1154)-v(1109)*v(361)
        v(981)=-(v(1109)*v(371))+v(980)
        v(983)=-(v(1109)*v(381))+v(982)
        v(984)=-(v(28)*(v(391)*v(436)+v(427)*v(674)))
        v(985)=v(1000)-v(1109)*v(402)
        v(986)=v(1017)-v(1109)*v(412)
        v(993)=v(661)-v(28)*(v(352)*v(436)+v(423)*v(675))+v(100)*v
     &  (753)
        v(995)=v(1155)-v(1109)*v(362)
        v(997)=-(v(1109)*v(372))+v(996)
        v(999)=-(v(1109)*v(382))+v(998)
        v(1001)=v(1000)-v(1109)*v(393)
        v(1002)=-(v(28)*(v(403)*v(436)+v(428)*v(675)))
        v(1003)=v(1019)-v(1109)*v(413)
        v(1010)=v(662)-v(28)*(v(353)*v(436)+v(423)*v(676))+v(100)*v
     &  (754)
        v(1012)=v(1156)-v(1109)*v(363)
        v(1014)=v(1013)-v(1109)*v(373)
        v(1016)=v(1015)-v(1109)*v(383)
        v(1018)=v(1017)-v(1109)*v(394)
        v(1020)=v(1019)-v(1109)*v(404)
        v(1021)=-(v(28)*(v(414)*v(436)+v(429)*v(676)))
        DDF(1,1)=v(803)+v(816)
        DDF(1,2)=v(804)+v(818)
        DDF(1,3)=v(805)+v(819)
        DDF(1,4)=v(806)+v(820)
        DDF(1,5)=v(807)+v(821)
        DDF(1,6)=v(808)+v(822)
        DDF(1,7)=v(809)+v(823)
        DDF(1,8)=v(810)+v(825)
        DDF(1,9)=v(811)+v(826)
        DDF(1,10)=v(812)+v(827)
        DDF(1,11)=v(813)+v(828)
        DDF(1,12)=v(814)+v(829)
        DDF(1,13)=v(815)+v(830)
        DDF(2,1)=v(786)+v(816)
        DDF(2,2)=v(788)+v(818)
        DDF(2,3)=v(790)+v(819)
        DDF(2,4)=v(791)+v(820)
        DDF(2,5)=v(792)+v(821)
        DDF(2,6)=v(794)+v(822)
        DDF(2,7)=v(795)+v(823)
        DDF(2,8)=v(796)+v(825)
        DDF(2,9)=v(797)+v(826)
        DDF(2,10)=v(798)+v(827)
        DDF(2,11)=v(799)+v(828)
        DDF(2,12)=v(800)+v(829)
        DDF(2,13)=v(801)+v(830)
        DDF(3,1)=v(769)+v(816)
        DDF(3,2)=v(771)+v(818)
        DDF(3,3)=v(774)+v(819)
        DDF(3,4)=v(775)+v(820)
        DDF(3,5)=v(776)+v(821)
        DDF(3,6)=v(777)+v(822)
        DDF(3,7)=v(778)+v(823)
        DDF(3,8)=v(779)+v(825)
        DDF(3,9)=v(780)+v(826)
        DDF(3,10)=v(781)+v(827)
        DDF(3,11)=v(782)+v(828)
        DDF(3,12)=v(783)+v(829)
        DDF(3,13)=v(784)+v(830)
        DDF(4,1)=v(1112)*(v(1157)+v(19)*v(706))+sigma(5)*v(720)
        DDF(4,2)=v(1112)*(v(1158)+v(19)*v(708))+sigma(5)*v(721)
        DDF(4,3)=sigma(5)*v(722)+v(1112)*(v(19)*v(709)+v(757)+v(851))
        DDF(4,4)=sigma(5)*v(723)+v(1112)*(v(19)*v(710)+v(758))+2d0*v
     &  (876)
        DDF(4,5)=v(451)+sigma(5)*v(724)+v(1112)*(v(19)*v(711)+v(759))
        DDF(4,6)=sigma(5)*v(725)+v(1112)*(v(19)*v(712)+v(760))
        DDF(4,7)=sigma(5)*v(726)+v(1112)*(v(19)*v(713)+v(761))
        DDF(4,8)=sigma(5)*v(727)+v(1112)*(v(19)*v(714)+v(762))
        DDF(4,9)=sigma(5)*v(728)+v(1112)*(v(19)*v(715)+v(763))
        DDF(4,10)=sigma(5)*v(729)+v(1112)*(v(19)*v(716)+v(764))
        DDF(4,11)=sigma(5)*v(730)+v(1112)*(v(19)*v(717)+v(765))
        DDF(4,12)=sigma(5)*v(731)+v(1112)*(v(19)*v(718)+v(766))
        DDF(4,13)=sigma(5)*v(732)+v(1112)*(v(19)*v(719)+v(767))
        DDF(5,1)=sigma(4)*v(720)+v(1113)*(v(17)*v(706)+v(755)+v(851))
        DDF(5,2)=v(1113)*(v(1158)+v(17)*v(708))+sigma(4)*v(721)
        DDF(5,3)=v(1113)*(v(1159)+v(17)*v(709))+sigma(4)*v(722)
        DDF(5,4)=v(451)+sigma(4)*v(723)+v(1113)*(v(17)*v(710)+v(758))
        DDF(5,5)=sigma(4)*v(724)+v(1113)*(v(17)*v(711)+v(759))+2d0*v
     &  (891)
        DDF(5,6)=sigma(4)*v(725)+v(1113)*(v(17)*v(712)+v(760))
        DDF(5,7)=sigma(4)*v(726)+v(1113)*(v(17)*v(713)+v(761))
        DDF(5,8)=sigma(4)*v(727)+v(1113)*(v(17)*v(714)+v(762))
        DDF(5,9)=sigma(4)*v(728)+v(1113)*(v(17)*v(715)+v(763))
        DDF(5,10)=sigma(4)*v(729)+v(1113)*(v(17)*v(716)+v(764))
        DDF(5,11)=sigma(4)*v(730)+v(1113)*(v(17)*v(717)+v(765))
        DDF(5,12)=sigma(4)*v(731)+v(1113)*(v(17)*v(718)+v(766))
        DDF(5,13)=sigma(4)*v(732)+v(1113)*(v(17)*v(719)+v(767))
        DDF(6,1)=v(1122)*v(691)+v(1114)*(v(1157)+v(18)*v(706))
        DDF(6,2)=v(1122)*v(693)+v(1114)*(v(18)*v(708)+v(756)+v(851))
        DDF(6,3)=v(1122)*v(694)+v(1114)*(v(1159)+v(18)*v(709))
        DDF(6,4)=sigma(5)*(v(455)+sigma(4)*v(695))+v(1114)*(v(18)*v
     &  (710)+v(758))
        DDF(6,5)=sigma(4)*(v(455)+sigma(5)*v(696))+v(1114)*(v(18)*v
     &  (711)+v(759))
        DDF(6,6)=v(1122)*v(697)+v(1114)*(v(18)*v(712)+v(760))+2d0*v
     &  (906)
        DDF(6,7)=v(1122)*v(698)+v(1114)*(v(18)*v(713)+v(761))
        DDF(6,8)=v(1122)*v(699)+v(1114)*(v(18)*v(714)+v(762))
        DDF(6,9)=v(1122)*v(700)+v(1114)*(v(18)*v(715)+v(763))
        DDF(6,10)=v(1122)*v(701)+v(1114)*(v(18)*v(716)+v(764))
        DDF(6,11)=v(1122)*v(702)+v(1114)*(v(18)*v(717)+v(765))
        DDF(6,12)=v(1122)*v(703)+v(1114)*(v(18)*v(718)+v(766))
        DDF(6,13)=v(1122)*v(704)+v(1114)*(v(18)*v(719)+v(767))
        DDF(7,1)=-(v(1143)*v(423))+v(436)*(-(v(423)*v(521))-v(28)*v
     &  (561))+v(647)+v(100)*v(738)
        DDF(7,2)=-(v(1144)*v(423))+v(436)*(-(v(423)*v(523))-v(28)*v
     &  (563))+v(650)+v(100)*v(742)
        DDF(7,3)=-(v(1145)*v(423))+v(436)*(-(v(423)*v(524))-v(28)*v
     &  (564))+v(651)+v(100)*v(744)
        DDF(7,4)=-(v(1146)*v(423))+v(436)*(-(v(423)*v(525))-v(28)*v
     &  (565))+v(652)+v(100)*v(745)
        DDF(7,5)=-(v(1147)*v(423))+v(436)*(-(v(423)*v(526))-v(28)*v
     &  (566))+v(653)+v(100)*v(746)
        DDF(7,6)=-(v(1148)*v(423))+v(436)*(-(v(423)*v(527))-v(28)*v
     &  (567))+v(654)+v(100)*v(747)
        DDF(7,7)=-(v(1149)*v(423))+2d0*v(1037)*v(438)+v(436)*(-(v(423
     &  )*v(528))-v(28)*v(568))+v(655)+v(100)*v(748)
        DDF(7,8)=v(934)
        DDF(7,9)=v(948)
        DDF(7,10)=v(962)
        DDF(7,11)=v(977)
        DDF(7,12)=v(993)
        DDF(7,13)=v(1010)
        DDF(8,1)=-(v(1143)*v(424))+v(436)*(-(v(424)*v(521))-v(28)*v
     &  (569))-v(647)
        DDF(8,2)=-(v(1144)*v(424))+v(436)*(-(v(424)*v(523))-v(28)*v
     &  (571))-v(650)
        DDF(8,3)=-(v(1145)*v(424))+v(436)*(-(v(424)*v(524))-v(28)*v
     &  (572))-v(651)
        DDF(8,4)=-(v(1146)*v(424))+v(436)*(-(v(424)*v(525))-v(28)*v
     &  (573))-v(652)
        DDF(8,5)=-(v(1147)*v(424))+v(436)*(-(v(424)*v(526))-v(28)*v
     &  (574))-v(653)
        DDF(8,6)=-(v(1148)*v(424))+v(436)*(-(v(424)*v(527))-v(28)*v
     &  (575))-v(654)
        DDF(8,7)=-(v(1149)*v(424))+v(436)*(-(v(424)*v(528))-v(28)*v
     &  (576))-v(655)
        DDF(8,8)=v(936)
        DDF(8,9)=v(950)
        DDF(8,10)=v(964)
        DDF(8,11)=v(979)
        DDF(8,12)=v(995)
        DDF(8,13)=v(1012)
        DDF(9,1)=-(v(1143)*v(425))+v(436)*(-(v(425)*v(521))-v(28)*v
     &  (577))
        DDF(9,2)=-(v(1144)*v(425))+v(436)*(-(v(425)*v(523))-v(28)*v
     &  (579))
        DDF(9,3)=-(v(1145)*v(425))+v(436)*(-(v(425)*v(524))-v(28)*v
     &  (580))
        DDF(9,4)=-(v(1146)*v(425))+v(436)*(-(v(425)*v(525))-v(28)*v
     &  (581))
        DDF(9,5)=-(v(1147)*v(425))+v(436)*(-(v(425)*v(526))-v(28)*v
     &  (582))
        DDF(9,6)=-(v(1148)*v(425))+v(436)*(-(v(425)*v(527))-v(28)*v
     &  (583))
        DDF(9,7)=-(v(1149)*v(425))+v(436)*(-(v(425)*v(528))-v(28)*v
     &  (584))
        DDF(9,8)=v(937)
        DDF(9,9)=v(951)
        DDF(9,10)=v(966)
        DDF(9,11)=v(981)
        DDF(9,12)=v(997)
        DDF(9,13)=v(1014)
        DDF(10,1)=-(v(1143)*v(426))+v(436)*(-(v(426)*v(521))-v(28)*v
     &  (585))
        DDF(10,2)=-(v(1144)*v(426))+v(436)*(-(v(426)*v(523))-v(28)*v
     &  (587))
        DDF(10,3)=-(v(1145)*v(426))+v(436)*(-(v(426)*v(524))-v(28)*v
     &  (588))
        DDF(10,4)=-(v(1146)*v(426))+v(436)*(-(v(426)*v(525))-v(28)*v
     &  (589))
        DDF(10,5)=-(v(1147)*v(426))+v(436)*(-(v(426)*v(526))-v(28)*v
     &  (590))
        DDF(10,6)=-(v(1148)*v(426))+v(436)*(-(v(426)*v(527))-v(28)*v
     &  (591))
        DDF(10,7)=-(v(1149)*v(426))+v(436)*(-(v(426)*v(528))-v(28)*v
     &  (592))
        DDF(10,8)=v(938)
        DDF(10,9)=v(952)
        DDF(10,10)=v(967)
        DDF(10,11)=v(983)
        DDF(10,12)=v(999)
        DDF(10,13)=v(1016)
        DDF(11,1)=-(v(1143)*v(427))+v(436)*(-(v(427)*v(521))-v(28)*v
     &  (593))
        DDF(11,2)=-(v(1144)*v(427))+v(436)*(-(v(427)*v(523))-v(28)*v
     &  (595))
        DDF(11,3)=-(v(1145)*v(427))+v(436)*(-(v(427)*v(524))-v(28)*v
     &  (596))
        DDF(11,4)=-(v(1146)*v(427))+v(436)*(-(v(427)*v(525))-v(28)*v
     &  (597))
        DDF(11,5)=-(v(1147)*v(427))+v(436)*(-(v(427)*v(526))-v(28)*v
     &  (598))
        DDF(11,6)=-(v(1148)*v(427))+v(436)*(-(v(427)*v(527))-v(28)*v
     &  (599))
        DDF(11,7)=-(v(1149)*v(427))+v(436)*(-(v(427)*v(528))-v(28)*v
     &  (600))
        DDF(11,8)=v(939)
        DDF(11,9)=v(953)
        DDF(11,10)=v(968)
        DDF(11,11)=v(984)
        DDF(11,12)=v(1001)
        DDF(11,13)=v(1018)
        DDF(12,1)=-(v(1143)*v(428))+v(436)*(-(v(428)*v(521))-v(28)*v
     &  (601))
        DDF(12,2)=-(v(1144)*v(428))+v(436)*(-(v(428)*v(523))-v(28)*v
     &  (603))
        DDF(12,3)=-(v(1145)*v(428))+v(436)*(-(v(428)*v(524))-v(28)*v
     &  (604))
        DDF(12,4)=-(v(1146)*v(428))+v(436)*(-(v(428)*v(525))-v(28)*v
     &  (605))
        DDF(12,5)=-(v(1147)*v(428))+v(436)*(-(v(428)*v(526))-v(28)*v
     &  (606))
        DDF(12,6)=-(v(1148)*v(428))+v(436)*(-(v(428)*v(527))-v(28)*v
     &  (607))
        DDF(12,7)=-(v(1149)*v(428))+v(436)*(-(v(428)*v(528))-v(28)*v
     &  (608))
        DDF(12,8)=v(940)
        DDF(12,9)=v(954)
        DDF(12,10)=v(969)
        DDF(12,11)=v(985)
        DDF(12,12)=v(1002)
        DDF(12,13)=v(1020)
        DDF(13,1)=v(436)*(-(v(429)*v(521))-v(28)*v(609))+v(1160)*v
     &  (663)
        DDF(13,2)=-(v(1144)*v(429))+v(436)*(-(v(429)*v(523))-v(28)*v
     &  (611))
        DDF(13,3)=-(v(1145)*v(429))+v(436)*(-(v(429)*v(524))-v(28)*v
     &  (612))
        DDF(13,4)=-(v(1146)*v(429))+v(436)*(-(v(429)*v(525))-v(28)*v
     &  (613))
        DDF(13,5)=-(v(1147)*v(429))+v(436)*(-(v(429)*v(526))-v(28)*v
     &  (614))
        DDF(13,6)=-(v(1148)*v(429))+v(436)*(-(v(429)*v(527))-v(28)*v
     &  (615))
        DDF(13,7)=-(v(1149)*v(429))+v(436)*(-(v(429)*v(528))-v(28)*v
     &  (616))
        DDF(13,8)=v(941)
        DDF(13,9)=v(955)
        DDF(13,10)=v(970)
        DDF(13,11)=v(986)
        DDF(13,12)=v(1003)
        DDF(13,13)=v(1021)
       ELSE
        DDF(1,1)=0d0
        DDF(1,2)=0d0
        DDF(1,3)=0d0
        DDF(1,4)=0d0
        DDF(1,5)=0d0
        DDF(1,6)=0d0
        DDF(1,7)=0d0
        DDF(1,8)=0d0
        DDF(1,9)=0d0
        DDF(1,10)=0d0
        DDF(1,11)=0d0
        DDF(1,12)=0d0
        DDF(1,13)=0d0
        DDF(2,2)=0d0
        DDF(2,3)=0d0
        DDF(2,4)=0d0
        DDF(2,5)=0d0
        DDF(2,6)=0d0
        DDF(2,7)=0d0
        DDF(2,8)=0d0
        DDF(2,9)=0d0
        DDF(2,10)=0d0
        DDF(2,11)=0d0
        DDF(2,12)=0d0
        DDF(2,13)=0d0
        DDF(3,3)=0d0
        DDF(3,4)=0d0
        DDF(3,5)=0d0
        DDF(3,6)=0d0
        DDF(3,7)=0d0
        DDF(3,8)=0d0
        DDF(3,9)=0d0
        DDF(3,10)=0d0
        DDF(3,11)=0d0
        DDF(3,12)=0d0
        DDF(3,13)=0d0
        DDF(4,4)=0d0
        DDF(4,5)=0d0
        DDF(4,6)=0d0
        DDF(4,7)=0d0
        DDF(4,8)=0d0
        DDF(4,9)=0d0
        DDF(4,10)=0d0
        DDF(4,11)=0d0
        DDF(4,12)=0d0
        DDF(4,13)=0d0
        DDF(5,5)=0d0
        DDF(5,6)=0d0
        DDF(5,7)=0d0
        DDF(5,8)=0d0
        DDF(5,9)=0d0
        DDF(5,10)=0d0
        DDF(5,11)=0d0
        DDF(5,12)=0d0
        DDF(5,13)=0d0
        DDF(6,6)=0d0
        DDF(6,7)=0d0
        DDF(6,8)=0d0
        DDF(6,9)=0d0
        DDF(6,10)=0d0
        DDF(6,11)=0d0
        DDF(6,12)=0d0
        DDF(6,13)=0d0
        DDF(7,7)=0d0
        DDF(7,8)=0d0
        DDF(7,9)=0d0
        DDF(7,10)=0d0
        DDF(7,11)=0d0
        DDF(7,12)=0d0
        DDF(7,13)=0d0
        DDF(8,8)=0d0
        DDF(8,9)=0d0
        DDF(8,10)=0d0
        DDF(8,11)=0d0
        DDF(8,12)=0d0
        DDF(8,13)=0d0
        DDF(9,9)=0d0
        DDF(9,10)=0d0
        DDF(9,11)=0d0
        DDF(9,12)=0d0
        DDF(9,13)=0d0
        DDF(10,10)=0d0
        DDF(10,11)=0d0
        DDF(10,12)=0d0
        DDF(10,13)=0d0
        DDF(11,11)=0d0
        DDF(11,12)=0d0
        DDF(11,13)=0d0
        DDF(12,12)=0d0
        DDF(12,13)=0d0
        DDF(13,13)=0d0
        DO i01=2,13
          DO i02=1,i01-1
            DDF(i01,i02)=DDF(i02,i01)
          ENDDO
        ENDDO
       ENDIF
      ELSE
       DF(1)=0d0
       DF(2)=0d0
       DF(3)=0d0
       DF(4)=0d0
       DF(5)=0d0
       DF(6)=0d0
       DF(7)=0d0
       DF(8)=0d0
       DF(9)=0d0
       DF(10)=0d0
       DF(11)=0d0
       DF(12)=0d0
       DF(13)=0d0
      ENDIF
      END
