#
#**************************
# SET UP THE INITIAL DATA *
#**************************
#   Problem :
#   *********
#   A minimum energy problem in quantum physics.
#   Source:
#   Leonor Cruzeiro-Hasson (private communication).
#   SIF input: Nick Gould, Feb 1991.
#   classification OOR2-MY-V-1
#   N is the number of discretization points.
#IE N                   10
#IE N                   50
#IE N                   100
#   problem constants
#   Define useful parameters
#   first contributions to the objective function
#   first contributions to the nonlinear constraint
#   general contributions to the objective function
#   general contributions to the nonlinear constraint
#   last contributions to the objective function
#   last contributions to the nonlinear constraint
#   Solution
#LO SOLTN(10)          -3.7565286
#LO SOLTN(50)          -4.2295940
#LO SOLTN(100)         -4.2887306
	param n := 200;
	param j := 1.549429;
	param k := 650.0;
	param x1 := 62.0;
	param x2 := 62.0;
	param nm1 := -1 + (200);
	param startp := sqrt(1.0 / (400.0));
	param a := -1.0 * (62.0);
	param b := (62.0) - (62.0);
	param c := 62.0;
	param ip1 := 1 + (199);
	param im1 := -1 + (199);

	var pr1 := 0.05;
	var pi1 := 0.05;
	var u1 := 0.1;
	var pr2 := 0.05;
	var pi2 := 0.05;
	var u2 := 0.2;
	var pr3 := 0.05;
	var pi3 := 0.05;
	var u3 := 0.3;
	var pr4 := 0.05;
	var pi4 := 0.05;
	var u4 := 0.4;
	var pr5 := 0.05;
	var pi5 := 0.05;
	var u5 := 0.5;
	var pr6 := 0.05;
	var pi6 := 0.05;
	var u6 := 0.6;
	var pr7 := 0.05;
	var pi7 := 0.05;
	var u7 := 0.7;
	var pr8 := 0.05;
	var pi8 := 0.05;
	var u8 := 0.7999999999999999;
	var pr9 := 0.05;
	var pi9 := 0.05;
	var u9 := 0.8999999999999999;
	var pr10 := 0.05;
	var pi10 := 0.05;
	var u10 := 0.9999999999999999;
	var pr11 := 0.05;
	var pi11 := 0.05;
	var u11 := 1.0999999999999999;
	var pr12 := 0.05;
	var pi12 := 0.05;
	var u12 := 1.2;
	var pr13 := 0.05;
	var pi13 := 0.05;
	var u13 := 1.3;
	var pr14 := 0.05;
	var pi14 := 0.05;
	var u14 := 1.4000000000000001;
	var pr15 := 0.05;
	var pi15 := 0.05;
	var u15 := 1.5000000000000002;
	var pr16 := 0.05;
	var pi16 := 0.05;
	var u16 := 1.6000000000000003;
	var pr17 := 0.05;
	var pi17 := 0.05;
	var u17 := 1.7000000000000004;
	var pr18 := 0.05;
	var pi18 := 0.05;
	var u18 := 1.8000000000000005;
	var pr19 := 0.05;
	var pi19 := 0.05;
	var u19 := 1.9000000000000006;
	var pr20 := 0.05;
	var pi20 := 0.05;
	var u20 := 2.0000000000000004;
	var pr21 := 0.05;
	var pi21 := 0.05;
	var u21 := 2.1000000000000005;
	var pr22 := 0.05;
	var pi22 := 0.05;
	var u22 := 2.2000000000000006;
	var pr23 := 0.05;
	var pi23 := 0.05;
	var u23 := 2.3000000000000007;
	var pr24 := 0.05;
	var pi24 := 0.05;
	var u24 := 2.400000000000001;
	var pr25 := 0.05;
	var pi25 := 0.05;
	var u25 := 2.500000000000001;
	var pr26 := 0.05;
	var pi26 := 0.05;
	var u26 := 2.600000000000001;
	var pr27 := 0.05;
	var pi27 := 0.05;
	var u27 := 2.700000000000001;
	var pr28 := 0.05;
	var pi28 := 0.05;
	var u28 := 2.800000000000001;
	var pr29 := 0.05;
	var pi29 := 0.05;
	var u29 := 2.9000000000000012;
	var pr30 := 0.05;
	var pi30 := 0.05;
	var u30 := 3.0000000000000013;
	var pr31 := 0.05;
	var pi31 := 0.05;
	var u31 := 3.1000000000000014;
	var pr32 := 0.05;
	var pi32 := 0.05;
	var u32 := 3.2000000000000015;
	var pr33 := 0.05;
	var pi33 := 0.05;
	var u33 := 3.3000000000000016;
	var pr34 := 0.05;
	var pi34 := 0.05;
	var u34 := 3.4000000000000017;
	var pr35 := 0.05;
	var pi35 := 0.05;
	var u35 := 3.5000000000000018;
	var pr36 := 0.05;
	var pi36 := 0.05;
	var u36 := 3.600000000000002;
	var pr37 := 0.05;
	var pi37 := 0.05;
	var u37 := 3.700000000000002;
	var pr38 := 0.05;
	var pi38 := 0.05;
	var u38 := 3.800000000000002;
	var pr39 := 0.05;
	var pi39 := 0.05;
	var u39 := 3.900000000000002;
	var pr40 := 0.05;
	var pi40 := 0.05;
	var u40 := 4.000000000000002;
	var pr41 := 0.05;
	var pi41 := 0.05;
	var u41 := 4.100000000000001;
	var pr42 := 0.05;
	var pi42 := 0.05;
	var u42 := 4.200000000000001;
	var pr43 := 0.05;
	var pi43 := 0.05;
	var u43 := 4.300000000000001;
	var pr44 := 0.05;
	var pi44 := 0.05;
	var u44 := 4.4;
	var pr45 := 0.05;
	var pi45 := 0.05;
	var u45 := 4.5;
	var pr46 := 0.05;
	var pi46 := 0.05;
	var u46 := 4.6;
	var pr47 := 0.05;
	var pi47 := 0.05;
	var u47 := 4.699999999999999;
	var pr48 := 0.05;
	var pi48 := 0.05;
	var u48 := 4.799999999999999;
	var pr49 := 0.05;
	var pi49 := 0.05;
	var u49 := 4.899999999999999;
	var pr50 := 0.05;
	var pi50 := 0.05;
	var u50 := 4.999999999999998;
	var pr51 := 0.05;
	var pi51 := 0.05;
	var u51 := 5.099999999999998;
	var pr52 := 0.05;
	var pi52 := 0.05;
	var u52 := 5.1999999999999975;
	var pr53 := 0.05;
	var pi53 := 0.05;
	var u53 := 5.299999999999997;
	var pr54 := 0.05;
	var pi54 := 0.05;
	var u54 := 5.399999999999997;
	var pr55 := 0.05;
	var pi55 := 0.05;
	var u55 := 5.4999999999999964;
	var pr56 := 0.05;
	var pi56 := 0.05;
	var u56 := 5.599999999999996;
	var pr57 := 0.05;
	var pi57 := 0.05;
	var u57 := 5.699999999999996;
	var pr58 := 0.05;
	var pi58 := 0.05;
	var u58 := 5.799999999999995;
	var pr59 := 0.05;
	var pi59 := 0.05;
	var u59 := 5.899999999999995;
	var pr60 := 0.05;
	var pi60 := 0.05;
	var u60 := 5.999999999999995;
	var pr61 := 0.05;
	var pi61 := 0.05;
	var u61 := 6.099999999999994;
	var pr62 := 0.05;
	var pi62 := 0.05;
	var u62 := 6.199999999999994;
	var pr63 := 0.05;
	var pi63 := 0.05;
	var u63 := 6.299999999999994;
	var pr64 := 0.05;
	var pi64 := 0.05;
	var u64 := 6.399999999999993;
	var pr65 := 0.05;
	var pi65 := 0.05;
	var u65 := 6.499999999999993;
	var pr66 := 0.05;
	var pi66 := 0.05;
	var u66 := 6.5999999999999925;
	var pr67 := 0.05;
	var pi67 := 0.05;
	var u67 := 6.699999999999992;
	var pr68 := 0.05;
	var pi68 := 0.05;
	var u68 := 6.799999999999992;
	var pr69 := 0.05;
	var pi69 := 0.05;
	var u69 := 6.8999999999999915;
	var pr70 := 0.05;
	var pi70 := 0.05;
	var u70 := 6.999999999999991;
	var pr71 := 0.05;
	var pi71 := 0.05;
	var u71 := 7.099999999999991;
	var pr72 := 0.05;
	var pi72 := 0.05;
	var u72 := 7.19999999999999;
	var pr73 := 0.05;
	var pi73 := 0.05;
	var u73 := 7.29999999999999;
	var pr74 := 0.05;
	var pi74 := 0.05;
	var u74 := 7.39999999999999;
	var pr75 := 0.05;
	var pi75 := 0.05;
	var u75 := 7.499999999999989;
	var pr76 := 0.05;
	var pi76 := 0.05;
	var u76 := 7.599999999999989;
	var pr77 := 0.05;
	var pi77 := 0.05;
	var u77 := 7.699999999999989;
	var pr78 := 0.05;
	var pi78 := 0.05;
	var u78 := 7.799999999999988;
	var pr79 := 0.05;
	var pi79 := 0.05;
	var u79 := 7.899999999999988;
	var pr80 := 0.05;
	var pi80 := 0.05;
	var u80 := 7.999999999999988;
	var pr81 := 0.05;
	var pi81 := 0.05;
	var u81 := 8.099999999999987;
	var pr82 := 0.05;
	var pi82 := 0.05;
	var u82 := 8.199999999999987;
	var pr83 := 0.05;
	var pi83 := 0.05;
	var u83 := 8.299999999999986;
	var pr84 := 0.05;
	var pi84 := 0.05;
	var u84 := 8.399999999999986;
	var pr85 := 0.05;
	var pi85 := 0.05;
	var u85 := 8.499999999999986;
	var pr86 := 0.05;
	var pi86 := 0.05;
	var u86 := 8.599999999999985;
	var pr87 := 0.05;
	var pi87 := 0.05;
	var u87 := 8.699999999999985;
	var pr88 := 0.05;
	var pi88 := 0.05;
	var u88 := 8.799999999999985;
	var pr89 := 0.05;
	var pi89 := 0.05;
	var u89 := 8.899999999999984;
	var pr90 := 0.05;
	var pi90 := 0.05;
	var u90 := 8.999999999999984;
	var pr91 := 0.05;
	var pi91 := 0.05;
	var u91 := 9.099999999999984;
	var pr92 := 0.05;
	var pi92 := 0.05;
	var u92 := 9.199999999999983;
	var pr93 := 0.05;
	var pi93 := 0.05;
	var u93 := 9.299999999999983;
	var pr94 := 0.05;
	var pi94 := 0.05;
	var u94 := 9.399999999999983;
	var pr95 := 0.05;
	var pi95 := 0.05;
	var u95 := 9.499999999999982;
	var pr96 := 0.05;
	var pi96 := 0.05;
	var u96 := 9.599999999999982;
	var pr97 := 0.05;
	var pi97 := 0.05;
	var u97 := 9.699999999999982;
	var pr98 := 0.05;
	var pi98 := 0.05;
	var u98 := 9.799999999999981;
	var pr99 := 0.05;
	var pi99 := 0.05;
	var u99 := 9.89999999999998;
	var pr100 := 0.05;
	var pi100 := 0.05;
	var u100 := 9.99999999999998;
	var pr101 := 0.05;
	var pi101 := 0.05;
	var u101 := 10.09999999999998;
	var pr102 := 0.05;
	var pi102 := 0.05;
	var u102 := 10.19999999999998;
	var pr103 := 0.05;
	var pi103 := 0.05;
	var u103 := 10.29999999999998;
	var pr104 := 0.05;
	var pi104 := 0.05;
	var u104 := 10.399999999999979;
	var pr105 := 0.05;
	var pi105 := 0.05;
	var u105 := 10.499999999999979;
	var pr106 := 0.05;
	var pi106 := 0.05;
	var u106 := 10.599999999999978;
	var pr107 := 0.05;
	var pi107 := 0.05;
	var u107 := 10.699999999999978;
	var pr108 := 0.05;
	var pi108 := 0.05;
	var u108 := 10.799999999999978;
	var pr109 := 0.05;
	var pi109 := 0.05;
	var u109 := 10.899999999999977;
	var pr110 := 0.05;
	var pi110 := 0.05;
	var u110 := 10.999999999999977;
	var pr111 := 0.05;
	var pi111 := 0.05;
	var u111 := 11.099999999999977;
	var pr112 := 0.05;
	var pi112 := 0.05;
	var u112 := 11.199999999999976;
	var pr113 := 0.05;
	var pi113 := 0.05;
	var u113 := 11.299999999999976;
	var pr114 := 0.05;
	var pi114 := 0.05;
	var u114 := 11.399999999999975;
	var pr115 := 0.05;
	var pi115 := 0.05;
	var u115 := 11.499999999999975;
	var pr116 := 0.05;
	var pi116 := 0.05;
	var u116 := 11.599999999999975;
	var pr117 := 0.05;
	var pi117 := 0.05;
	var u117 := 11.699999999999974;
	var pr118 := 0.05;
	var pi118 := 0.05;
	var u118 := 11.799999999999974;
	var pr119 := 0.05;
	var pi119 := 0.05;
	var u119 := 11.899999999999974;
	var pr120 := 0.05;
	var pi120 := 0.05;
	var u120 := 11.999999999999973;
	var pr121 := 0.05;
	var pi121 := 0.05;
	var u121 := 12.099999999999973;
	var pr122 := 0.05;
	var pi122 := 0.05;
	var u122 := 12.199999999999973;
	var pr123 := 0.05;
	var pi123 := 0.05;
	var u123 := 12.299999999999972;
	var pr124 := 0.05;
	var pi124 := 0.05;
	var u124 := 12.399999999999972;
	var pr125 := 0.05;
	var pi125 := 0.05;
	var u125 := 12.499999999999972;
	var pr126 := 0.05;
	var pi126 := 0.05;
	var u126 := 12.599999999999971;
	var pr127 := 0.05;
	var pi127 := 0.05;
	var u127 := 12.69999999999997;
	var pr128 := 0.05;
	var pi128 := 0.05;
	var u128 := 12.79999999999997;
	var pr129 := 0.05;
	var pi129 := 0.05;
	var u129 := 12.89999999999997;
	var pr130 := 0.05;
	var pi130 := 0.05;
	var u130 := 12.99999999999997;
	var pr131 := 0.05;
	var pi131 := 0.05;
	var u131 := 13.09999999999997;
	var pr132 := 0.05;
	var pi132 := 0.05;
	var u132 := 13.199999999999969;
	var pr133 := 0.05;
	var pi133 := 0.05;
	var u133 := 13.299999999999969;
	var pr134 := 0.05;
	var pi134 := 0.05;
	var u134 := 13.399999999999968;
	var pr135 := 0.05;
	var pi135 := 0.05;
	var u135 := 13.499999999999968;
	var pr136 := 0.05;
	var pi136 := 0.05;
	var u136 := 13.599999999999968;
	var pr137 := 0.05;
	var pi137 := 0.05;
	var u137 := 13.699999999999967;
	var pr138 := 0.05;
	var pi138 := 0.05;
	var u138 := 13.799999999999967;
	var pr139 := 0.05;
	var pi139 := 0.05;
	var u139 := 13.899999999999967;
	var pr140 := 0.05;
	var pi140 := 0.05;
	var u140 := 13.999999999999966;
	var pr141 := 0.05;
	var pi141 := 0.05;
	var u141 := 14.099999999999966;
	var pr142 := 0.05;
	var pi142 := 0.05;
	var u142 := 14.199999999999966;
	var pr143 := 0.05;
	var pi143 := 0.05;
	var u143 := 14.299999999999965;
	var pr144 := 0.05;
	var pi144 := 0.05;
	var u144 := 14.399999999999965;
	var pr145 := 0.05;
	var pi145 := 0.05;
	var u145 := 14.499999999999964;
	var pr146 := 0.05;
	var pi146 := 0.05;
	var u146 := 14.599999999999964;
	var pr147 := 0.05;
	var pi147 := 0.05;
	var u147 := 14.699999999999964;
	var pr148 := 0.05;
	var pi148 := 0.05;
	var u148 := 14.799999999999963;
	var pr149 := 0.05;
	var pi149 := 0.05;
	var u149 := 14.899999999999963;
	var pr150 := 0.05;
	var pi150 := 0.05;
	var u150 := 14.999999999999963;
	var pr151 := 0.05;
	var pi151 := 0.05;
	var u151 := 15.099999999999962;
	var pr152 := 0.05;
	var pi152 := 0.05;
	var u152 := 15.199999999999962;
	var pr153 := 0.05;
	var pi153 := 0.05;
	var u153 := 15.299999999999962;
	var pr154 := 0.05;
	var pi154 := 0.05;
	var u154 := 15.399999999999961;
	var pr155 := 0.05;
	var pi155 := 0.05;
	var u155 := 15.499999999999961;
	var pr156 := 0.05;
	var pi156 := 0.05;
	var u156 := 15.59999999999996;
	var pr157 := 0.05;
	var pi157 := 0.05;
	var u157 := 15.69999999999996;
	var pr158 := 0.05;
	var pi158 := 0.05;
	var u158 := 15.79999999999996;
	var pr159 := 0.05;
	var pi159 := 0.05;
	var u159 := 15.89999999999996;
	var pr160 := 0.05;
	var pi160 := 0.05;
	var u160 := 15.99999999999996;
	var pr161 := 0.05;
	var pi161 := 0.05;
	var u161 := 16.09999999999996;
	var pr162 := 0.05;
	var pi162 := 0.05;
	var u162 := 16.19999999999996;
	var pr163 := 0.05;
	var pi163 := 0.05;
	var u163 := 16.29999999999996;
	var pr164 := 0.05;
	var pi164 := 0.05;
	var u164 := 16.399999999999963;
	var pr165 := 0.05;
	var pi165 := 0.05;
	var u165 := 16.499999999999964;
	var pr166 := 0.05;
	var pi166 := 0.05;
	var u166 := 16.599999999999966;
	var pr167 := 0.05;
	var pi167 := 0.05;
	var u167 := 16.699999999999967;
	var pr168 := 0.05;
	var pi168 := 0.05;
	var u168 := 16.79999999999997;
	var pr169 := 0.05;
	var pi169 := 0.05;
	var u169 := 16.89999999999997;
	var pr170 := 0.05;
	var pi170 := 0.05;
	var u170 := 16.99999999999997;
	var pr171 := 0.05;
	var pi171 := 0.05;
	var u171 := 17.099999999999973;
	var pr172 := 0.05;
	var pi172 := 0.05;
	var u172 := 17.199999999999974;
	var pr173 := 0.05;
	var pi173 := 0.05;
	var u173 := 17.299999999999976;
	var pr174 := 0.05;
	var pi174 := 0.05;
	var u174 := 17.399999999999977;
	var pr175 := 0.05;
	var pi175 := 0.05;
	var u175 := 17.49999999999998;
	var pr176 := 0.05;
	var pi176 := 0.05;
	var u176 := 17.59999999999998;
	var pr177 := 0.05;
	var pi177 := 0.05;
	var u177 := 17.69999999999998;
	var pr178 := 0.05;
	var pi178 := 0.05;
	var u178 := 17.799999999999983;
	var pr179 := 0.05;
	var pi179 := 0.05;
	var u179 := 17.899999999999984;
	var pr180 := 0.05;
	var pi180 := 0.05;
	var u180 := 17.999999999999986;
	var pr181 := 0.05;
	var pi181 := 0.05;
	var u181 := 18.099999999999987;
	var pr182 := 0.05;
	var pi182 := 0.05;
	var u182 := 18.19999999999999;
	var pr183 := 0.05;
	var pi183 := 0.05;
	var u183 := 18.29999999999999;
	var pr184 := 0.05;
	var pi184 := 0.05;
	var u184 := 18.39999999999999;
	var pr185 := 0.05;
	var pi185 := 0.05;
	var u185 := 18.499999999999993;
	var pr186 := 0.05;
	var pi186 := 0.05;
	var u186 := 18.599999999999994;
	var pr187 := 0.05;
	var pi187 := 0.05;
	var u187 := 18.699999999999996;
	var pr188 := 0.05;
	var pi188 := 0.05;
	var u188 := 18.799999999999997;
	var pr189 := 0.05;
	var pi189 := 0.05;
	var u189 := 18.9;
	var pr190 := 0.05;
	var pi190 := 0.05;
	var u190 := 19.0;
	var pr191 := 0.05;
	var pi191 := 0.05;
	var u191 := 19.1;
	var pr192 := 0.05;
	var pi192 := 0.05;
	var u192 := 19.200000000000003;
	var pr193 := 0.05;
	var pi193 := 0.05;
	var u193 := 19.300000000000004;
	var pr194 := 0.05;
	var pi194 := 0.05;
	var u194 := 19.400000000000006;
	var pr195 := 0.05;
	var pi195 := 0.05;
	var u195 := 19.500000000000007;
	var pr196 := 0.05;
	var pi196 := 0.05;
	var u196 := 19.60000000000001;
	var pr197 := 0.05;
	var pi197 := 0.05;
	var u197 := 19.70000000000001;
	var pr198 := 0.05;
	var pi198 := 0.05;
	var u198 := 19.80000000000001;
	var pr199 := 0.05;
	var pi199 := 0.05;
	var u199 := 19.900000000000013;
	var pr200 := 0.05;
	var pi200 := 0.05;
	var u200 := 20.000000000000014;

minimize obj:
	 - 1.549429*(pr1) * (pr200+pr2) - 1.549429*(pi1) * (pi200+pi2) + 
	(-62.0*u200+0.0*u1+62.0*u2) * (pr1*pr1+pi1*pi1) + 650.0*(u200-u1) * (u200-u1) - 
	1.549429*(pr2) * (pr1+pr3) - 1.549429*(pi2) * (pi1+pi3) + 
	(-62.0*u1+0.0*u2+62.0*u3) * (pr2*pr2+pi2*pi2) + 650.0*(u1-u2) * (u1-u2) - 
	1.549429*(pr3) * (pr2+pr4) - 1.549429*(pi3) * (pi2+pi4) + 
	(-62.0*u2+0.0*u3+62.0*u4) * (pr3*pr3+pi3*pi3) + 650.0*(u2-u3) * (u2-u3) - 
	1.549429*(pr4) * (pr3+pr5) - 1.549429*(pi4) * (pi3+pi5) + 
	(-62.0*u3+0.0*u4+62.0*u5) * (pr4*pr4+pi4*pi4) + 650.0*(u3-u4) * (u3-u4) - 
	1.549429*(pr5) * (pr4+pr6) - 1.549429*(pi5) * (pi4+pi6) + 
	(-62.0*u4+0.0*u5+62.0*u6) * (pr5*pr5+pi5*pi5) + 650.0*(u4-u5) * (u4-u5) - 
	1.549429*(pr6) * (pr5+pr7) - 1.549429*(pi6) * (pi5+pi7) + 
	(-62.0*u5+0.0*u6+62.0*u7) * (pr6*pr6+pi6*pi6) + 650.0*(u5-u6) * (u5-u6) - 
	1.549429*(pr7) * (pr6+pr8) - 1.549429*(pi7) * (pi6+pi8) + 
	(-62.0*u6+0.0*u7+62.0*u8) * (pr7*pr7+pi7*pi7) + 650.0*(u6-u7) * (u6-u7) - 
	1.549429*(pr8) * (pr7+pr9) - 1.549429*(pi8) * (pi7+pi9) + 
	(-62.0*u7+0.0*u8+62.0*u9) * (pr8*pr8+pi8*pi8) + 650.0*(u7-u8) * (u7-u8) - 
	1.549429*(pr9) * (pr8+pr10) - 1.549429*(pi9) * (pi8+pi10) + 
	(-62.0*u8+0.0*u9+62.0*u10) * (pr9*pr9+pi9*pi9) + 650.0*(u8-u9) * (u8-u9) - 
	1.549429*(pr10) * (pr9+pr11) - 1.549429*(pi10) * (pi9+pi11) + 
	(-62.0*u9+0.0*u10+62.0*u11) * (pr10*pr10+pi10*pi10) + 650.0*(u9-u10) * (u9-u10) 
	- 1.549429*(pr11) * (pr10+pr12) - 1.549429*(pi11) * (pi10+pi12) + 
	(-62.0*u10+0.0*u11+62.0*u12) * (pr11*pr11+pi11*pi11) + 650.0*(u10-u11) * 
	(u10-u11) - 1.549429*(pr12) * (pr11+pr13) - 1.549429*(pi12) * (pi11+pi13) + 
	(-62.0*u11+0.0*u12+62.0*u13) * (pr12*pr12+pi12*pi12) + 650.0*(u11-u12) * 
	(u11-u12) - 1.549429*(pr13) * (pr12+pr14) - 1.549429*(pi13) * (pi12+pi14) + 
	(-62.0*u12+0.0*u13+62.0*u14) * (pr13*pr13+pi13*pi13) + 650.0*(u12-u13) * 
	(u12-u13) - 1.549429*(pr14) * (pr13+pr15) - 1.549429*(pi14) * (pi13+pi15) + 
	(-62.0*u13+0.0*u14+62.0*u15) * (pr14*pr14+pi14*pi14) + 650.0*(u13-u14) * 
	(u13-u14) - 1.549429*(pr15) * (pr14+pr16) - 1.549429*(pi15) * (pi14+pi16) + 
	(-62.0*u14+0.0*u15+62.0*u16) * (pr15*pr15+pi15*pi15) + 650.0*(u14-u15) * 
	(u14-u15) - 1.549429*(pr16) * (pr15+pr17) - 1.549429*(pi16) * (pi15+pi17) + 
	(-62.0*u15+0.0*u16+62.0*u17) * (pr16*pr16+pi16*pi16) + 650.0*(u15-u16) * 
	(u15-u16) - 1.549429*(pr17) * (pr16+pr18) - 1.549429*(pi17) * (pi16+pi18) + 
	(-62.0*u16+0.0*u17+62.0*u18) * (pr17*pr17+pi17*pi17) + 650.0*(u16-u17) * 
	(u16-u17) - 1.549429*(pr18) * (pr17+pr19) - 1.549429*(pi18) * (pi17+pi19) + 
	(-62.0*u17+0.0*u18+62.0*u19) * (pr18*pr18+pi18*pi18) + 650.0*(u17-u18) * 
	(u17-u18) - 1.549429*(pr19) * (pr18+pr20) - 1.549429*(pi19) * (pi18+pi20) + 
	(-62.0*u18+0.0*u19+62.0*u20) * (pr19*pr19+pi19*pi19) + 650.0*(u18-u19) * 
	(u18-u19) - 1.549429*(pr20) * (pr19+pr21) - 1.549429*(pi20) * (pi19+pi21) + 
	(-62.0*u19+0.0*u20+62.0*u21) * (pr20*pr20+pi20*pi20) + 650.0*(u19-u20) * 
	(u19-u20) - 1.549429*(pr21) * (pr20+pr22) - 1.549429*(pi21) * (pi20+pi22) + 
	(-62.0*u20+0.0*u21+62.0*u22) * (pr21*pr21+pi21*pi21) + 650.0*(u20-u21) * 
	(u20-u21) - 1.549429*(pr22) * (pr21+pr23) - 1.549429*(pi22) * (pi21+pi23) + 
	(-62.0*u21+0.0*u22+62.0*u23) * (pr22*pr22+pi22*pi22) + 650.0*(u21-u22) * 
	(u21-u22) - 1.549429*(pr23) * (pr22+pr24) - 1.549429*(pi23) * (pi22+pi24) + 
	(-62.0*u22+0.0*u23+62.0*u24) * (pr23*pr23+pi23*pi23) + 650.0*(u22-u23) * 
	(u22-u23) - 1.549429*(pr24) * (pr23+pr25) - 1.549429*(pi24) * (pi23+pi25) + 
	(-62.0*u23+0.0*u24+62.0*u25) * (pr24*pr24+pi24*pi24) + 650.0*(u23-u24) * 
	(u23-u24) - 1.549429*(pr25) * (pr24+pr26) - 1.549429*(pi25) * (pi24+pi26) + 
	(-62.0*u24+0.0*u25+62.0*u26) * (pr25*pr25+pi25*pi25) + 650.0*(u24-u25) * 
	(u24-u25) - 1.549429*(pr26) * (pr25+pr27) - 1.549429*(pi26) * (pi25+pi27) + 
	(-62.0*u25+0.0*u26+62.0*u27) * (pr26*pr26+pi26*pi26) + 650.0*(u25-u26) * 
	(u25-u26) - 1.549429*(pr27) * (pr26+pr28) - 1.549429*(pi27) * (pi26+pi28) + 
	(-62.0*u26+0.0*u27+62.0*u28) * (pr27*pr27+pi27*pi27) + 650.0*(u26-u27) * 
	(u26-u27) - 1.549429*(pr28) * (pr27+pr29) - 1.549429*(pi28) * (pi27+pi29) + 
	(-62.0*u27+0.0*u28+62.0*u29) * (pr28*pr28+pi28*pi28) + 650.0*(u27-u28) * 
	(u27-u28) - 1.549429*(pr29) * (pr28+pr30) - 1.549429*(pi29) * (pi28+pi30) + 
	(-62.0*u28+0.0*u29+62.0*u30) * (pr29*pr29+pi29*pi29) + 650.0*(u28-u29) * 
	(u28-u29) - 1.549429*(pr30) * (pr29+pr31) - 1.549429*(pi30) * (pi29+pi31) + 
	(-62.0*u29+0.0*u30+62.0*u31) * (pr30*pr30+pi30*pi30) + 650.0*(u29-u30) * 
	(u29-u30) - 1.549429*(pr31) * (pr30+pr32) - 1.549429*(pi31) * (pi30+pi32) + 
	(-62.0*u30+0.0*u31+62.0*u32) * (pr31*pr31+pi31*pi31) + 650.0*(u30-u31) * 
	(u30-u31) - 1.549429*(pr32) * (pr31+pr33) - 1.549429*(pi32) * (pi31+pi33) + 
	(-62.0*u31+0.0*u32+62.0*u33) * (pr32*pr32+pi32*pi32) + 650.0*(u31-u32) * 
	(u31-u32) - 1.549429*(pr33) * (pr32+pr34) - 1.549429*(pi33) * (pi32+pi34) + 
	(-62.0*u32+0.0*u33+62.0*u34) * (pr33*pr33+pi33*pi33) + 650.0*(u32-u33) * 
	(u32-u33) - 1.549429*(pr34) * (pr33+pr35) - 1.549429*(pi34) * (pi33+pi35) + 
	(-62.0*u33+0.0*u34+62.0*u35) * (pr34*pr34+pi34*pi34) + 650.0*(u33-u34) * 
	(u33-u34) - 1.549429*(pr35) * (pr34+pr36) - 1.549429*(pi35) * (pi34+pi36) + 
	(-62.0*u34+0.0*u35+62.0*u36) * (pr35*pr35+pi35*pi35) + 650.0*(u34-u35) * 
	(u34-u35) - 1.549429*(pr36) * (pr35+pr37) - 1.549429*(pi36) * (pi35+pi37) + 
	(-62.0*u35+0.0*u36+62.0*u37) * (pr36*pr36+pi36*pi36) + 650.0*(u35-u36) * 
	(u35-u36) - 1.549429*(pr37) * (pr36+pr38) - 1.549429*(pi37) * (pi36+pi38) + 
	(-62.0*u36+0.0*u37+62.0*u38) * (pr37*pr37+pi37*pi37) + 650.0*(u36-u37) * 
	(u36-u37) - 1.549429*(pr38) * (pr37+pr39) - 1.549429*(pi38) * (pi37+pi39) + 
	(-62.0*u37+0.0*u38+62.0*u39) * (pr38*pr38+pi38*pi38) + 650.0*(u37-u38) * 
	(u37-u38) - 1.549429*(pr39) * (pr38+pr40) - 1.549429*(pi39) * (pi38+pi40) + 
	(-62.0*u38+0.0*u39+62.0*u40) * (pr39*pr39+pi39*pi39) + 650.0*(u38-u39) * 
	(u38-u39) - 1.549429*(pr40) * (pr39+pr41) - 1.549429*(pi40) * (pi39+pi41) + 
	(-62.0*u39+0.0*u40+62.0*u41) * (pr40*pr40+pi40*pi40) + 650.0*(u39-u40) * 
	(u39-u40) - 1.549429*(pr41) * (pr40+pr42) - 1.549429*(pi41) * (pi40+pi42) + 
	(-62.0*u40+0.0*u41+62.0*u42) * (pr41*pr41+pi41*pi41) + 650.0*(u40-u41) * 
	(u40-u41) - 1.549429*(pr42) * (pr41+pr43) - 1.549429*(pi42) * (pi41+pi43) + 
	(-62.0*u41+0.0*u42+62.0*u43) * (pr42*pr42+pi42*pi42) + 650.0*(u41-u42) * 
	(u41-u42) - 1.549429*(pr43) * (pr42+pr44) - 1.549429*(pi43) * (pi42+pi44) + 
	(-62.0*u42+0.0*u43+62.0*u44) * (pr43*pr43+pi43*pi43) + 650.0*(u42-u43) * 
	(u42-u43) - 1.549429*(pr44) * (pr43+pr45) - 1.549429*(pi44) * (pi43+pi45) + 
	(-62.0*u43+0.0*u44+62.0*u45) * (pr44*pr44+pi44*pi44) + 650.0*(u43-u44) * 
	(u43-u44) - 1.549429*(pr45) * (pr44+pr46) - 1.549429*(pi45) * (pi44+pi46) + 
	(-62.0*u44+0.0*u45+62.0*u46) * (pr45*pr45+pi45*pi45) + 650.0*(u44-u45) * 
	(u44-u45) - 1.549429*(pr46) * (pr45+pr47) - 1.549429*(pi46) * (pi45+pi47) + 
	(-62.0*u45+0.0*u46+62.0*u47) * (pr46*pr46+pi46*pi46) + 650.0*(u45-u46) * 
	(u45-u46) - 1.549429*(pr47) * (pr46+pr48) - 1.549429*(pi47) * (pi46+pi48) + 
	(-62.0*u46+0.0*u47+62.0*u48) * (pr47*pr47+pi47*pi47) + 650.0*(u46-u47) * 
	(u46-u47) - 1.549429*(pr48) * (pr47+pr49) - 1.549429*(pi48) * (pi47+pi49) + 
	(-62.0*u47+0.0*u48+62.0*u49) * (pr48*pr48+pi48*pi48) + 650.0*(u47-u48) * 
	(u47-u48) - 1.549429*(pr49) * (pr48+pr50) - 1.549429*(pi49) * (pi48+pi50) + 
	(-62.0*u48+0.0*u49+62.0*u50) * (pr49*pr49+pi49*pi49) + 650.0*(u48-u49) * 
	(u48-u49) - 1.549429*(pr50) * (pr49+pr51) - 1.549429*(pi50) * (pi49+pi51) + 
	(-62.0*u49+0.0*u50+62.0*u51) * (pr50*pr50+pi50*pi50) + 650.0*(u49-u50) * 
	(u49-u50) - 1.549429*(pr51) * (pr50+pr52) - 1.549429*(pi51) * (pi50+pi52) + 
	(-62.0*u50+0.0*u51+62.0*u52) * (pr51*pr51+pi51*pi51) + 650.0*(u50-u51) * 
	(u50-u51) - 1.549429*(pr52) * (pr51+pr53) - 1.549429*(pi52) * (pi51+pi53) + 
	(-62.0*u51+0.0*u52+62.0*u53) * (pr52*pr52+pi52*pi52) + 650.0*(u51-u52) * 
	(u51-u52) - 1.549429*(pr53) * (pr52+pr54) - 1.549429*(pi53) * (pi52+pi54) + 
	(-62.0*u52+0.0*u53+62.0*u54) * (pr53*pr53+pi53*pi53) + 650.0*(u52-u53) * 
	(u52-u53) - 1.549429*(pr54) * (pr53+pr55) - 1.549429*(pi54) * (pi53+pi55) + 
	(-62.0*u53+0.0*u54+62.0*u55) * (pr54*pr54+pi54*pi54) + 650.0*(u53-u54) * 
	(u53-u54) - 1.549429*(pr55) * (pr54+pr56) - 1.549429*(pi55) * (pi54+pi56) + 
	(-62.0*u54+0.0*u55+62.0*u56) * (pr55*pr55+pi55*pi55) + 650.0*(u54-u55) * 
	(u54-u55) - 1.549429*(pr56) * (pr55+pr57) - 1.549429*(pi56) * (pi55+pi57) + 
	(-62.0*u55+0.0*u56+62.0*u57) * (pr56*pr56+pi56*pi56) + 650.0*(u55-u56) * 
	(u55-u56) - 1.549429*(pr57) * (pr56+pr58) - 1.549429*(pi57) * (pi56+pi58) + 
	(-62.0*u56+0.0*u57+62.0*u58) * (pr57*pr57+pi57*pi57) + 650.0*(u56-u57) * 
	(u56-u57) - 1.549429*(pr58) * (pr57+pr59) - 1.549429*(pi58) * (pi57+pi59) + 
	(-62.0*u57+0.0*u58+62.0*u59) * (pr58*pr58+pi58*pi58) + 650.0*(u57-u58) * 
	(u57-u58) - 1.549429*(pr59) * (pr58+pr60) - 1.549429*(pi59) * (pi58+pi60) + 
	(-62.0*u58+0.0*u59+62.0*u60) * (pr59*pr59+pi59*pi59) + 650.0*(u58-u59) * 
	(u58-u59) - 1.549429*(pr60) * (pr59+pr61) - 1.549429*(pi60) * (pi59+pi61) + 
	(-62.0*u59+0.0*u60+62.0*u61) * (pr60*pr60+pi60*pi60) + 650.0*(u59-u60) * 
	(u59-u60) - 1.549429*(pr61) * (pr60+pr62) - 1.549429*(pi61) * (pi60+pi62) + 
	(-62.0*u60+0.0*u61+62.0*u62) * (pr61*pr61+pi61*pi61) + 650.0*(u60-u61) * 
	(u60-u61) - 1.549429*(pr62) * (pr61+pr63) - 1.549429*(pi62) * (pi61+pi63) + 
	(-62.0*u61+0.0*u62+62.0*u63) * (pr62*pr62+pi62*pi62) + 650.0*(u61-u62) * 
	(u61-u62) - 1.549429*(pr63) * (pr62+pr64) - 1.549429*(pi63) * (pi62+pi64) + 
	(-62.0*u62+0.0*u63+62.0*u64) * (pr63*pr63+pi63*pi63) + 650.0*(u62-u63) * 
	(u62-u63) - 1.549429*(pr64) * (pr63+pr65) - 1.549429*(pi64) * (pi63+pi65) + 
	(-62.0*u63+0.0*u64+62.0*u65) * (pr64*pr64+pi64*pi64) + 650.0*(u63-u64) * 
	(u63-u64) - 1.549429*(pr65) * (pr64+pr66) - 1.549429*(pi65) * (pi64+pi66) + 
	(-62.0*u64+0.0*u65+62.0*u66) * (pr65*pr65+pi65*pi65) + 650.0*(u64-u65) * 
	(u64-u65) - 1.549429*(pr66) * (pr65+pr67) - 1.549429*(pi66) * (pi65+pi67) + 
	(-62.0*u65+0.0*u66+62.0*u67) * (pr66*pr66+pi66*pi66) + 650.0*(u65-u66) * 
	(u65-u66) - 1.549429*(pr67) * (pr66+pr68) - 1.549429*(pi67) * (pi66+pi68) + 
	(-62.0*u66+0.0*u67+62.0*u68) * (pr67*pr67+pi67*pi67) + 650.0*(u66-u67) * 
	(u66-u67) - 1.549429*(pr68) * (pr67+pr69) - 1.549429*(pi68) * (pi67+pi69) + 
	(-62.0*u67+0.0*u68+62.0*u69) * (pr68*pr68+pi68*pi68) + 650.0*(u67-u68) * 
	(u67-u68) - 1.549429*(pr69) * (pr68+pr70) - 1.549429*(pi69) * (pi68+pi70) + 
	(-62.0*u68+0.0*u69+62.0*u70) * (pr69*pr69+pi69*pi69) + 650.0*(u68-u69) * 
	(u68-u69) - 1.549429*(pr70) * (pr69+pr71) - 1.549429*(pi70) * (pi69+pi71) + 
	(-62.0*u69+0.0*u70+62.0*u71) * (pr70*pr70+pi70*pi70) + 650.0*(u69-u70) * 
	(u69-u70) - 1.549429*(pr71) * (pr70+pr72) - 1.549429*(pi71) * (pi70+pi72) + 
	(-62.0*u70+0.0*u71+62.0*u72) * (pr71*pr71+pi71*pi71) + 650.0*(u70-u71) * 
	(u70-u71) - 1.549429*(pr72) * (pr71+pr73) - 1.549429*(pi72) * (pi71+pi73) + 
	(-62.0*u71+0.0*u72+62.0*u73) * (pr72*pr72+pi72*pi72) + 650.0*(u71-u72) * 
	(u71-u72) - 1.549429*(pr73) * (pr72+pr74) - 1.549429*(pi73) * (pi72+pi74) + 
	(-62.0*u72+0.0*u73+62.0*u74) * (pr73*pr73+pi73*pi73) + 650.0*(u72-u73) * 
	(u72-u73) - 1.549429*(pr74) * (pr73+pr75) - 1.549429*(pi74) * (pi73+pi75) + 
	(-62.0*u73+0.0*u74+62.0*u75) * (pr74*pr74+pi74*pi74) + 650.0*(u73-u74) * 
	(u73-u74) - 1.549429*(pr75) * (pr74+pr76) - 1.549429*(pi75) * (pi74+pi76) + 
	(-62.0*u74+0.0*u75+62.0*u76) * (pr75*pr75+pi75*pi75) + 650.0*(u74-u75) * 
	(u74-u75) - 1.549429*(pr76) * (pr75+pr77) - 1.549429*(pi76) * (pi75+pi77) + 
	(-62.0*u75+0.0*u76+62.0*u77) * (pr76*pr76+pi76*pi76) + 650.0*(u75-u76) * 
	(u75-u76) - 1.549429*(pr77) * (pr76+pr78) - 1.549429*(pi77) * (pi76+pi78) + 
	(-62.0*u76+0.0*u77+62.0*u78) * (pr77*pr77+pi77*pi77) + 650.0*(u76-u77) * 
	(u76-u77) - 1.549429*(pr78) * (pr77+pr79) - 1.549429*(pi78) * (pi77+pi79) + 
	(-62.0*u77+0.0*u78+62.0*u79) * (pr78*pr78+pi78*pi78) + 650.0*(u77-u78) * 
	(u77-u78) - 1.549429*(pr79) * (pr78+pr80) - 1.549429*(pi79) * (pi78+pi80) + 
	(-62.0*u78+0.0*u79+62.0*u80) * (pr79*pr79+pi79*pi79) + 650.0*(u78-u79) * 
	(u78-u79) - 1.549429*(pr80) * (pr79+pr81) - 1.549429*(pi80) * (pi79+pi81) + 
	(-62.0*u79+0.0*u80+62.0*u81) * (pr80*pr80+pi80*pi80) + 650.0*(u79-u80) * 
	(u79-u80) - 1.549429*(pr81) * (pr80+pr82) - 1.549429*(pi81) * (pi80+pi82) + 
	(-62.0*u80+0.0*u81+62.0*u82) * (pr81*pr81+pi81*pi81) + 650.0*(u80-u81) * 
	(u80-u81) - 1.549429*(pr82) * (pr81+pr83) - 1.549429*(pi82) * (pi81+pi83) + 
	(-62.0*u81+0.0*u82+62.0*u83) * (pr82*pr82+pi82*pi82) + 650.0*(u81-u82) * 
	(u81-u82) - 1.549429*(pr83) * (pr82+pr84) - 1.549429*(pi83) * (pi82+pi84) + 
	(-62.0*u82+0.0*u83+62.0*u84) * (pr83*pr83+pi83*pi83) + 650.0*(u82-u83) * 
	(u82-u83) - 1.549429*(pr84) * (pr83+pr85) - 1.549429*(pi84) * (pi83+pi85) + 
	(-62.0*u83+0.0*u84+62.0*u85) * (pr84*pr84+pi84*pi84) + 650.0*(u83-u84) * 
	(u83-u84) - 1.549429*(pr85) * (pr84+pr86) - 1.549429*(pi85) * (pi84+pi86) + 
	(-62.0*u84+0.0*u85+62.0*u86) * (pr85*pr85+pi85*pi85) + 650.0*(u84-u85) * 
	(u84-u85) - 1.549429*(pr86) * (pr85+pr87) - 1.549429*(pi86) * (pi85+pi87) + 
	(-62.0*u85+0.0*u86+62.0*u87) * (pr86*pr86+pi86*pi86) + 650.0*(u85-u86) * 
	(u85-u86) - 1.549429*(pr87) * (pr86+pr88) - 1.549429*(pi87) * (pi86+pi88) + 
	(-62.0*u86+0.0*u87+62.0*u88) * (pr87*pr87+pi87*pi87) + 650.0*(u86-u87) * 
	(u86-u87) - 1.549429*(pr88) * (pr87+pr89) - 1.549429*(pi88) * (pi87+pi89) + 
	(-62.0*u87+0.0*u88+62.0*u89) * (pr88*pr88+pi88*pi88) + 650.0*(u87-u88) * 
	(u87-u88) - 1.549429*(pr89) * (pr88+pr90) - 1.549429*(pi89) * (pi88+pi90) + 
	(-62.0*u88+0.0*u89+62.0*u90) * (pr89*pr89+pi89*pi89) + 650.0*(u88-u89) * 
	(u88-u89) - 1.549429*(pr90) * (pr89+pr91) - 1.549429*(pi90) * (pi89+pi91) + 
	(-62.0*u89+0.0*u90+62.0*u91) * (pr90*pr90+pi90*pi90) + 650.0*(u89-u90) * 
	(u89-u90) - 1.549429*(pr91) * (pr90+pr92) - 1.549429*(pi91) * (pi90+pi92) + 
	(-62.0*u90+0.0*u91+62.0*u92) * (pr91*pr91+pi91*pi91) + 650.0*(u90-u91) * 
	(u90-u91) - 1.549429*(pr92) * (pr91+pr93) - 1.549429*(pi92) * (pi91+pi93) + 
	(-62.0*u91+0.0*u92+62.0*u93) * (pr92*pr92+pi92*pi92) + 650.0*(u91-u92) * 
	(u91-u92) - 1.549429*(pr93) * (pr92+pr94) - 1.549429*(pi93) * (pi92+pi94) + 
	(-62.0*u92+0.0*u93+62.0*u94) * (pr93*pr93+pi93*pi93) + 650.0*(u92-u93) * 
	(u92-u93) - 1.549429*(pr94) * (pr93+pr95) - 1.549429*(pi94) * (pi93+pi95) + 
	(-62.0*u93+0.0*u94+62.0*u95) * (pr94*pr94+pi94*pi94) + 650.0*(u93-u94) * 
	(u93-u94) - 1.549429*(pr95) * (pr94+pr96) - 1.549429*(pi95) * (pi94+pi96) + 
	(-62.0*u94+0.0*u95+62.0*u96) * (pr95*pr95+pi95*pi95) + 650.0*(u94-u95) * 
	(u94-u95) - 1.549429*(pr96) * (pr95+pr97) - 1.549429*(pi96) * (pi95+pi97) + 
	(-62.0*u95+0.0*u96+62.0*u97) * (pr96*pr96+pi96*pi96) + 650.0*(u95-u96) * 
	(u95-u96) - 1.549429*(pr97) * (pr96+pr98) - 1.549429*(pi97) * (pi96+pi98) + 
	(-62.0*u96+0.0*u97+62.0*u98) * (pr97*pr97+pi97*pi97) + 650.0*(u96-u97) * 
	(u96-u97) - 1.549429*(pr98) * (pr97+pr99) - 1.549429*(pi98) * (pi97+pi99) + 
	(-62.0*u97+0.0*u98+62.0*u99) * (pr98*pr98+pi98*pi98) + 650.0*(u97-u98) * 
	(u97-u98) - 1.549429*(pr99) * (pr98+pr100) - 1.549429*(pi99) * (pi98+pi100) + 
	(-62.0*u98+0.0*u99+62.0*u100) * (pr99*pr99+pi99*pi99) + 650.0*(u98-u99) * 
	(u98-u99) - 1.549429*(pr100) * (pr99+pr101) - 1.549429*(pi100) * (pi99+pi101) + 
	(-62.0*u99+0.0*u100+62.0*u101) * (pr100*pr100+pi100*pi100) + 650.0*(u99-u100) * 
	(u99-u100) - 1.549429*(pr101) * (pr100+pr102) - 1.549429*(pi101) * 
	(pi100+pi102) + (-62.0*u100+0.0*u101+62.0*u102) * (pr101*pr101+pi101*pi101) + 
	650.0*(u100-u101) * (u100-u101) - 1.549429*(pr102) * (pr101+pr103) - 
	1.549429*(pi102) * (pi101+pi103) + (-62.0*u101+0.0*u102+62.0*u103) * 
	(pr102*pr102+pi102*pi102) + 650.0*(u101-u102) * (u101-u102) - 1.549429*(pr103) 
	* (pr102+pr104) - 1.549429*(pi103) * (pi102+pi104) + 
	(-62.0*u102+0.0*u103+62.0*u104) * (pr103*pr103+pi103*pi103) + 650.0*(u102-u103) 
	* (u102-u103) - 1.549429*(pr104) * (pr103+pr105) - 1.549429*(pi104) * 
	(pi103+pi105) + (-62.0*u103+0.0*u104+62.0*u105) * (pr104*pr104+pi104*pi104) + 
	650.0*(u103-u104) * (u103-u104) - 1.549429*(pr105) * (pr104+pr106) - 
	1.549429*(pi105) * (pi104+pi106) + (-62.0*u104+0.0*u105+62.0*u106) * 
	(pr105*pr105+pi105*pi105) + 650.0*(u104-u105) * (u104-u105) - 1.549429*(pr106) 
	* (pr105+pr107) - 1.549429*(pi106) * (pi105+pi107) + 
	(-62.0*u105+0.0*u106+62.0*u107) * (pr106*pr106+pi106*pi106) + 650.0*(u105-u106) 
	* (u105-u106) - 1.549429*(pr107) * (pr106+pr108) - 1.549429*(pi107) * 
	(pi106+pi108) + (-62.0*u106+0.0*u107+62.0*u108) * (pr107*pr107+pi107*pi107) + 
	650.0*(u106-u107) * (u106-u107) - 1.549429*(pr108) * (pr107+pr109) - 
	1.549429*(pi108) * (pi107+pi109) + (-62.0*u107+0.0*u108+62.0*u109) * 
	(pr108*pr108+pi108*pi108) + 650.0*(u107-u108) * (u107-u108) - 1.549429*(pr109) 
	* (pr108+pr110) - 1.549429*(pi109) * (pi108+pi110) + 
	(-62.0*u108+0.0*u109+62.0*u110) * (pr109*pr109+pi109*pi109) + 650.0*(u108-u109) 
	* (u108-u109) - 1.549429*(pr110) * (pr109+pr111) - 1.549429*(pi110) * 
	(pi109+pi111) + (-62.0*u109+0.0*u110+62.0*u111) * (pr110*pr110+pi110*pi110) + 
	650.0*(u109-u110) * (u109-u110) - 1.549429*(pr111) * (pr110+pr112) - 
	1.549429*(pi111) * (pi110+pi112) + (-62.0*u110+0.0*u111+62.0*u112) * 
	(pr111*pr111+pi111*pi111) + 650.0*(u110-u111) * (u110-u111) - 1.549429*(pr112) 
	* (pr111+pr113) - 1.549429*(pi112) * (pi111+pi113) + 
	(-62.0*u111+0.0*u112+62.0*u113) * (pr112*pr112+pi112*pi112) + 650.0*(u111-u112) 
	* (u111-u112) - 1.549429*(pr113) * (pr112+pr114) - 1.549429*(pi113) * 
	(pi112+pi114) + (-62.0*u112+0.0*u113+62.0*u114) * (pr113*pr113+pi113*pi113) + 
	650.0*(u112-u113) * (u112-u113) - 1.549429*(pr114) * (pr113+pr115) - 
	1.549429*(pi114) * (pi113+pi115) + (-62.0*u113+0.0*u114+62.0*u115) * 
	(pr114*pr114+pi114*pi114) + 650.0*(u113-u114) * (u113-u114) - 1.549429*(pr115) 
	* (pr114+pr116) - 1.549429*(pi115) * (pi114+pi116) + 
	(-62.0*u114+0.0*u115+62.0*u116) * (pr115*pr115+pi115*pi115) + 650.0*(u114-u115) 
	* (u114-u115) - 1.549429*(pr116) * (pr115+pr117) - 1.549429*(pi116) * 
	(pi115+pi117) + (-62.0*u115+0.0*u116+62.0*u117) * (pr116*pr116+pi116*pi116) + 
	650.0*(u115-u116) * (u115-u116) - 1.549429*(pr117) * (pr116+pr118) - 
	1.549429*(pi117) * (pi116+pi118) + (-62.0*u116+0.0*u117+62.0*u118) * 
	(pr117*pr117+pi117*pi117) + 650.0*(u116-u117) * (u116-u117) - 1.549429*(pr118) 
	* (pr117+pr119) - 1.549429*(pi118) * (pi117+pi119) + 
	(-62.0*u117+0.0*u118+62.0*u119) * (pr118*pr118+pi118*pi118) + 650.0*(u117-u118) 
	* (u117-u118) - 1.549429*(pr119) * (pr118+pr120) - 1.549429*(pi119) * 
	(pi118+pi120) + (-62.0*u118+0.0*u119+62.0*u120) * (pr119*pr119+pi119*pi119) + 
	650.0*(u118-u119) * (u118-u119) - 1.549429*(pr120) * (pr119+pr121) - 
	1.549429*(pi120) * (pi119+pi121) + (-62.0*u119+0.0*u120+62.0*u121) * 
	(pr120*pr120+pi120*pi120) + 650.0*(u119-u120) * (u119-u120) - 1.549429*(pr121) 
	* (pr120+pr122) - 1.549429*(pi121) * (pi120+pi122) + 
	(-62.0*u120+0.0*u121+62.0*u122) * (pr121*pr121+pi121*pi121) + 650.0*(u120-u121) 
	* (u120-u121) - 1.549429*(pr122) * (pr121+pr123) - 1.549429*(pi122) * 
	(pi121+pi123) + (-62.0*u121+0.0*u122+62.0*u123) * (pr122*pr122+pi122*pi122) + 
	650.0*(u121-u122) * (u121-u122) - 1.549429*(pr123) * (pr122+pr124) - 
	1.549429*(pi123) * (pi122+pi124) + (-62.0*u122+0.0*u123+62.0*u124) * 
	(pr123*pr123+pi123*pi123) + 650.0*(u122-u123) * (u122-u123) - 1.549429*(pr124) 
	* (pr123+pr125) - 1.549429*(pi124) * (pi123+pi125) + 
	(-62.0*u123+0.0*u124+62.0*u125) * (pr124*pr124+pi124*pi124) + 650.0*(u123-u124) 
	* (u123-u124) - 1.549429*(pr125) * (pr124+pr126) - 1.549429*(pi125) * 
	(pi124+pi126) + (-62.0*u124+0.0*u125+62.0*u126) * (pr125*pr125+pi125*pi125) + 
	650.0*(u124-u125) * (u124-u125) - 1.549429*(pr126) * (pr125+pr127) - 
	1.549429*(pi126) * (pi125+pi127) + (-62.0*u125+0.0*u126+62.0*u127) * 
	(pr126*pr126+pi126*pi126) + 650.0*(u125-u126) * (u125-u126) - 1.549429*(pr127) 
	* (pr126+pr128) - 1.549429*(pi127) * (pi126+pi128) + 
	(-62.0*u126+0.0*u127+62.0*u128) * (pr127*pr127+pi127*pi127) + 650.0*(u126-u127) 
	* (u126-u127) - 1.549429*(pr128) * (pr127+pr129) - 1.549429*(pi128) * 
	(pi127+pi129) + (-62.0*u127+0.0*u128+62.0*u129) * (pr128*pr128+pi128*pi128) + 
	650.0*(u127-u128) * (u127-u128) - 1.549429*(pr129) * (pr128+pr130) - 
	1.549429*(pi129) * (pi128+pi130) + (-62.0*u128+0.0*u129+62.0*u130) * 
	(pr129*pr129+pi129*pi129) + 650.0*(u128-u129) * (u128-u129) - 1.549429*(pr130) 
	* (pr129+pr131) - 1.549429*(pi130) * (pi129+pi131) + 
	(-62.0*u129+0.0*u130+62.0*u131) * (pr130*pr130+pi130*pi130) + 650.0*(u129-u130) 
	* (u129-u130) - 1.549429*(pr131) * (pr130+pr132) - 1.549429*(pi131) * 
	(pi130+pi132) + (-62.0*u130+0.0*u131+62.0*u132) * (pr131*pr131+pi131*pi131) + 
	650.0*(u130-u131) * (u130-u131) - 1.549429*(pr132) * (pr131+pr133) - 
	1.549429*(pi132) * (pi131+pi133) + (-62.0*u131+0.0*u132+62.0*u133) * 
	(pr132*pr132+pi132*pi132) + 650.0*(u131-u132) * (u131-u132) - 1.549429*(pr133) 
	* (pr132+pr134) - 1.549429*(pi133) * (pi132+pi134) + 
	(-62.0*u132+0.0*u133+62.0*u134) * (pr133*pr133+pi133*pi133) + 650.0*(u132-u133) 
	* (u132-u133) - 1.549429*(pr134) * (pr133+pr135) - 1.549429*(pi134) * 
	(pi133+pi135) + (-62.0*u133+0.0*u134+62.0*u135) * (pr134*pr134+pi134*pi134) + 
	650.0*(u133-u134) * (u133-u134) - 1.549429*(pr135) * (pr134+pr136) - 
	1.549429*(pi135) * (pi134+pi136) + (-62.0*u134+0.0*u135+62.0*u136) * 
	(pr135*pr135+pi135*pi135) + 650.0*(u134-u135) * (u134-u135) - 1.549429*(pr136) 
	* (pr135+pr137) - 1.549429*(pi136) * (pi135+pi137) + 
	(-62.0*u135+0.0*u136+62.0*u137) * (pr136*pr136+pi136*pi136) + 650.0*(u135-u136) 
	* (u135-u136) - 1.549429*(pr137) * (pr136+pr138) - 1.549429*(pi137) * 
	(pi136+pi138) + (-62.0*u136+0.0*u137+62.0*u138) * (pr137*pr137+pi137*pi137) + 
	650.0*(u136-u137) * (u136-u137) - 1.549429*(pr138) * (pr137+pr139) - 
	1.549429*(pi138) * (pi137+pi139) + (-62.0*u137+0.0*u138+62.0*u139) * 
	(pr138*pr138+pi138*pi138) + 650.0*(u137-u138) * (u137-u138) - 1.549429*(pr139) 
	* (pr138+pr140) - 1.549429*(pi139) * (pi138+pi140) + 
	(-62.0*u138+0.0*u139+62.0*u140) * (pr139*pr139+pi139*pi139) + 650.0*(u138-u139) 
	* (u138-u139) - 1.549429*(pr140) * (pr139+pr141) - 1.549429*(pi140) * 
	(pi139+pi141) + (-62.0*u139+0.0*u140+62.0*u141) * (pr140*pr140+pi140*pi140) + 
	650.0*(u139-u140) * (u139-u140) - 1.549429*(pr141) * (pr140+pr142) - 
	1.549429*(pi141) * (pi140+pi142) + (-62.0*u140+0.0*u141+62.0*u142) * 
	(pr141*pr141+pi141*pi141) + 650.0*(u140-u141) * (u140-u141) - 1.549429*(pr142) 
	* (pr141+pr143) - 1.549429*(pi142) * (pi141+pi143) + 
	(-62.0*u141+0.0*u142+62.0*u143) * (pr142*pr142+pi142*pi142) + 650.0*(u141-u142) 
	* (u141-u142) - 1.549429*(pr143) * (pr142+pr144) - 1.549429*(pi143) * 
	(pi142+pi144) + (-62.0*u142+0.0*u143+62.0*u144) * (pr143*pr143+pi143*pi143) + 
	650.0*(u142-u143) * (u142-u143) - 1.549429*(pr144) * (pr143+pr145) - 
	1.549429*(pi144) * (pi143+pi145) + (-62.0*u143+0.0*u144+62.0*u145) * 
	(pr144*pr144+pi144*pi144) + 650.0*(u143-u144) * (u143-u144) - 1.549429*(pr145) 
	* (pr144+pr146) - 1.549429*(pi145) * (pi144+pi146) + 
	(-62.0*u144+0.0*u145+62.0*u146) * (pr145*pr145+pi145*pi145) + 650.0*(u144-u145) 
	* (u144-u145) - 1.549429*(pr146) * (pr145+pr147) - 1.549429*(pi146) * 
	(pi145+pi147) + (-62.0*u145+0.0*u146+62.0*u147) * (pr146*pr146+pi146*pi146) + 
	650.0*(u145-u146) * (u145-u146) - 1.549429*(pr147) * (pr146+pr148) - 
	1.549429*(pi147) * (pi146+pi148) + (-62.0*u146+0.0*u147+62.0*u148) * 
	(pr147*pr147+pi147*pi147) + 650.0*(u146-u147) * (u146-u147) - 1.549429*(pr148) 
	* (pr147+pr149) - 1.549429*(pi148) * (pi147+pi149) + 
	(-62.0*u147+0.0*u148+62.0*u149) * (pr148*pr148+pi148*pi148) + 650.0*(u147-u148) 
	* (u147-u148) - 1.549429*(pr149) * (pr148+pr150) - 1.549429*(pi149) * 
	(pi148+pi150) + (-62.0*u148+0.0*u149+62.0*u150) * (pr149*pr149+pi149*pi149) + 
	650.0*(u148-u149) * (u148-u149) - 1.549429*(pr150) * (pr149+pr151) - 
	1.549429*(pi150) * (pi149+pi151) + (-62.0*u149+0.0*u150+62.0*u151) * 
	(pr150*pr150+pi150*pi150) + 650.0*(u149-u150) * (u149-u150) - 1.549429*(pr151) 
	* (pr150+pr152) - 1.549429*(pi151) * (pi150+pi152) + 
	(-62.0*u150+0.0*u151+62.0*u152) * (pr151*pr151+pi151*pi151) + 650.0*(u150-u151) 
	* (u150-u151) - 1.549429*(pr152) * (pr151+pr153) - 1.549429*(pi152) * 
	(pi151+pi153) + (-62.0*u151+0.0*u152+62.0*u153) * (pr152*pr152+pi152*pi152) + 
	650.0*(u151-u152) * (u151-u152) - 1.549429*(pr153) * (pr152+pr154) - 
	1.549429*(pi153) * (pi152+pi154) + (-62.0*u152+0.0*u153+62.0*u154) * 
	(pr153*pr153+pi153*pi153) + 650.0*(u152-u153) * (u152-u153) - 1.549429*(pr154) 
	* (pr153+pr155) - 1.549429*(pi154) * (pi153+pi155) + 
	(-62.0*u153+0.0*u154+62.0*u155) * (pr154*pr154+pi154*pi154) + 650.0*(u153-u154) 
	* (u153-u154) - 1.549429*(pr155) * (pr154+pr156) - 1.549429*(pi155) * 
	(pi154+pi156) + (-62.0*u154+0.0*u155+62.0*u156) * (pr155*pr155+pi155*pi155) + 
	650.0*(u154-u155) * (u154-u155) - 1.549429*(pr156) * (pr155+pr157) - 
	1.549429*(pi156) * (pi155+pi157) + (-62.0*u155+0.0*u156+62.0*u157) * 
	(pr156*pr156+pi156*pi156) + 650.0*(u155-u156) * (u155-u156) - 1.549429*(pr157) 
	* (pr156+pr158) - 1.549429*(pi157) * (pi156+pi158) + 
	(-62.0*u156+0.0*u157+62.0*u158) * (pr157*pr157+pi157*pi157) + 650.0*(u156-u157) 
	* (u156-u157) - 1.549429*(pr158) * (pr157+pr159) - 1.549429*(pi158) * 
	(pi157+pi159) + (-62.0*u157+0.0*u158+62.0*u159) * (pr158*pr158+pi158*pi158) + 
	650.0*(u157-u158) * (u157-u158) - 1.549429*(pr159) * (pr158+pr160) - 
	1.549429*(pi159) * (pi158+pi160) + (-62.0*u158+0.0*u159+62.0*u160) * 
	(pr159*pr159+pi159*pi159) + 650.0*(u158-u159) * (u158-u159) - 1.549429*(pr160) 
	* (pr159+pr161) - 1.549429*(pi160) * (pi159+pi161) + 
	(-62.0*u159+0.0*u160+62.0*u161) * (pr160*pr160+pi160*pi160) + 650.0*(u159-u160) 
	* (u159-u160) - 1.549429*(pr161) * (pr160+pr162) - 1.549429*(pi161) * 
	(pi160+pi162) + (-62.0*u160+0.0*u161+62.0*u162) * (pr161*pr161+pi161*pi161) + 
	650.0*(u160-u161) * (u160-u161) - 1.549429*(pr162) * (pr161+pr163) - 
	1.549429*(pi162) * (pi161+pi163) + (-62.0*u161+0.0*u162+62.0*u163) * 
	(pr162*pr162+pi162*pi162) + 650.0*(u161-u162) * (u161-u162) - 1.549429*(pr163) 
	* (pr162+pr164) - 1.549429*(pi163) * (pi162+pi164) + 
	(-62.0*u162+0.0*u163+62.0*u164) * (pr163*pr163+pi163*pi163) + 650.0*(u162-u163) 
	* (u162-u163) - 1.549429*(pr164) * (pr163+pr165) - 1.549429*(pi164) * 
	(pi163+pi165) + (-62.0*u163+0.0*u164+62.0*u165) * (pr164*pr164+pi164*pi164) + 
	650.0*(u163-u164) * (u163-u164) - 1.549429*(pr165) * (pr164+pr166) - 
	1.549429*(pi165) * (pi164+pi166) + (-62.0*u164+0.0*u165+62.0*u166) * 
	(pr165*pr165+pi165*pi165) + 650.0*(u164-u165) * (u164-u165) - 1.549429*(pr166) 
	* (pr165+pr167) - 1.549429*(pi166) * (pi165+pi167) + 
	(-62.0*u165+0.0*u166+62.0*u167) * (pr166*pr166+pi166*pi166) + 650.0*(u165-u166) 
	* (u165-u166) - 1.549429*(pr167) * (pr166+pr168) - 1.549429*(pi167) * 
	(pi166+pi168) + (-62.0*u166+0.0*u167+62.0*u168) * (pr167*pr167+pi167*pi167) + 
	650.0*(u166-u167) * (u166-u167) - 1.549429*(pr168) * (pr167+pr169) - 
	1.549429*(pi168) * (pi167+pi169) + (-62.0*u167+0.0*u168+62.0*u169) * 
	(pr168*pr168+pi168*pi168) + 650.0*(u167-u168) * (u167-u168) - 1.549429*(pr169) 
	* (pr168+pr170) - 1.549429*(pi169) * (pi168+pi170) + 
	(-62.0*u168+0.0*u169+62.0*u170) * (pr169*pr169+pi169*pi169) + 650.0*(u168-u169) 
	* (u168-u169) - 1.549429*(pr170) * (pr169+pr171) - 1.549429*(pi170) * 
	(pi169+pi171) + (-62.0*u169+0.0*u170+62.0*u171) * (pr170*pr170+pi170*pi170) + 
	650.0*(u169-u170) * (u169-u170) - 1.549429*(pr171) * (pr170+pr172) - 
	1.549429*(pi171) * (pi170+pi172) + (-62.0*u170+0.0*u171+62.0*u172) * 
	(pr171*pr171+pi171*pi171) + 650.0*(u170-u171) * (u170-u171) - 1.549429*(pr172) 
	* (pr171+pr173) - 1.549429*(pi172) * (pi171+pi173) + 
	(-62.0*u171+0.0*u172+62.0*u173) * (pr172*pr172+pi172*pi172) + 650.0*(u171-u172) 
	* (u171-u172) - 1.549429*(pr173) * (pr172+pr174) - 1.549429*(pi173) * 
	(pi172+pi174) + (-62.0*u172+0.0*u173+62.0*u174) * (pr173*pr173+pi173*pi173) + 
	650.0*(u172-u173) * (u172-u173) - 1.549429*(pr174) * (pr173+pr175) - 
	1.549429*(pi174) * (pi173+pi175) + (-62.0*u173+0.0*u174+62.0*u175) * 
	(pr174*pr174+pi174*pi174) + 650.0*(u173-u174) * (u173-u174) - 1.549429*(pr175) 
	* (pr174+pr176) - 1.549429*(pi175) * (pi174+pi176) + 
	(-62.0*u174+0.0*u175+62.0*u176) * (pr175*pr175+pi175*pi175) + 650.0*(u174-u175) 
	* (u174-u175) - 1.549429*(pr176) * (pr175+pr177) - 1.549429*(pi176) * 
	(pi175+pi177) + (-62.0*u175+0.0*u176+62.0*u177) * (pr176*pr176+pi176*pi176) + 
	650.0*(u175-u176) * (u175-u176) - 1.549429*(pr177) * (pr176+pr178) - 
	1.549429*(pi177) * (pi176+pi178) + (-62.0*u176+0.0*u177+62.0*u178) * 
	(pr177*pr177+pi177*pi177) + 650.0*(u176-u177) * (u176-u177) - 1.549429*(pr178) 
	* (pr177+pr179) - 1.549429*(pi178) * (pi177+pi179) + 
	(-62.0*u177+0.0*u178+62.0*u179) * (pr178*pr178+pi178*pi178) + 650.0*(u177-u178) 
	* (u177-u178) - 1.549429*(pr179) * (pr178+pr180) - 1.549429*(pi179) * 
	(pi178+pi180) + (-62.0*u178+0.0*u179+62.0*u180) * (pr179*pr179+pi179*pi179) + 
	650.0*(u178-u179) * (u178-u179) - 1.549429*(pr180) * (pr179+pr181) - 
	1.549429*(pi180) * (pi179+pi181) + (-62.0*u179+0.0*u180+62.0*u181) * 
	(pr180*pr180+pi180*pi180) + 650.0*(u179-u180) * (u179-u180) - 1.549429*(pr181) 
	* (pr180+pr182) - 1.549429*(pi181) * (pi180+pi182) + 
	(-62.0*u180+0.0*u181+62.0*u182) * (pr181*pr181+pi181*pi181) + 650.0*(u180-u181) 
	* (u180-u181) - 1.549429*(pr182) * (pr181+pr183) - 1.549429*(pi182) * 
	(pi181+pi183) + (-62.0*u181+0.0*u182+62.0*u183) * (pr182*pr182+pi182*pi182) + 
	650.0*(u181-u182) * (u181-u182) - 1.549429*(pr183) * (pr182+pr184) - 
	1.549429*(pi183) * (pi182+pi184) + (-62.0*u182+0.0*u183+62.0*u184) * 
	(pr183*pr183+pi183*pi183) + 650.0*(u182-u183) * (u182-u183) - 1.549429*(pr184) 
	* (pr183+pr185) - 1.549429*(pi184) * (pi183+pi185) + 
	(-62.0*u183+0.0*u184+62.0*u185) * (pr184*pr184+pi184*pi184) + 650.0*(u183-u184) 
	* (u183-u184) - 1.549429*(pr185) * (pr184+pr186) - 1.549429*(pi185) * 
	(pi184+pi186) + (-62.0*u184+0.0*u185+62.0*u186) * (pr185*pr185+pi185*pi185) + 
	650.0*(u184-u185) * (u184-u185) - 1.549429*(pr186) * (pr185+pr187) - 
	1.549429*(pi186) * (pi185+pi187) + (-62.0*u185+0.0*u186+62.0*u187) * 
	(pr186*pr186+pi186*pi186) + 650.0*(u185-u186) * (u185-u186) - 1.549429*(pr187) 
	* (pr186+pr188) - 1.549429*(pi187) * (pi186+pi188) + 
	(-62.0*u186+0.0*u187+62.0*u188) * (pr187*pr187+pi187*pi187) + 650.0*(u186-u187) 
	* (u186-u187) - 1.549429*(pr188) * (pr187+pr189) - 1.549429*(pi188) * 
	(pi187+pi189) + (-62.0*u187+0.0*u188+62.0*u189) * (pr188*pr188+pi188*pi188) + 
	650.0*(u187-u188) * (u187-u188) - 1.549429*(pr189) * (pr188+pr190) - 
	1.549429*(pi189) * (pi188+pi190) + (-62.0*u188+0.0*u189+62.0*u190) * 
	(pr189*pr189+pi189*pi189) + 650.0*(u188-u189) * (u188-u189) - 1.549429*(pr190) 
	* (pr189+pr191) - 1.549429*(pi190) * (pi189+pi191) + 
	(-62.0*u189+0.0*u190+62.0*u191) * (pr190*pr190+pi190*pi190) + 650.0*(u189-u190) 
	* (u189-u190) - 1.549429*(pr191) * (pr190+pr192) - 1.549429*(pi191) * 
	(pi190+pi192) + (-62.0*u190+0.0*u191+62.0*u192) * (pr191*pr191+pi191*pi191) + 
	650.0*(u190-u191) * (u190-u191) - 1.549429*(pr192) * (pr191+pr193) - 
	1.549429*(pi192) * (pi191+pi193) + (-62.0*u191+0.0*u192+62.0*u193) * 
	(pr192*pr192+pi192*pi192) + 650.0*(u191-u192) * (u191-u192) - 1.549429*(pr193) 
	* (pr192+pr194) - 1.549429*(pi193) * (pi192+pi194) + 
	(-62.0*u192+0.0*u193+62.0*u194) * (pr193*pr193+pi193*pi193) + 650.0*(u192-u193) 
	* (u192-u193) - 1.549429*(pr194) * (pr193+pr195) - 1.549429*(pi194) * 
	(pi193+pi195) + (-62.0*u193+0.0*u194+62.0*u195) * (pr194*pr194+pi194*pi194) + 
	650.0*(u193-u194) * (u193-u194) - 1.549429*(pr195) * (pr194+pr196) - 
	1.549429*(pi195) * (pi194+pi196) + (-62.0*u194+0.0*u195+62.0*u196) * 
	(pr195*pr195+pi195*pi195) + 650.0*(u194-u195) * (u194-u195) - 1.549429*(pr196) 
	* (pr195+pr197) - 1.549429*(pi196) * (pi195+pi197) + 
	(-62.0*u195+0.0*u196+62.0*u197) * (pr196*pr196+pi196*pi196) + 650.0*(u195-u196) 
	* (u195-u196) - 1.549429*(pr197) * (pr196+pr198) - 1.549429*(pi197) * 
	(pi196+pi198) + (-62.0*u196+0.0*u197+62.0*u198) * (pr197*pr197+pi197*pi197) + 
	650.0*(u196-u197) * (u196-u197) - 1.549429*(pr198) * (pr197+pr199) - 
	1.549429*(pi198) * (pi197+pi199) + (-62.0*u197+0.0*u198+62.0*u199) * 
	(pr198*pr198+pi198*pi198) + 650.0*(u197-u198) * (u197-u198) - 1.549429*(pr199) 
	* (pr198+pr200) - 1.549429*(pi199) * (pi198+pi200) + 
	(-62.0*u198+0.0*u199+62.0*u200) * (pr199*pr199+pi199*pi199) + 650.0*(u198-u199) 
	* (u198-u199) - 1.549429*(pr200) * (pr199+pr1) - 1.549429*(pi200) * (pi199+pi1) 
	+ (-62.0*u199+0.0*u200+62.0*u1) * (pr200*pr200+pi200*pi200) + 650.0*(u199-u200) 
	* (u199-u200);

subject to con:
	pr1 * pr1 + pi1 * pi1 + pr2 * pr2 + pi2 * pi2 + pr3 * pr3 + pi3 * pi3 + pr4 * 
	pr4 + pi4 * pi4 + pr5 * pr5 + pi5 * pi5 + pr6 * pr6 + pi6 * pi6 + pr7 * pr7 + 
	pi7 * pi7 + pr8 * pr8 + pi8 * pi8 + pr9 * pr9 + pi9 * pi9 + pr10 * pr10 + pi10 
	* pi10 + pr11 * pr11 + pi11 * pi11 + pr12 * pr12 + pi12 * pi12 + pr13 * pr13 + 
	pi13 * pi13 + pr14 * pr14 + pi14 * pi14 + pr15 * pr15 + pi15 * pi15 + pr16 * 
	pr16 + pi16 * pi16 + pr17 * pr17 + pi17 * pi17 + pr18 * pr18 + pi18 * pi18 + 
	pr19 * pr19 + pi19 * pi19 + pr20 * pr20 + pi20 * pi20 + pr21 * pr21 + pi21 * 
	pi21 + pr22 * pr22 + pi22 * pi22 + pr23 * pr23 + pi23 * pi23 + pr24 * pr24 + 
	pi24 * pi24 + pr25 * pr25 + pi25 * pi25 + pr26 * pr26 + pi26 * pi26 + pr27 * 
	pr27 + pi27 * pi27 + pr28 * pr28 + pi28 * pi28 + pr29 * pr29 + pi29 * pi29 + 
	pr30 * pr30 + pi30 * pi30 + pr31 * pr31 + pi31 * pi31 + pr32 * pr32 + pi32 * 
	pi32 + pr33 * pr33 + pi33 * pi33 + pr34 * pr34 + pi34 * pi34 + pr35 * pr35 + 
	pi35 * pi35 + pr36 * pr36 + pi36 * pi36 + pr37 * pr37 + pi37 * pi37 + pr38 * 
	pr38 + pi38 * pi38 + pr39 * pr39 + pi39 * pi39 + pr40 * pr40 + pi40 * pi40 + 
	pr41 * pr41 + pi41 * pi41 + pr42 * pr42 + pi42 * pi42 + pr43 * pr43 + pi43 * 
	pi43 + pr44 * pr44 + pi44 * pi44 + pr45 * pr45 + pi45 * pi45 + pr46 * pr46 + 
	pi46 * pi46 + pr47 * pr47 + pi47 * pi47 + pr48 * pr48 + pi48 * pi48 + pr49 * 
	pr49 + pi49 * pi49 + pr50 * pr50 + pi50 * pi50 + pr51 * pr51 + pi51 * pi51 + 
	pr52 * pr52 + pi52 * pi52 + pr53 * pr53 + pi53 * pi53 + pr54 * pr54 + pi54 * 
	pi54 + pr55 * pr55 + pi55 * pi55 + pr56 * pr56 + pi56 * pi56 + pr57 * pr57 + 
	pi57 * pi57 + pr58 * pr58 + pi58 * pi58 + pr59 * pr59 + pi59 * pi59 + pr60 * 
	pr60 + pi60 * pi60 + pr61 * pr61 + pi61 * pi61 + pr62 * pr62 + pi62 * pi62 + 
	pr63 * pr63 + pi63 * pi63 + pr64 * pr64 + pi64 * pi64 + pr65 * pr65 + pi65 * 
	pi65 + pr66 * pr66 + pi66 * pi66 + pr67 * pr67 + pi67 * pi67 + pr68 * pr68 + 
	pi68 * pi68 + pr69 * pr69 + pi69 * pi69 + pr70 * pr70 + pi70 * pi70 + pr71 * 
	pr71 + pi71 * pi71 + pr72 * pr72 + pi72 * pi72 + pr73 * pr73 + pi73 * pi73 + 
	pr74 * pr74 + pi74 * pi74 + pr75 * pr75 + pi75 * pi75 + pr76 * pr76 + pi76 * 
	pi76 + pr77 * pr77 + pi77 * pi77 + pr78 * pr78 + pi78 * pi78 + pr79 * pr79 + 
	pi79 * pi79 + pr80 * pr80 + pi80 * pi80 + pr81 * pr81 + pi81 * pi81 + pr82 * 
	pr82 + pi82 * pi82 + pr83 * pr83 + pi83 * pi83 + pr84 * pr84 + pi84 * pi84 + 
	pr85 * pr85 + pi85 * pi85 + pr86 * pr86 + pi86 * pi86 + pr87 * pr87 + pi87 * 
	pi87 + pr88 * pr88 + pi88 * pi88 + pr89 * pr89 + pi89 * pi89 + pr90 * pr90 + 
	pi90 * pi90 + pr91 * pr91 + pi91 * pi91 + pr92 * pr92 + pi92 * pi92 + pr93 * 
	pr93 + pi93 * pi93 + pr94 * pr94 + pi94 * pi94 + pr95 * pr95 + pi95 * pi95 + 
	pr96 * pr96 + pi96 * pi96 + pr97 * pr97 + pi97 * pi97 + pr98 * pr98 + pi98 * 
	pi98 + pr99 * pr99 + pi99 * pi99 + pr100 * pr100 + pi100 * pi100 + pr101 * 
	pr101 + pi101 * pi101 + pr102 * pr102 + pi102 * pi102 + pr103 * pr103 + pi103 * 
	pi103 + pr104 * pr104 + pi104 * pi104 + pr105 * pr105 + pi105 * pi105 + pr106 * 
	pr106 + pi106 * pi106 + pr107 * pr107 + pi107 * pi107 + pr108 * pr108 + pi108 * 
	pi108 + pr109 * pr109 + pi109 * pi109 + pr110 * pr110 + pi110 * pi110 + pr111 * 
	pr111 + pi111 * pi111 + pr112 * pr112 + pi112 * pi112 + pr113 * pr113 + pi113 * 
	pi113 + pr114 * pr114 + pi114 * pi114 + pr115 * pr115 + pi115 * pi115 + pr116 * 
	pr116 + pi116 * pi116 + pr117 * pr117 + pi117 * pi117 + pr118 * pr118 + pi118 * 
	pi118 + pr119 * pr119 + pi119 * pi119 + pr120 * pr120 + pi120 * pi120 + pr121 * 
	pr121 + pi121 * pi121 + pr122 * pr122 + pi122 * pi122 + pr123 * pr123 + pi123 * 
	pi123 + pr124 * pr124 + pi124 * pi124 + pr125 * pr125 + pi125 * pi125 + pr126 * 
	pr126 + pi126 * pi126 + pr127 * pr127 + pi127 * pi127 + pr128 * pr128 + pi128 * 
	pi128 + pr129 * pr129 + pi129 * pi129 + pr130 * pr130 + pi130 * pi130 + pr131 * 
	pr131 + pi131 * pi131 + pr132 * pr132 + pi132 * pi132 + pr133 * pr133 + pi133 * 
	pi133 + pr134 * pr134 + pi134 * pi134 + pr135 * pr135 + pi135 * pi135 + pr136 * 
	pr136 + pi136 * pi136 + pr137 * pr137 + pi137 * pi137 + pr138 * pr138 + pi138 * 
	pi138 + pr139 * pr139 + pi139 * pi139 + pr140 * pr140 + pi140 * pi140 + pr141 * 
	pr141 + pi141 * pi141 + pr142 * pr142 + pi142 * pi142 + pr143 * pr143 + pi143 * 
	pi143 + pr144 * pr144 + pi144 * pi144 + pr145 * pr145 + pi145 * pi145 + pr146 * 
	pr146 + pi146 * pi146 + pr147 * pr147 + pi147 * pi147 + pr148 * pr148 + pi148 * 
	pi148 + pr149 * pr149 + pi149 * pi149 + pr150 * pr150 + pi150 * pi150 + pr151 * 
	pr151 + pi151 * pi151 + pr152 * pr152 + pi152 * pi152 + pr153 * pr153 + pi153 * 
	pi153 + pr154 * pr154 + pi154 * pi154 + pr155 * pr155 + pi155 * pi155 + pr156 * 
	pr156 + pi156 * pi156 + pr157 * pr157 + pi157 * pi157 + pr158 * pr158 + pi158 * 
	pi158 + pr159 * pr159 + pi159 * pi159 + pr160 * pr160 + pi160 * pi160 + pr161 * 
	pr161 + pi161 * pi161 + pr162 * pr162 + pi162 * pi162 + pr163 * pr163 + pi163 * 
	pi163 + pr164 * pr164 + pi164 * pi164 + pr165 * pr165 + pi165 * pi165 + pr166 * 
	pr166 + pi166 * pi166 + pr167 * pr167 + pi167 * pi167 + pr168 * pr168 + pi168 * 
	pi168 + pr169 * pr169 + pi169 * pi169 + pr170 * pr170 + pi170 * pi170 + pr171 * 
	pr171 + pi171 * pi171 + pr172 * pr172 + pi172 * pi172 + pr173 * pr173 + pi173 * 
	pi173 + pr174 * pr174 + pi174 * pi174 + pr175 * pr175 + pi175 * pi175 + pr176 * 
	pr176 + pi176 * pi176 + pr177 * pr177 + pi177 * pi177 + pr178 * pr178 + pi178 * 
	pi178 + pr179 * pr179 + pi179 * pi179 + pr180 * pr180 + pi180 * pi180 + pr181 * 
	pr181 + pi181 * pi181 + pr182 * pr182 + pi182 * pi182 + pr183 * pr183 + pi183 * 
	pi183 + pr184 * pr184 + pi184 * pi184 + pr185 * pr185 + pi185 * pi185 + pr186 * 
	pr186 + pi186 * pi186 + pr187 * pr187 + pi187 * pi187 + pr188 * pr188 + pi188 * 
	pi188 + pr189 * pr189 + pi189 * pi189 + pr190 * pr190 + pi190 * pi190 + pr191 * 
	pr191 + pi191 * pi191 + pr192 * pr192 + pi192 * pi192 + pr193 * pr193 + pi193 * 
	pi193 + pr194 * pr194 + pi194 * pi194 + pr195 * pr195 + pi195 * pi195 + pr196 * 
	pr196 + pi196 * pi196 + pr197 * pr197 + pi197 * pi197 + pr198 * pr198 + pi198 * 
	pi198 + pr199 * pr199 + pi199 * pi199 + pr200 * pr200 + pi200 * pi200 - 1.0 = 0;

solve;
	display pr1;
	display pi1;
	display u1;
	display pr2;
	display pi2;
	display u2;
	display pr3;
	display pi3;
	display u3;
	display pr4;
	display pi4;
	display u4;
	display pr5;
	display pi5;
	display u5;
	display pr6;
	display pi6;
	display u6;
	display pr7;
	display pi7;
	display u7;
	display pr8;
	display pi8;
	display u8;
	display pr9;
	display pi9;
	display u9;
	display pr10;
	display pi10;
	display u10;
	display pr11;
	display pi11;
	display u11;
	display pr12;
	display pi12;
	display u12;
	display pr13;
	display pi13;
	display u13;
	display pr14;
	display pi14;
	display u14;
	display pr15;
	display pi15;
	display u15;
	display pr16;
	display pi16;
	display u16;
	display pr17;
	display pi17;
	display u17;
	display pr18;
	display pi18;
	display u18;
	display pr19;
	display pi19;
	display u19;
	display pr20;
	display pi20;
	display u20;
	display pr21;
	display pi21;
	display u21;
	display pr22;
	display pi22;
	display u22;
	display pr23;
	display pi23;
	display u23;
	display pr24;
	display pi24;
	display u24;
	display pr25;
	display pi25;
	display u25;
	display pr26;
	display pi26;
	display u26;
	display pr27;
	display pi27;
	display u27;
	display pr28;
	display pi28;
	display u28;
	display pr29;
	display pi29;
	display u29;
	display pr30;
	display pi30;
	display u30;
	display pr31;
	display pi31;
	display u31;
	display pr32;
	display pi32;
	display u32;
	display pr33;
	display pi33;
	display u33;
	display pr34;
	display pi34;
	display u34;
	display pr35;
	display pi35;
	display u35;
	display pr36;
	display pi36;
	display u36;
	display pr37;
	display pi37;
	display u37;
	display pr38;
	display pi38;
	display u38;
	display pr39;
	display pi39;
	display u39;
	display pr40;
	display pi40;
	display u40;
	display pr41;
	display pi41;
	display u41;
	display pr42;
	display pi42;
	display u42;
	display pr43;
	display pi43;
	display u43;
	display pr44;
	display pi44;
	display u44;
	display pr45;
	display pi45;
	display u45;
	display pr46;
	display pi46;
	display u46;
	display pr47;
	display pi47;
	display u47;
	display pr48;
	display pi48;
	display u48;
	display pr49;
	display pi49;
	display u49;
	display pr50;
	display pi50;
	display u50;
	display pr51;
	display pi51;
	display u51;
	display pr52;
	display pi52;
	display u52;
	display pr53;
	display pi53;
	display u53;
	display pr54;
	display pi54;
	display u54;
	display pr55;
	display pi55;
	display u55;
	display pr56;
	display pi56;
	display u56;
	display pr57;
	display pi57;
	display u57;
	display pr58;
	display pi58;
	display u58;
	display pr59;
	display pi59;
	display u59;
	display pr60;
	display pi60;
	display u60;
	display pr61;
	display pi61;
	display u61;
	display pr62;
	display pi62;
	display u62;
	display pr63;
	display pi63;
	display u63;
	display pr64;
	display pi64;
	display u64;
	display pr65;
	display pi65;
	display u65;
	display pr66;
	display pi66;
	display u66;
	display pr67;
	display pi67;
	display u67;
	display pr68;
	display pi68;
	display u68;
	display pr69;
	display pi69;
	display u69;
	display pr70;
	display pi70;
	display u70;
	display pr71;
	display pi71;
	display u71;
	display pr72;
	display pi72;
	display u72;
	display pr73;
	display pi73;
	display u73;
	display pr74;
	display pi74;
	display u74;
	display pr75;
	display pi75;
	display u75;
	display pr76;
	display pi76;
	display u76;
	display pr77;
	display pi77;
	display u77;
	display pr78;
	display pi78;
	display u78;
	display pr79;
	display pi79;
	display u79;
	display pr80;
	display pi80;
	display u80;
	display pr81;
	display pi81;
	display u81;
	display pr82;
	display pi82;
	display u82;
	display pr83;
	display pi83;
	display u83;
	display pr84;
	display pi84;
	display u84;
	display pr85;
	display pi85;
	display u85;
	display pr86;
	display pi86;
	display u86;
	display pr87;
	display pi87;
	display u87;
	display pr88;
	display pi88;
	display u88;
	display pr89;
	display pi89;
	display u89;
	display pr90;
	display pi90;
	display u90;
	display pr91;
	display pi91;
	display u91;
	display pr92;
	display pi92;
	display u92;
	display pr93;
	display pi93;
	display u93;
	display pr94;
	display pi94;
	display u94;
	display pr95;
	display pi95;
	display u95;
	display pr96;
	display pi96;
	display u96;
	display pr97;
	display pi97;
	display u97;
	display pr98;
	display pi98;
	display u98;
	display pr99;
	display pi99;
	display u99;
	display pr100;
	display pi100;
	display u100;
	display pr101;
	display pi101;
	display u101;
	display pr102;
	display pi102;
	display u102;
	display pr103;
	display pi103;
	display u103;
	display pr104;
	display pi104;
	display u104;
	display pr105;
	display pi105;
	display u105;
	display pr106;
	display pi106;
	display u106;
	display pr107;
	display pi107;
	display u107;
	display pr108;
	display pi108;
	display u108;
	display pr109;
	display pi109;
	display u109;
	display pr110;
	display pi110;
	display u110;
	display pr111;
	display pi111;
	display u111;
	display pr112;
	display pi112;
	display u112;
	display pr113;
	display pi113;
	display u113;
	display pr114;
	display pi114;
	display u114;
	display pr115;
	display pi115;
	display u115;
	display pr116;
	display pi116;
	display u116;
	display pr117;
	display pi117;
	display u117;
	display pr118;
	display pi118;
	display u118;
	display pr119;
	display pi119;
	display u119;
	display pr120;
	display pi120;
	display u120;
	display pr121;
	display pi121;
	display u121;
	display pr122;
	display pi122;
	display u122;
	display pr123;
	display pi123;
	display u123;
	display pr124;
	display pi124;
	display u124;
	display pr125;
	display pi125;
	display u125;
	display pr126;
	display pi126;
	display u126;
	display pr127;
	display pi127;
	display u127;
	display pr128;
	display pi128;
	display u128;
	display pr129;
	display pi129;
	display u129;
	display pr130;
	display pi130;
	display u130;
	display pr131;
	display pi131;
	display u131;
	display pr132;
	display pi132;
	display u132;
	display pr133;
	display pi133;
	display u133;
	display pr134;
	display pi134;
	display u134;
	display pr135;
	display pi135;
	display u135;
	display pr136;
	display pi136;
	display u136;
	display pr137;
	display pi137;
	display u137;
	display pr138;
	display pi138;
	display u138;
	display pr139;
	display pi139;
	display u139;
	display pr140;
	display pi140;
	display u140;
	display pr141;
	display pi141;
	display u141;
	display pr142;
	display pi142;
	display u142;
	display pr143;
	display pi143;
	display u143;
	display pr144;
	display pi144;
	display u144;
	display pr145;
	display pi145;
	display u145;
	display pr146;
	display pi146;
	display u146;
	display pr147;
	display pi147;
	display u147;
	display pr148;
	display pi148;
	display u148;
	display pr149;
	display pi149;
	display u149;
	display pr150;
	display pi150;
	display u150;
	display pr151;
	display pi151;
	display u151;
	display pr152;
	display pi152;
	display u152;
	display pr153;
	display pi153;
	display u153;
	display pr154;
	display pi154;
	display u154;
	display pr155;
	display pi155;
	display u155;
	display pr156;
	display pi156;
	display u156;
	display pr157;
	display pi157;
	display u157;
	display pr158;
	display pi158;
	display u158;
	display pr159;
	display pi159;
	display u159;
	display pr160;
	display pi160;
	display u160;
	display pr161;
	display pi161;
	display u161;
	display pr162;
	display pi162;
	display u162;
	display pr163;
	display pi163;
	display u163;
	display pr164;
	display pi164;
	display u164;
	display pr165;
	display pi165;
	display u165;
	display pr166;
	display pi166;
	display u166;
	display pr167;
	display pi167;
	display u167;
	display pr168;
	display pi168;
	display u168;
	display pr169;
	display pi169;
	display u169;
	display pr170;
	display pi170;
	display u170;
	display pr171;
	display pi171;
	display u171;
	display pr172;
	display pi172;
	display u172;
	display pr173;
	display pi173;
	display u173;
	display pr174;
	display pi174;
	display u174;
	display pr175;
	display pi175;
	display u175;
	display pr176;
	display pi176;
	display u176;
	display pr177;
	display pi177;
	display u177;
	display pr178;
	display pi178;
	display u178;
	display pr179;
	display pi179;
	display u179;
	display pr180;
	display pi180;
	display u180;
	display pr181;
	display pi181;
	display u181;
	display pr182;
	display pi182;
	display u182;
	display pr183;
	display pi183;
	display u183;
	display pr184;
	display pi184;
	display u184;
	display pr185;
	display pi185;
	display u185;
	display pr186;
	display pi186;
	display u186;
	display pr187;
	display pi187;
	display u187;
	display pr188;
	display pi188;
	display u188;
	display pr189;
	display pi189;
	display u189;
	display pr190;
	display pi190;
	display u190;
	display pr191;
	display pi191;
	display u191;
	display pr192;
	display pi192;
	display u192;
	display pr193;
	display pi193;
	display u193;
	display pr194;
	display pi194;
	display u194;
	display pr195;
	display pi195;
	display u195;
	display pr196;
	display pi196;
	display u196;
	display pr197;
	display pi197;
	display u197;
	display pr198;
	display pi198;
	display u198;
	display pr199;
	display pi199;
	display u199;
	display pr200;
	display pi200;
	display u200;
display obj;