%**************************************************************
%* AceGen    5.101 Windows (28 May 13)                        *
%*           Co. J. Korelc  2007            19 Mar 14 10:42:26*
%**************************************************************
% User     : Full professional version
% Notebook : BP_yield_function.nb
% Evaluation time                 : 41 s    Mode  : Optimal
% Number of formulae              : 873     Method: Automatic
% Subroutine                      : BPYieldFunction size :15006
% Total size of Mathematica  code : 15006 subexpressions
% Total size of Matlab code      : 41628 bytes

%*********************** F U N C T I O N **************************
function[F,DF,DDF,iter]=BPYieldFunction(sigma,data,dorder);
persistent v;
if size(v)<5005
  v=zeros(5005,'double');
end;
v(1121)=0.7071067811865476e0*sigma(6);
v(1039)=-0.5235987755982988e0*data(5);
v(1038)=(-sigma(1)-sigma(2)-sigma(3))/3e0;
v(1037)=Power(data(8),2);
v(1035)=data(7);
v(1036)=2e0*v(1035);
v(1034)=data(6);
v(1086)=Power(v(1034),3);
v(1033)=data(4);
v(1032)=data(3);
v(1031)=data(2);
v(1030)=data(1);
v(1029)=Power(sigma(6),2);
v(1104)=-v(1029)/2e0;
v(1028)=sigma(5)*v(1121);
v(1027)=Power(sigma(5),2);
v(1103)=-v(1027)/2e0;
v(1026)=Power(sigma(4),2);
v(96)=2e0*v(1030);
v(39)=(v(1030)*v(1030));
v(1077)=v(1035)*v(39);
v(37)=v(1030)+v(1031);
v(279)=-2e0+v(1032);
v(50)=-1e0+v(1032);
v(45)=1e0-v(1033);
v(1040)=2e0*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);
v(17)=sigma(1)+v(1038);
v(18)=sigma(2)+v(1038);
v(1105)=-v(1026)/2e0+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)))/2e0+v(1037)*v(39);
v(1106)=sqrt(v(20));
v(449)=1/Power(v(1106),3);
v(439)=0.25980762113533157e1*v(449);
v(437)=3e0*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)/2e0;
v(27)=-v(1031)+v(26);
v(31)=v(1038)-v(27);
v(133)=2e0*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)=Power(v(29),3);
v(1072)=v(1034)*v(29);
v(242)=(v(29)*v(29));
v(1091)=-(v(234)*v(242));
v(232)=1e0+v(1091);
v(168)=1/Sqrt(v(232));
v(1090)=(2e0/3e0)*v(168);
v(82)=v(1039)+acos(v(1072))/3e0;
v(178)=tan(v(82));
v(30)=sec(v(82));
v(236)=1/Power(v(30),2);
v(1044)=2e0*v(236);
v(32)=atan2(v(25),v(31));
v(43)=sin(v(32));
v(1045)=v(1044)*(v(43)*v(43));
v(41)=cos(v(32));
v(42)=v(41)/v(37);
v(57)=v(1040)*v(42);
v(33)=v(26);
for 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.20943951023931953e1*v(36);
 v(295)=(4e0/3e0)*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)=2e0*v(38);
 v(151)=1e0-v(1042);
 v(152)=v(1041)*v(151);
 v(149)=v(1042)*v(45);
 v(53)=sign(v(38));
 v(1063)=v(38)*v(53);
 v(67)=-(v(1063)*v(50));
 v(54)=v(42)*v(53);
 v(49)=abs(v(38));
 v(278)=Power(v(49),-3e0+v(1032));
 v(147)=Power(v(49),v(279));
 v(1048)=v(147)*v(54);
 v(146)=log(v(49));
 v(136)=v(278)*v(279)*v(67);
 v(56)=v(1048)*v(50);
 v(48)=Power(v(49),v(50));
 v(1043)=-1e0+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)-2e0*v(56))+2e0*v(57)*v(58);
 v(64)=-(v(40)/v(52));
 v(33)=v(33)+v(64);
 if(abs(v(64))<0.1e-9 && abs(v(40))<0.1e-7)
  v(1097)=v(38)*v(66);
  v(1079)=-(v(43)/v(37));
  v(1078)=-2e0*(v(281)*v(47)+v(57)*v(85));
  v(1075)=-2e0*v(44);
  v(1074)=v(31)/v(185);
  v(1073)=2e0*v(47);
  v(1068)=v(41)*v(69);
  v(1064)=1/Power(v(37),3);
  v(1062)=v(146)*v(48);
  v(1061)=v(147)*v(50);
  v(1067)=(-2e0*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)=1/Power(v(185),2);
  v(1052)=1/Power(v(232),0.15e1);
  v(1050)=sin(v(82));
  v(1047)=1/Power(v(37),2);
  v(1046)=cos(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)=2e0*v(1068)-v(43)*(v(287)+v(285)*v(71));
  v(293)=v(1058)*v(288);
  v(1080)=2e0*v(293);
  v(290)=v(285)+v(293);
  v(292)=(v(290)-2e0*v(1047)*(v(275)*v(41)+v(284)*v(76)))/2e0;
  v(268)=0.3490658503988659e0*v(168)*v(83);
  v(253)=-(v(1071)*v(85));
  v(249)=4e0*v(253)*v(45);
  v(243)=((2e0/9e0)*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)=(2e0/3e0)*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)=2e0*v(138);
  v(251)=v(1062)*v(38);
  v(121)=-(v(1047)*v(76));
  v(104)=v(1058)/2e0;
  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.5e0+v(1060);
  v(114)=0.5e0-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)=(-2e0*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)/2e0;
  v(205)=v(1058)*v(183);
  v(217)=v(205)/2e0;
  v(204)=v(1058)*v(180);
  v(216)=v(204)/2e0;
  v(202)=v(1058)*v(177);
  v(203)=v(163)+v(202);
  v(215)=v(121)*v(153)+v(203)/2e0;
  v(200)=v(1058)*v(176);
  v(201)=v(162)+v(200);
  v(214)=v(121)*v(150)+v(201)/2e0;
  v(194)=v(1058)*v(170);
  v(210)=v(194)/2e0;
  v(212)=2e0*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)/(2e0*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)/2e0+v(212)-v(1047)*(v(114)*v(275)+v(144)*v(76));
  v(195)=v(1058)*v(171)+v(197);
  v(78)=-2e0*v(223)*v(25);
  v(75)=v(287)+v(78);
  v(77)=v(75)/2e0+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.10471975511965977e1*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.5483113556160755e0*v(83);
  v(294)=v(1074)*v(288);
  v(297)=v(1034)*v(1076);
  v(1089)=2e0*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)=-2e0*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)+(-2e0*v(1053)*v(1081)+1/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)+2e0*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)=-(2e0*v(121)*v(142)+v(159)+2e0*v(112)*v(160)+v(195)+2e0*v(212)+2e0*v(120)*v(249)+2e0*v(247 ...
   )*v(85)+4e0*v(300)*v(86)+2e0*v(63)*(v(86)*v(86))+v(246)*v(247)*v(96))/(2e0*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)+2e0*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)-4e0*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)+2e0*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)-2e0*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)+2e0*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)-2e0*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)-2e0*v(310)*v(95)-v(63)*(v(95)*v(95)))/v(52);
  iter=i34;
  break;
 else;
 end;
 if(i34==15)
  iter=-1e0;
  break;
 else;
 end;
end;
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=-1e0+v(28)/v(430);
if(dorder>=1)
 v(1123)=0.15e1/v(25);
 v(1122)=sigma(4)*sigma(5);
 v(1114)=2e0*sigma(6);
 v(1113)=2e0*sigma(5);
 v(1112)=2e0*sigma(4);
 v(1111)=v(17)*v(19);
 v(1110)=1/(2e0*v(28));
 v(1108)=1/(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)=1/Power(v(20),2);
 v(734)=(-0.3897114317029974e1*v(477))/v(1106);
 v(733)=v(24)*v(734);
 v(529)=v(1110)*v(31);
 v(458)=-(v(529)/v(430));
 v(436)=1/Power(v(430),2);
 v(1109)=v(28)*v(436);
 v(434)=-(v(1109)*v(422));
 v(1107)=v(434)*v(439);
 v(705)=-0.12990381056766578e1*v(434);
 v(690)=0.18371173070873836e1*v(434);
 v(455)=v(449)*v(690);
 v(452)=v(449)*v(705);
 v(451)=sigma(6)*v(455);
 v(438)=0.15e1*v(1108)+v(434)*v(733)+v(436)*v(735);
 v(453)=v(438)/2e0;
 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))/3e0;
 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(dorder==2)
  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.15e1/v(430);
  v(1129)=0.15e1/v(28);
  v(1117)=3e0*sigma(6);
  v(1116)=3e0*sigma(5);
  v(1115)=3e0*sigma(4);
  v(851)=(2e0/3e0)*v(452);
  v(836)=-v(452)/3e0;
  v(824)=v(1110)/v(430);
  v(817)=v(1108)/3e0;
  v(789)=-v(19)/3e0;
  v(787)=(2e0/3e0)*v(438);
  v(773)=-v(18)/3e0;
  v(772)=-v(17)/3e0;
  v(770)=-v(438)/3e0;
  v(741)=(0.75e0*v(28))/v(437);
  v(740)=(0.5625e1*v(434))/Power(v(20),3);
  v(739)=(0.17320508075688773e1*v(24))/v(1106);
  v(1132)=v(739)*v(740);
  v(707)=-0.12990381056766578e1*v(449);
  v(692)=0.18371173070873836e1*v(449);
  v(664)=1/Power(v(430),3);
  v(1124)=-2e0*v(664);
  v(656)=v(1108)/4e0;
  v(649)=-v(31)/(2e0*v(430));
  v(648)=v(1108)/6e0;
  v(639)=1/Power(v(28),2);
  v(618)=-v(420)/3e0;
  v(610)=-v(405)/3e0;
  v(602)=-v(395)/3e0;
  v(594)=-v(384)/3e0;
  v(586)=-v(374)/3e0;
  v(578)=-v(364)/3e0;
  v(570)=-v(354)/3e0;
  v(562)=-v(344)/3e0;
  v(554)=-v(332)/3e0;
  v(546)=-v(322)/3e0;
  v(538)=-v(311)/3e0;
  v(522)=(-2e0/3e0)*v(31);
  v(507)=v(1029)/6e0;
  v(506)=v(1027)/6e0;
  v(505)=v(19)*v(773);
  v(1119)=v(505)+v(506);
  v(504)=-v(1111)/3e0;
  v(1120)=v(504)+v(507);
  v(502)=v(1026)/6e0;
  v(501)=v(18)*v(772);
  v(1118)=v(501)+v(502);
  v(479)=0.3849001794597505e0*v(734);
  v(474)=(-2e0/3e0)*v(18);
  v(473)=(-2e0/3e0)*v(17);
  v(471)=(-2e0/3e0)*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.25980762113533157e1*v(478);
  v(487)=0.25980762113533157e1*v(480);
  v(488)=0.25980762113533157e1*v(481);
  v(489)=0.25980762113533157e1*v(482);
  v(490)=0.25980762113533157e1*v(483);
  v(491)=0.25980762113533157e1*v(484);
  v(492)=0.25980762113533157e1*v(485);
  v(493)=3e0*v(17);
  v(494)=3e0*v(18);
  v(495)=3e0*v(19);
  v(499)=3e0*v(100);
  v(500)=-v(1027)/3e0+v(1118)+v(1120)-v(18)*v(471);
  v(503)=-v(1029)/3e0+v(1118)+v(1119)-v(17)*v(471);
  v(508)=-v(1026)/3e0+v(1119)+v(1120)-v(18)*v(473);
  v(510)=v(1028)-(v(1112)*v(19))/2e0;
  v(511)=sigma(4)*v(1121)-(v(1113)*v(17))/2e0;
  v(512)=0.7071067811865476e0*v(1122)-(v(1114)*v(18))/2e0;
  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.15e1*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.15e1*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.15e1*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.15e1*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.15e1*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.15e1*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.15e1*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.15e1*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)/2e0;
  v(1157)=v(755)+v(836);
  v(756)=v(742)/2e0;
  v(1158)=v(756)+v(836);
  v(757)=v(744)/2e0;
  v(1159)=v(757)+v(836);
  v(758)=v(745)/2e0;
  v(759)=v(746)/2e0;
  v(760)=v(747)/2e0;
  v(761)=v(748)/2e0;
  v(762)=v(749)/2e0;
  v(763)=v(750)/2e0;
  v(764)=v(751)/2e0;
  v(765)=v(752)/2e0;
  v(766)=v(753)/2e0;
  v(767)=v(754)/2e0;
  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))/3e0;
  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))/3e0;
  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))/3e0;
  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))/3e0;
  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))/3e0;
  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))/3e0;
  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))/3e0;
  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))/3e0;
  v(826)=(v(1150)*v(633)+v(28)*(v(317)*v(436)+v(420)*v(672))-v(780)-v(797)-v(811))/3e0;
  v(827)=(v(1150)*v(634)+v(28)*(v(318)*v(436)+v(420)*v(673))-v(781)-v(798)-v(812))/3e0;
  v(828)=(v(1150)*v(635)+v(28)*(v(319)*v(436)+v(420)*v(674))-v(782)-v(799)-v(813))/3e0;
  v(829)=(v(1150)*v(636)+v(28)*(v(320)*v(436)+v(420)*v(675))-v(783)-v(800)-v(814))/3e0;
  v(830)=(v(1150)*v(637)+v(28)*(v(321)*v(436)+v(420)*v(676))-v(784)-v(801)-v(815))/3e0;
  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))+2e0*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))+2e0*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))+2e0*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))+2e0*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)=0;
  DDF(1,2)=0;
  DDF(1,3)=0;
  DDF(1,4)=0;
  DDF(1,5)=0;
  DDF(1,6)=0;
  DDF(1,7)=0;
  DDF(1,8)=0;
  DDF(1,9)=0;
  DDF(1,10)=0;
  DDF(1,11)=0;
  DDF(1,12)=0;
  DDF(1,13)=0;
  DDF(2,2)=0;
  DDF(2,3)=0;
  DDF(2,4)=0;
  DDF(2,5)=0;
  DDF(2,6)=0;
  DDF(2,7)=0;
  DDF(2,8)=0;
  DDF(2,9)=0;
  DDF(2,10)=0;
  DDF(2,11)=0;
  DDF(2,12)=0;
  DDF(2,13)=0;
  DDF(3,3)=0;
  DDF(3,4)=0;
  DDF(3,5)=0;
  DDF(3,6)=0;
  DDF(3,7)=0;
  DDF(3,8)=0;
  DDF(3,9)=0;
  DDF(3,10)=0;
  DDF(3,11)=0;
  DDF(3,12)=0;
  DDF(3,13)=0;
  DDF(4,4)=0;
  DDF(4,5)=0;
  DDF(4,6)=0;
  DDF(4,7)=0;
  DDF(4,8)=0;
  DDF(4,9)=0;
  DDF(4,10)=0;
  DDF(4,11)=0;
  DDF(4,12)=0;
  DDF(4,13)=0;
  DDF(5,5)=0;
  DDF(5,6)=0;
  DDF(5,7)=0;
  DDF(5,8)=0;
  DDF(5,9)=0;
  DDF(5,10)=0;
  DDF(5,11)=0;
  DDF(5,12)=0;
  DDF(5,13)=0;
  DDF(6,6)=0;
  DDF(6,7)=0;
  DDF(6,8)=0;
  DDF(6,9)=0;
  DDF(6,10)=0;
  DDF(6,11)=0;
  DDF(6,12)=0;
  DDF(6,13)=0;
  DDF(7,7)=0;
  DDF(7,8)=0;
  DDF(7,9)=0;
  DDF(7,10)=0;
  DDF(7,11)=0;
  DDF(7,12)=0;
  DDF(7,13)=0;
  DDF(8,8)=0;
  DDF(8,9)=0;
  DDF(8,10)=0;
  DDF(8,11)=0;
  DDF(8,12)=0;
  DDF(8,13)=0;
  DDF(9,9)=0;
  DDF(9,10)=0;
  DDF(9,11)=0;
  DDF(9,12)=0;
  DDF(9,13)=0;
  DDF(10,10)=0;
  DDF(10,11)=0;
  DDF(10,12)=0;
  DDF(10,13)=0;
  DDF(11,11)=0;
  DDF(11,12)=0;
  DDF(11,13)=0;
  DDF(12,12)=0;
  DDF(12,13)=0;
  DDF(13,13)=0;
  for i01=1:13:1; for i02=1:i01:1;
      DDF(i01,i02)=DDF(i02,i01);
    end;
  end;
 end;
else;
 DF(1)=0;
 DF(2)=0;
 DF(3)=0;
 DF(4)=0;
 DF(5)=0;
 DF(6)=0;
 DF(7)=0;
 DF(8)=0;
 DF(9)=0;
 DF(10)=0;
 DF(11)=0;
 DF(12)=0;
 DF(13)=0;
end;


function [x]=SMSKDelta(i,j)
if (i==j) , x=1; else x=0; end;
end
function [x]=SMSDeltaPart(a,i,j,k)
l=round(i/j);
if (mod(i,j) ~= 0 | l>k) , x=0; else x=a(l); end;
end
function [x]=Power(a,b)
x=a^b;
end

end