(* ::Package:: *) (************************************************************************) (* This file was generated automatically by the Mathematica front end. *) (* It contains Initialization cells from a Notebook file, which *) (* typically will have the same name as this file except ending in *) (* ".nb" instead of ".m". *) (* *) (* This file is intended to be loaded into the Mathematica kernel using *) (* the package loading commands Get or Needs. Doing so is equivalent *) (* to using the Evaluate Initialization Cells menu command in the front *) (* end. *) (* *) (* DO NOT EDIT THIS FILE. This entire file is regenerated *) (* automatically each time the parent Notebook file is saved in the *) (* Mathematica front end. Any changes you make to this file will be *) (* overwritten. *) (************************************************************************) (* ::Title:: *) (*Algorithms*) (* ::Section:: *) (*Introduction*) (* ::Text:: *) (*This file contains the algorithms and examples from the manuscript: "Searching for a best LAD-solution of an overdetermined system of linear equations motivated by searching for a best LAD-hyperplane on the basis of given data" written by K. Sabo, R. Scitovski and I. Vazler, Department of Mathematics, University of Osijek.*) (*The cells containing the algorithms don't need to be manually executed in order to run the examples, Mathematica should offer to run them automatically the first time a cell is run. A more detailed explanation of the examples and their illustrations can be found in the article.*) (* ::Section:: *) (*Algorithms*) (* ::Subsection:: *) (*Subroutines necessary for the algorithms *) (* ::Subsubsection::Closed:: *) (*QR decomposition*) (* ::Text:: *) (*input :*) (* matrix A*) (*output:*) (* orthogonal matrix Q*) (* upper triangular matrix R*) QR[A_]:=Module[{m=Length[A],n=Length[A[[1]]],Q,R=N[A],a,gamma,u,H,i,j}, Q=IdentityMatrix[m]; Do[ a=R[[All,i]]; gamma=Norm[a[[i;;m]]]; u=Table[0.,{m}]; u[[i]]=a[[i]]+Sign[Sign[a[[i]]]+.5]*gamma; Do[u[[j]]=a[[j]],{j,i+1,m}]; u/=Sqrt[2*gamma*(gamma+Abs[a[[i]]])]; H=IdentityMatrix[m]-2*Transpose[{u}].{u}; R=H.R; Q=Q.H; ,{i,Min[{n,m}]}]; Chop[{Q,R}] ] (* ::Subsubsection::Closed:: *) (*QR update - coulumn deletion*) (* ::Text:: *) (*input :*) (* orthogonal matrix Q*) (* upper triangular matrix R*) (* integer k*) (*output:*) (* orthogonal matrix QQ*) (* upper triangular matrix RR*) Izbaci[Q_,R_,k_]:=Module[{m=Length[R],n=Length[R[[1]]],RR=Transpose[Delete[Transpose[R],k]],QQ=Q,a,b,G}, Do[ If[RR[[i,i]]!=0||RR[[i+1,i]]!=0, a=RR[[i,i]]/Sqrt[RR[[i,i]]^2+RR[[i+1,i]]^2]; b=RR[[i+1,i]]/Sqrt[RR[[i,i]]^2+RR[[i+1,i]]^2]; G=IdentityMatrix[m]; G[[i;;i+1,i;;i+1]]={{a,b},{b,-a}}; RR=Chop[G.RR]; QQ=QQ.G; ] ,{i,k,n-1}]; Chop[N[{QQ,RR}]] ] (* ::Subsubsection::Closed:: *) (*QR update - coulumn insertion*) (* ::Text:: *) (*input :*) (* orthogonal matrix Q*) (* upper triangular matrix R*) (* vector vekt*) (*output:*) (* orthogonal matrix QQ*) (* upper triangular matrix RR*) Ubaci[Q_,R_,vekt_]:=Module[{m=Length[R],n=Length[R[[1]]]+1,RR=Transpose[Append[Transpose[R],Transpose[Q].vekt]],QQ=Q,a,b,G}, Do[ If[RR[[i+1,n]]!=0, a=RR[[i,n]]/Sqrt[RR[[i,n]]^2+RR[[i+1,n]]^2]; b=RR[[i+1,n]]/Sqrt[RR[[i,n]]^2+RR[[i+1,n]]^2]; G=IdentityMatrix[m]; G[[i;;i+1,i;;i+1]]={{a,b},{b,-a}}; RR=Chop[G.RR]; QQ=QQ.G; ] ,{i,m-1,n,-1}]; Chop[N[{QQ,RR}]] ] (* ::Subsubsection::Closed:: *) (*Upper triangular matrix inverse*) (* ::Text:: *) (*input :*) (* upper triangular matrix R*) (*output:*) (* upper triangular matrix InvR*) InvGT[R_]:=Module[{InvR,i,j,k,n=Length[R]}, InvR=IdentityMatrix[n]; Do[ Do[ InvR[[i,j]]=(InvR[[i,j]]-Sum[R[[k,j]]*InvR[[i,k]],{k,i,j-1}])/R[[j,j]]; ,{j,i,n}]; ,{i,n,1,-1}]; InvR ] (* ::Subsubsection::Closed:: *) (*QR decomposition w/ pivoting*) (* ::Text:: *) (*input :*) (* matrix A*) (*output:*) (* orthogonal matrix Q*) (* upper triangular matrix R*) (* permutation vector I*) QRnez[A_]:=Module[{m=Length[A],n=Length[A[[1]]],Q,R=N[A],a,gamma,u,H,i,j,I,naj}, I=Range[n]; Q=IdentityMatrix[m]; Do[ naj={i,Norm[R[[i;;m,i]]]}; Do[If[naj[[2]]1, Do[If[naj[[2]]0&&di[[ii]]ToString[i]],{i,Length[I0]}]; uvjeti=Flatten[Join[{X0.ll==h},Table[-1<=ll[[i]]<=1,{i,Length[ll]}]]]; If[Depth[FindInstance[And@@uvjeti,ll]]>2, Return[a]; , I0bezI2IB=Complement[Complement[I0,I2],IB]; While[Length[I0bezI2IB]>0, (*step 6*) ni=I0bezI2IB[[1]]; AppendTo[I2,ni]; IBpot=Prepend[IB,ni]; Bpot=X[[IBpot]]; rje=QRnez2[Transpose[Bpot]]; {Q,R}=Izbaci[rje[[1]],rje[[2]],n+1]; IB=IBpot[[rje[[3,1;;n]]]]; B=X[[IB]]; zB=z[[IB]]; BInv=Q.Transpose[InvGT[R]]; BInvT=Transpose[BInv]; j0=0; I0bezIB=Complement[I0,IB]; Do[ If[Abs[BInvT[[j]].h]>1+Sum[Abs[BInvT[[j]].X[[I0bezIB[[i]]]]],{i,Length[I0bezIB]}], j0=IB[[j]];k=j; Break[]; ]; ,{j,n}]; If[j0!=0, Goto[step3]; , I0bezI2IB=Complement[Complement[I0,I2],IB]; ]; ]; Print["Error!"]; Return[]; ]; ]; (*step 4*) {Q,R}=Izbaci[Q,R,k]; {Q,R}=Ubaci[Q,R,X[[ni]]]; IB=Append[Delete[IB,k],ni]; B=Transpose[Q.R]; BInv=Q.Transpose[InvGT[R]]; BInvT=Transpose[BInv]; I0={}; Do[If[Chop[z[[i]]-a.X[[i]],10.^(-prec)]==0,AppendTo[I0,i]],{i,m}]; IbezI0=Complement[Range[m],I0]; h=Table[0,{n}]+Sum[Sign[z[[IbezI0[[i]]]]-a.X[[IbezI0[[i]]]]]*X[[IbezI0[[i]]]],{i,Length[IbezI0]}]; ]; Print["Maximum iterations reached. Current point a=",a]; Return[a]; ] (* ::Section:: *) (*Examples*) (* ::Subsection::Closed:: *) (*The minimizing function*) F[al_,X_,z_]:=Total[Abs[z-X.al]] (* ::Subsection::Closed:: *) (*Algorithm 1 examples*) (* ::Subsubsection:: *) (*Example 1*) (* ::Input:: *) (*SeedRandom[3];*) (*m=50;*) (*\[CapitalLambda]=RandomInteger[{0,2m},{m,3}]//N*) (*X=\[CapitalLambda][[All,{1,2}]];*) (*z=\[CapitalLambda][[All,3]];*) (* ::Input:: *) (*rj=Alg1[\[CapitalLambda]];*) (*a=Table[ToExpression["a"<>ToString[i]],{i,Length[X[[1]]]}];*) (*Print["NMinimize:",NMinimize[F[a,X,z],a]];*) (*Print["Algoritam 1:",{F[rj,X,z],rj}];*) (* ::Subsubsection:: *) (*Example 2 (Sabo)*) (* ::Input:: *) (*\[CapitalLambda]={{8,4,4},{2,-1,1},{4,3,0},{-2,6,6}};*) (*X=\[CapitalLambda][[All,{1,2}]];*) (*z=\[CapitalLambda][[All,3]];*) (* ::Input:: *) (*rj=Alg1[\[CapitalLambda]];*) (*a=Table[ToExpression["a"<>ToString[i]],{i,Length[X[[1]]]}];*) (*Print["NMinimize:",NMinimize[F[a,X,z],a]];*) (*Print["Algoritam 1:",{F[rj,X,z],rj}];*) (* ::Subsection::Closed:: *) (*Algorithm 2 examples*) (* ::Subsubsection::Closed:: *) (*Example 1*) (* ::Input:: *) (*SeedRandom[3];*) (*m=50;*) (*\[CapitalLambda]=RandomInteger[{0,2m},{m,3}]//N;*) (*X=\[CapitalLambda][[All,{1,2}]]*) (*z=\[CapitalLambda][[All,3]]*) (* ::Input:: *) (*rj=Alg2[X,z][[-1,2]];*) (*a=Table[ToExpression["a"<>ToString[i]],{i,Length[X[[1]]]}];*) (*Print["NMinimize:",NMinimize[F[a,X,z],a]];*) (*Print["Algoritam 2:",{F[rj,X,z],rj}];*) (* ::Subsubsection::Closed:: *) (*Example 2 (Sabo)*) (* ::Input:: *) (*X={{8,4},{2,-1},{4,3},{-2,6}};*) (*z={4,1,0,6};*) (* ::Text:: *) (*In this illustration of the problem the numbers near the edges of the image denote the ordinal number of the equation and the numbers on the intersections are the values of the function which is minimized.*) (*Below the picture is the result of minimization using Mathematica-s function NMinimize.*) (* ::Input:: *) (*m=Length[X];*) (*TT={};*) (*Do[bla=NSolve[X[[{i,j}]].{x1,x2}==z[[{i,j}]],{x1,x2}];*) (*If[Depth[bla]==4&&Length[bla[[1]]]==2,AppendTo[TT,{x1,x2}/.bla[[1]]]];*) (*,{i,m-1},{j,i+1,m}];*) (*TT=Union[Round[TT,0.000001]];*) (*pol={TT[[1]]};*) (*Do[*) (*If[Chop[F[TT[[i]],X,z]-F[pol[[-1]],X,z]]==0,AppendTo[pol,TT[[i]]]];*) (*If[Chop[F[TT[[i]],X,z]-F[pol[[-1]],X,z]]<0,pol={TT[[i]]};];*) (*,{i,Length[TT]}];*) (*slpol=Graphics[{RGBColor[1,1,0],Opacity[.5],Polygon[pol]}];*) (*slpod=Show[{*) (*Table[ContourPlot[X[[i]].{x1,x2}==z[[i]],{x1,-1,2},{x2,-2.5,2},Frame->False,ColorFunction->(Gray&),ContourLabels->({RGBColor[0,0,1],Text[Style[i,Bold],{#1+.1,#2}]}&)],{i,m}],*) (*slpol,*) (*Graphics[Table[Text[Round[F[TT[[i]],X,z],.1],TT[[i]],{0,-1}],{i,Length[TT]}]]},ImageSize->700,AspectRatio->1/2,Axes->True,TicksStyle->Directive[Orange,Opacity[.5]],AxesStyle->Directive[Orange,Opacity[.5]]]*) (*Print["NMinimize:",NMinimize[F[{a1,a2},X,z],{a1,a2}]]*) (* ::Text:: *) (*In the table below we illustrate how the algorithm works starting from four different starting points.*) (*In every iteration we show the current point and the function value in it.*) (* ::Input:: *) (*AAA={Alg2[X,z,{2,3},10],Alg2[X,z,{2,4},10],Alg2[X,z,{3,4},10],Alg2[X,z,{2,1},10]};*) (*TableForm[Table[Table[{{AAA[[i,j,2,1]],AAA[[i,j,2,2]]},F[AAA[[i,j,2]],X,z]},{j,Length[AAA[[i]]]}],{i,Length[AAA]}],TableHeadings->{{"1.","2.","3.","4."},{"Initial","1st iteration","2nd iteration"}},TableDepth->2]*) (*TableForm[*) (*Table[*) (*a=AAA[[k,-1,2]];*) (*I0={};*) (*Do[If[Chop[z[[i]]-X[[i]].a]==0,AppendTo[I0,i]],{i,Length[X]}];*) (*IbezI0=Complement[Range[Length[X]],I0];*) (*ha=Sum[Sign[z[[IbezI0[[i]]]]-X[[IbezI0[[i]]]].a]*X[[IbezI0[[i]]]],{i,Length[IbezI0]}];*) (*ll=Table[ToExpression["l"<>ToString[i]],{i,Length[I0]}];*) (*uvjeti=Flatten[Join[{Transpose[X[[I0]]].ll==ha},Table[-1<=ll[[i]]<=1,{i,Length[ll]}]]];*) (*FindInstance[And@@uvjeti,ll]*) (*,{k,4}],TableHeadings->{{"1.","2.","3.","4."},{"\!\(\*SuperscriptBox[\"\[Lambda]\", \"*\"]\)(\!\(\*SubscriptBox[\"I\", \"0\"]\))"}},TableDepth->2]*) (* ::Text:: *) (*Pictures below are illustrations of the algorithm iterations above.*) (* ::Input:: *) (*slike=GraphicsGrid[Table[aaa=AAA[[2i+j]];*) (*slpod=Show[Table[ContourPlot[X[[i]].{x1,x2}==z[[i]],{x1,-1,2},{x2,-2.5,2},Frame->False,ColorFunction->(Gray&)],{i,m}]];*) (*slF=Graphics[Table[Text[F[TT[[i]],X,z],TT[[i]],{0,-1}],{i,Length[TT]}]];*) (*plave=Graphics[{{RGBColor[0,0,1],Thickness[.003],Arrowheads[.02],Table[Arrow[{aaa[[i,2]],aaa[[i+1,2]]-.015(aaa[[i+1,2]]-aaa[[i,2]])/Norm[aaa[[i+1,2]]-aaa[[i,2]]]}],{i,Length[aaa]-1}]},{PointSize[.015],Point/@aaa[[All,2]]}}];*) (*Show[{slpod,slpol,plave,slF},ImageSize->400,Axes->True,TicksStyle->Directive[Orange,Opacity[.5]],AxesStyle->Directive[Orange,Opacity[.5]],AspectRatio->1/2],{i,0,1},{j,2}]]*) (* ::Subsubsection::Closed:: *) (*Example 4*) (* ::Text:: *) (*This example shows the degenerative case.*) (* ::Input:: *) (*X={{-2,-1},{-1,0},{0,-3},{0,-1},{1,-2},{1,-1},{1,1},{2,-3},{2,-2},{2,0}};*) (*z={-5,-2,-9,-1,-2,0,4,-1,2,2};*) (* ::Text:: *) (*This shows the result of minimization using Mathematica-s function NMinimize.*) (*In the table beneath we illustrate how the algorithm works starting from four different starting points.*) (*In every iteration we show the current point and the function value in it.*) (* ::Input:: *) (*Print["NMinimize:",NMinimize[F[{a1,a2},X,z],{a1,a2}]]*) (*AAA={Alg2[X,z,{9,4},10],Alg2[X,z,{4,10},10],Alg2[X,z,{3,9},10],Alg2[X,z,{3,10},10]};*) (*TableForm[Table[Table[{{AAA[[i,j,2,1]],AAA[[i,j,2,2]]},F[AAA[[i,j,2]],X,z]},{j,Length[AAA[[i]]]}],{i,Length[AAA]}],TableHeadings->{{"1.","2.","3.","4."},{"Initial","1st iteration","2nd iteration"}},TableDepth->2]*) (*TableForm[*) (*Table[*) (*a=AAA[[k,-1,2]];*) (*I0={};*) (*Do[If[Chop[z[[i]]-X[[i]].a]==0,AppendTo[I0,i]],{i,Length[X]}];*) (*IbezI0=Complement[Range[Length[X]],I0];*) (*ha=Sum[Sign[z[[IbezI0[[i]]]]-X[[IbezI0[[i]]]].a]*X[[IbezI0[[i]]]],{i,Length[IbezI0]}];*) (*ll=Table[ToExpression["l"<>ToString[i]],{i,Length[I0]}];*) (*uvjeti=Flatten[Join[{Transpose[X[[I0]]].ll==ha},Table[-1<=ll[[i]]<=1,{i,Length[ll]}]]];*) (*FindInstance[And@@uvjeti,ll]*) (*,{k,4}],TableHeadings->{{"1.","2.","3.","4."},{"\!\(\*SuperscriptBox[\"\[Lambda]\", \"*\"]\)(\!\(\*SubscriptBox[\"I\", \"0\"]\))"}},TableDepth->2]*) (* ::Text:: *) (*In the illustration below the first picture depicts the problem.*) (*The other pictures show the algorithm iterations above (blue arrows) and the optimal paths to the solution (green arrows).*) (* ::Input:: *) (*m=Length[X];*) (*TT={};*) (*Do[bla=NSolve[X[[{i,j}]].{x1,x2}==z[[{i,j}]],{x1,x2}];*) (*If[Depth[bla]==4&&Length[bla[[1]]]==2,AppendTo[TT,{x1,x2}/.bla[[1]]]];*) (*,{i,m-1},{j,i+1,m}];*) (*TT=Union[Round[TT,10.^(-12)]];*) (*pol={TT[[1]]};*) (*Do[*) (*If[Chop[F[TT[[i]],X,z]-F[pol[[-1]],X,z]]==0,AppendTo[pol,TT[[i]]]];*) (*If[Chop[F[TT[[i]],X,z]-F[pol[[-1]],X,z]]<0,pol={TT[[i]]};];*) (*,{i,Length[TT]}];*) (*slpol=Graphics[{RGBColor[1,1,0],Opacity[.5],Polygon[pol]}];*) (*sl11=Show[{slpol,*) (*Table[If[i!=2&&i!=10,*) (*ContourPlot[X[[i]].{x1,x2}==z[[i]],{x1,-.25,5},{x2,-.25,4},Frame->False,ColorFunction->(Gray&),ContourLabels->({RGBColor[0,0,1],Text[Style[i,Bold],{#1+.1,#2}]}&)],*) (*ContourPlot[X[[i]].{x1,x2}==z[[i]],{x1,-.25,5},{x2,-.25,4},Frame->False,ColorFunction->(Gray&),ContourLabels->({RGBColor[0,0,1],Text[Style[i,Bold],{#1+.1,3.9}],Text[Style[i,Bold],{#1+.1,-.2}]}&)]*) (*],{i,m}],*) (*Graphics[Table[Text[Round[F[TT[[i]],X,z],.1],TT[[i]],{0,-1}],{i,Length[TT]}]]},AspectRatio->Automatic,Axes->True,TicksStyle->Directive[Orange,Opacity[.5]],AxesStyle->Directive[Orange,Opacity[.5]]];*) (*zelene={Graphics[{RGBColor[0,1,0],Thickness[.001],Arrowheads[.02],Arrow[{AAA[[1,1,2]],pol[[2]]}],Arrow[{AAA[[1,1,2]],pol[[3]]}]}],*) (*Graphics[{RGBColor[0,1,0],Thickness[.001],Arrowheads[.02],Arrow[{AAA[[2,1,2]],pol[[1]]}],Arrow[{AAA[[2,1,2]],pol[[2]]}]}],*) (*Graphics[{RGBColor[0,1,0],Thickness[.001],Arrowheads[.02],Arrow[{AAA[[3,1,2]],pol[[4]]}],Arrow[{AAA[[3,1,2]],pol[[3]]}]}],*) (*Graphics[{RGBColor[0,1,0],Thickness[.001],Arrowheads[.02],Arrow[{AAA[[4,1,2]],pol[[1]]}],Arrow[{AAA[[4,1,2]],pol[[4]]}]}]};*) (*sldod={ContourPlot[{X[[4]].{x1,x2}==z[[4]],X[[9]].{x1,x2}==z[[9]]},{x1,-.25,4.25},{x2,-.25,3.25},Frame->False,ColorFunction->(Yellow&)],ContourPlot[{X[[4]].{x1,x2}==z[[4]],X[[10]].{x1,x2}==z[[10]]},{x1,-.25,4.25},{x2,-.25,3.25},Frame->False,ColorFunction->(Yellow&)],Graphics[],Graphics[]};*) (*slike=Table[aaa=AAA[[2i+j]];*) (*slpod=Show[Table[ContourPlot[X[[i]].{x1,x2}==z[[i]],{x1,-.25,5},{x2,-.25,4},Frame->False,ColorFunction->(Gray&)],{i,m}]];*) (*slF=Graphics[Table[Text[F[TT[[i]],X,z],TT[[i]],{0,-1}],{i,Length[TT]}]];*) (*plave=Graphics[{{RGBColor[0,0,1],Thickness[.003],Arrowheads[.02],Table[Arrow[{aaa[[i,2]],aaa[[i+1,2]]-.015(aaa[[i+1,2]]-aaa[[i,2]])/Norm[aaa[[i+1,2]]-aaa[[i,2]]]}],{i,Length[aaa]-1}]},{PointSize[.015],Point/@aaa[[All,2]]}}];*) (*Show[{slpol,slpod,zelene[[2*i+j]],plave,slF},ImageSize->400,Axes->True,TicksStyle->Directive[Orange,Opacity[.5]],AxesStyle->Directive[Orange,Opacity[.5]],AspectRatio->Automatic],{i,0,1},{j,2}];*) (*slike[[1]]={sl11,aaa=AAA[[1]];*) (*slpod=Show[Table[ContourPlot[X[[i]].{x1,x2}==z[[i]],{x1,-.25,5},{x2,-.25,4},Frame->False,ColorFunction->(Gray&)],{i,m}]];*) (*slF=Graphics[Table[Text[F[TT[[i]],X,z],TT[[i]],{0,-1}],{i,Length[TT]}]];*) (*plave1=Graphics[{PointSize[.015],Point/@aaa[[All,2]]}];*) (*aaa=AAA[[2]];*) (*plave2=Graphics[{PointSize[.015],Point/@aaa[[All,2]]}];*) (*Show[{slpol,slpod,zelene[[1]],zelene[[2]],plave1,plave2,slF},ImageSize->400,Axes->True,TicksStyle->Directive[Orange,Opacity[.5]],AxesStyle->Directive[Orange,Opacity[.5]],AspectRatio->Automatic]};*) (*GraphicsGrid[slike,ImageSize->800]*) (* ::Text:: *) (*In the table below *) (*\!\(\*SuperscriptBox["\[Lambda]", "*"]\)(Subscript[I, 0]) for each of the characteristic solutions are calculated.*) (* ::Input:: *) (*TableForm[*) (*Table[{a=pol[[k]];*) (*a,*) (*I0={};*) (*Do[If[Chop[z[[i]]-X[[i]].a]==0,AppendTo[I0,i]],{i,Length[X]}];*) (*IbezI0=Complement[Range[Length[X]],I0];*) (*ha=Sum[Sign[z[[IbezI0[[i]]]]-X[[IbezI0[[i]]]].a]*X[[IbezI0[[i]]]],{i,Length[IbezI0]}];*) (*ll=Table[ToExpression["l"<>ToString[i]],{i,Length[I0]}];*) (*uvjeti=Flatten[Join[{Transpose[X[[I0]]].ll==ha},Table[-1<=ll[[i]]<=1,{i,Length[ll]}]]];*) (*Print["Intersecting lines:",I0,"\nConditions:",uvjeti];*) (*FindInstance[And@@uvjeti,ll]}*) (*,{k,4}],TableHeadings->{{"1.","2.","3.","4."},{"a","\!\(\*SuperscriptBox[\"\[Lambda]\", \"*\"]\)(\!\(\*SubscriptBox[\"I\", \"0\"]\))"}},TableDepth->2]*) (* ::Subsection:: *) (*Algorithm 3 examples*) (* ::Subsubsection::Closed:: *) (*Example 1*) (* ::Input:: *) (*X={{1,1,1},{3,3,1},{0,0,3},{3,3,3},{2,3,1}};*) (*z={1,2,0,4,2};*) (* ::Text:: *) (*Below is the comparison of Mathematica function NMinimize and our Algorithm 3.*) (* ::Input:: *) (*rj=Alg3[X,z];*) (*a=Table[ToExpression["a"<>ToString[i]],{i,Length[X[[1]]]}];*) (*Print["NMinimize:",NMinimize[F[a,X,z],a]];*) (*Print["Algoritam 3:",{F[rj,X,z],rj}];*) (* ::Subsubsection::Closed:: *) (*Example 2 (Sabo)*) (* ::Input:: *) (*X={{8,4},{2,-1},{4,3},{-2,6}};*) (*z={4,1,0,6};*) (* ::Text:: *) (*Below is the comparison of Mathematica function NMinimize and our Algorithm 3.*) (* ::Input:: *) (*rj=Alg3[X,z];*) (*a=Table[ToExpression["a"<>ToString[i]],{i,Length[X[[1]]]}];*) (*Print["NMinimize:",NMinimize[F[a,X,z],a]];*) (*Print["Algoritam 3:",{F[rj,X,z],rj}];*) (* ::Subsubsection::Closed:: *) (*Example 3 {Coleman}*) (* ::Text:: *) (*This example illustrates a LAD - problem of finding a polynomial approximation of the function.*) (* ::Input:: *) (*f[x_]:=Exp[x]+If[1.2<= x<= 1.4,5,0]*) (*slf=Plot[f[x],{x,1,2},PlotRange->{2,9.5},PlotStyle->{RGBColor[0,0,1],Opacity[.3]}];*) (*nnn=10;m=50; sig=0.5;SeedRandom[25]*) (*x=y=Table[1+i /m,{i,m}]; *) (*error=RandomReal[NormalDistribution[0,sig^2],m];*) (*Do[y[[i]]=f[x[[i]]]+error[[i]],{i,m}];*) (*tab=Table[{x[[i]],f[x[[i]]],error[[i]],y[[i]]},{i,m}];*) (*(*Transpose[tab]//TableForm *)*) (*pod=Table[{x[[i]],y[[i]]},{i,m}];*) (*slpod=ListPlot[pod,PlotRange->{2,9.5},PlotStyle->{AbsolutePointSize[4]}];*) (*Show[slpod,slf,ImageSize->300,AspectRatio->.5,Axes->True,TicksStyle->Directive[Orange,Opacity[.5]],AxesStyle->Directive[Orange,Opacity[.5]]]*) (**) (*(* Definiranje sustava *)*) (*X=Table[Table[x[[i]]^(j-1)//N,{j,nnn}],{i,m}];*) (*(* Print[Transpose[X]//TableForm]*)*) (*z=y;*) (**) (*rj=Alg3[X,z,,200,6];*) (*a=Table[ToExpression["a"<>ToString[i]],{i,Length[X[[1]]]}];*) (*Print["NMinimize:",NMinimize[F[a,X,z],a]];*) (*Print["Algoritam 3:",{F[rj,X,z],rj}];*) (* ::Subsubsection::Closed:: *) (*Example 4*) (* ::Text:: *) (*This is the same Example 4 from the Algorithm 2 examples.*) (* ::Input:: *) (*Graphics[{{RGBColor[1, 1, 0], Opacity[0.5], Polygon[{{1.666666666667, 1.666666666667}, {1.75, 1.5}, {2., 1.666666666667}, {2., 2.}}]}, GraphicsComplex[CompressedData["*) (*1:eJxdVG1oU1cYvjf31iRtI37tj0LVEllxiANd8Yf6vHOUGhGZ/SHbWnR+VGQ0*) (*s4KCK7gVoWX4p6wKCgoKVtkYKoLEOrXZJjjaDiwYtBjFZmtrqjX9ykfbJI33*) (*fc/xjwfCgZtz3vf5es/KfYdr6l2GYex3fryrNQC1L6CBUDAaCv6HX/yBdn/3*) (*fOrZZfs+PvY/zkdO/NBQNZ+WXa+ONS0ZxFRFXdvfKR/d4+PRQZQ7xwPtPlp7*) (*svX0H2uHEO/aGCo/5CNVd1jX99Hpa56lO14Po/K37GT/qVI6x30Cr3CL69SW*) (*0qdctiKO+M7OspY3JZTuP3Up1hFHi1O17VoJNXK70AiquZ2/hL5meHte4zHD*) (*6yvW/UZ1v2K679zyLB1F+aENTkcvbZd7o3jGZd56aUvx08uNm0bxsmmJw8xL*) (*K2oXrr8QeYvDNZmhmx95qUJ4JdAgfT26/piq3+yhVqbbNgahu8FDawT/OB47*) (*VYufukk0bh7HXMFZP7kpKXzGUeqgLqt163oTGq+briZ6939ycgLfMIwT8zTf*) (*CeVHYB5VruM1iZ5/eRXRV8J/Er+ynJeK6LnwmoQcX1REXkftro1TYDafP7B1*) (*v6TuZ9Nd0SeJqr5Gh7FF20SfJCQGUYsg+iTxV4oNtygoB5MQOzwWdYpvKWzl*) (*ex0uWi76pRDrYCIu4rRcr06D3dzZaZLgak9DaAVNkrY30xAYO0yNL4OCLIMk*) (*Tm0Z/Mj2Nhok8Loz+J4LVBr0pRSexm6Wva6AdRc4qNM49w8Ht6B5zuh9DsMs*) (*Q2AG6+VcXn+f1XtO35uFwpXVfWaVnu1ZiL11s1C8stj24qxzI4sG0WUGon8w*) (*Cznu132bszA1jhzbFMtC+TKN+oO8cvhMfM1oHXOqbygNGZeuHK5ILtIaZ17r*) (*k8K3ZSx4Hg4pB0kS4kc0D2kfSOLIJjZwDiJ/2RQWCfE5CGzflM5JQfs9ASmz*) (*q4DjQuh9Lg2S/BbGoPJi0GoxZkzzNEjdT2Cz5MWgxfJOJHCGc/KzQQ5458Ib*) (*Ve9PgxT+Eag8GSR2VMXxgHO2xaQRmf8h3JZ8mSS6+IdUTutNkuN9g1rf97kZ*) (*hMgRMckn8xVTOqZNlZvAgJqDL1zEU5jojWofXdTK8rREsUoIuTSfJzoHLlK+*) (*RbQeFs0I/keQZ2yFRRFVUOV6jUVqXnsgY7rdou/E54dQebbo6B5Gek/rYdFL*) (*qXdD+dZtUb8A+B1CJ2/RMuGzFzWSd5vkuWq6FRa4B2w9t3fCKtc2uUWAu2Gl*) (*j022ahSWz1GbLsqA9qr/jSKtX1/4w/f/w/PvAEJ7rdU=*) (*"], {{{}, {}, {Hue[0.67, 0.6, 0.6], Tooltip[{{GrayLevel[0.5], Line[{1, 2}]}, {GrayLevel[0.5], Line[{2, 3}]}, {GrayLevel[0.5], Line[{3, 4}]}, {GrayLevel[0.5], Line[{4, 5}]}, {GrayLevel[0.5], Line[{5, 6}]}, {GrayLevel[0.5], Line[{6, 7}]}, {GrayLevel[0.5], Line[{7, 8}]}, {GrayLevel[0.5], Line[{8, 9}]}, {GrayLevel[0.5], Line[{9, 10}]}, {GrayLevel[0.5], Line[{10, 11}]}, {GrayLevel[0.5], Line[{11, 12}]}, {GrayLevel[0.5], Line[{12, 13}]}, {GrayLevel[0.5], Line[{13, 14}]}, {GrayLevel[0.5], Line[{14, 15}]}, {GrayLevel[0.5], Line[{15, 16}]}, {GrayLevel[0.5], Line[{16, 17}]}, {GrayLevel[0.5], Line[{17, 18}]}, {GrayLevel[0.5], Line[{18, 19}]}, {GrayLevel[0.5], Line[{19, 20}]}, {GrayLevel[0.5], Line[{20, 21}]}, {GrayLevel[0.5], Line[{21, 22}]}, {GrayLevel[0.5], Line[{22, 23}]}, {GrayLevel[0.5], Line[{23, 24}]}, {GrayLevel[0.5], Line[{24, 25}]}, {GrayLevel[0.5], Line[{25, 26}]}, {GrayLevel[0.5], Line[{26, 27}]}, {GrayLevel[0.5], Line[{27, 28}]}, {GrayLevel[0.5], Line[{28, 29}]}, {GrayLevel[0.5], Line[{29, 30}]}, {GrayLevel[0.5], Line[{30, 31}]}, {GrayLevel[0.5], Line[{31, 32}]}, {GrayLevel[0.5], Line[{32, 33}]}, {GrayLevel[0.5], Line[{33, 34}]}, {GrayLevel[0.5], Line[{34, 35}]}, {GrayLevel[0.5], Line[{35, 36}]}, {GrayLevel[0.5], Line[{36, 37}]}, {GrayLevel[0.5], Line[{37, 38}]}, {GrayLevel[0.5], Line[{38, 39}]}, {GrayLevel[0.5], Line[{39, 40}]}, {GrayLevel[0.5], Line[{40, 41}]}, {GrayLevel[0.5], Line[{41, 42}]}, {GrayLevel[0.5], Line[{42, 43}]}, {GrayLevel[0.5], Line[{43, 44}]}, {GrayLevel[0.5], Line[{44, 45}]}, {GrayLevel[0.5], Line[{45, 46}]}, {GrayLevel[0.5], Line[{46, 47}]}, {GrayLevel[0.5], Line[{47, 48}]}, {GrayLevel[0.5], Line[{48, 49}]}, {GrayLevel[0.5], Line[{49, 50}]}, {GrayLevel[0.5], Line[{50, 51}]}, {GrayLevel[0.5], Line[{51, 52}]}, {GrayLevel[0.5], Line[{52, 53}]}, {GrayLevel[0.5], Line[{53, 54}]}, {GrayLevel[0.5], Line[{54, 55}]}, {GrayLevel[0.5], Line[{55, 56}]}, {GrayLevel[0.5], Line[{56, 57}]}, {GrayLevel[0.5], Line[{57, 58}]}, {GrayLevel[0.5], Line[{58, 59}]}, {GrayLevel[0.5], Line[{59, 60}]}, {GrayLevel[0.5], Line[{60, 61}]}, {GrayLevel[0.5], Line[{61, 62}]}, {GrayLevel[0.5], Line[{62, 63}]}, {GrayLevel[0.5], Line[{63, 64}]}, {GrayLevel[0.5], Line[{64, 65}]}, {GrayLevel[0.5], Line[{65, 66}]}, {GrayLevel[0.5], Line[{66, 67}]}, {GrayLevel[0.5], Line[{67, 68}]}, {GrayLevel[0.5], Line[{68, 69}]}, {GrayLevel[0.5], Line[{69, 70}]}, {GrayLevel[0.5], Line[{70, 71}]}, {GrayLevel[0.5], Line[{71, 72}]}, {GrayLevel[0.5], Line[{72, 73}]}, {GrayLevel[0.5], Line[{73, 74}]}, {GrayLevel[0.5], Line[{74, 75}]}, {GrayLevel[0.5], Line[{75, 76}]}, {GrayLevel[0.5], Line[{76, 77}]}, {GrayLevel[0.5], Line[{77, 78}]}, {GrayLevel[0.5], Line[{78, 79}]}, {GrayLevel[0.5], Line[{79, 80}]}, {GrayLevel[0.5], Line[{80, 81}]}, {GrayLevel[0.5], Line[{81, 82}]}, {GrayLevel[0.5], Line[{82, 83}]}, {GrayLevel[0.5], Line[{83, 84}]}, {GrayLevel[0.5], Line[{84, 85}]}, {GrayLevel[0.5], Line[{85, 86}]}, {GrayLevel[0.5], Line[{86, 87}]}, {GrayLevel[0.5], Line[{87, 88}]}, {GrayLevel[0.5], Line[{88, 89}]}, {GrayLevel[0.5], Line[{89, 90}]}, {GrayLevel[0.5], Line[{90, 91}]}, {GrayLevel[0.5], Line[{91, 92}]}, {GrayLevel[0.5], Line[{92, 93}]}, {GrayLevel[0.5], Line[{93, 94}]}, {GrayLevel[0.5], Line[{94, 95}]}, {GrayLevel[0.5], Line[{95, 96}]}, {GrayLevel[0.5], Line[{96, 97}]}}, (-2) $CellContext`x1 - $CellContext`x2 == -5]}}, {{RGBColor[0, 0, 1], Inset[Style[1, FontWeight -> Bold], {0.6758928571428576, 3.848214285714285}]}, {RGBColor[0, 0, 1], Inset[Style[1, FontWeight -> Bold], {2.706991525423729, -0.21398305084745758`}]}}}, VertexColors -> {0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5}], GraphicsComplex[{{2., -0.25}, {2., 0.05357142857142855}, {2., 0.3571428571428571}, {2., 0.6607142857142856}, {2., 0.9642857142857142}, {2., 1.2678571428571428`}, {2., 1.5714285714285714`}, {2., 1.875}, {2., 2.1785714285714284`}, {2., 2.4821428571428568`}, {2., 2.785714285714285}, {2., 3.0892857142857135`}, {2., 3.392857142857142}, {2., 3.6964285714285703`}, {2., 4.}, {2., 0.05357142857142855}, {2., 0.3571428571428571}}, {{{}, {}, {Hue[0.67, 0.6, 0.6], Tooltip[{{GrayLevel[0.5], Line[{1, 2}]}, {GrayLevel[0.5], Line[{2, 3}]}, {GrayLevel[0.5], Line[{3, 4}]}, {GrayLevel[0.5], Line[{4, 5}]}, {GrayLevel[0.5], Line[{5, 6}]}, {GrayLevel[0.5], Line[{6, 7}]}, {GrayLevel[0.5], Line[{7, 8}]}, {GrayLevel[0.5], Line[{8, 9}]}, {GrayLevel[0.5], Line[{9, 10}]}, {GrayLevel[0.5], Line[{10, 11}]}, {GrayLevel[0.5], Line[{11, 12}]}, {GrayLevel[0.5], Line[{12, 13}]}, {GrayLevel[0.5], Line[{13, 14}]}, {GrayLevel[0.5], Line[{14, 15}]}}, -$CellContext`x1 == -2]}}, {{RGBColor[0, 0, 1], Inset[Style[2, FontWeight -> Bold], {2.1, 3.9}], Inset[Style[2, FontWeight -> Bold], {2.1, -0.2}]}, {RGBColor[0, 0, 1], Inset[Style[2, FontWeight -> Bold], {2.1, 3.9}], Inset[Style[2, FontWeight -> Bold], {2.1, -0.2}]}}}, VertexColors -> {0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5}], GraphicsComplex[CompressedData["*) (*1:eJxd1L1Kw1AABeBr+mNrY5M0rhZ9Ax9AyH0DwcW5grNudvM6uekbODv2DawI*) (*gos6OAguvS2IgmhaLBarwdyG3nDOgXA58HESaJP13f3tPUcIsZNe5szycJmd*) (*Ffm72kxzY3uxr9Nc2Z6dHdvlmvGn0bzHh+00FxH6boR7t7Zn97sj/2i7M4zT*) (*PNme1L00z+R7th+b27f7th+dmAzIv5B/Jf9G/p38B/lP8kPb5abJyPbutcmI*) (*/Bf5Mfkx+W/yE/IT8j/kp+Sn5P/IJ+QT8kLO++zvcJ73nvm5dd5nh1pA33TQ*) (*txz0okD7BfS6gF4Vab9E+yX0okz7ZfS6jF4t0n4FfauCXlRpv4peV9GrJdqv*) (*0X4NvXBp30WvXfRqmfbrtF9HLzza99BrD730bffN63qWd/M2xx0fvQjQbwXo*) (*DwL0XfJeA/1GA71qoL8nH5OXIT1/iL4TohcrtvP3k/0/eyEG4Q==*) (*"], {{{}, {}, {Hue[0.67, 0.6, 0.6], Tooltip[{{GrayLevel[0.5], Line[{1, 2}]}, {GrayLevel[0.5], Line[{2, 3}]}, {GrayLevel[0.5], Line[{3, 4}]}, {GrayLevel[0.5], Line[{4, 5}]}, {GrayLevel[0.5], Line[{5, 6}]}, {GrayLevel[0.5], Line[{6, 7}]}, {GrayLevel[0.5], Line[{7, 8}]}, {GrayLevel[0.5], Line[{8, 9}]}, {GrayLevel[0.5], Line[{9, 10}]}, {GrayLevel[0.5], Line[{10, 11}]}, {GrayLevel[0.5], Line[{11, 12}]}, {GrayLevel[0.5], Line[{12, 13}]}, {GrayLevel[0.5], Line[{13, 14}]}, {GrayLevel[0.5], Line[{14, 15}]}, {GrayLevel[0.5], Line[{15, 16}]}, {GrayLevel[0.5], Line[{16, 17}]}, {GrayLevel[0.5], Line[{17, 18}]}, {GrayLevel[0.5], Line[{18, 19}]}, {GrayLevel[0.5], Line[{19, 20}]}, {GrayLevel[0.5], Line[{20, 21}]}, {GrayLevel[0.5], Line[{21, 22}]}, {GrayLevel[0.5], Line[{22, 23}]}, {GrayLevel[0.5], Line[{23, 24}]}, {GrayLevel[0.5], Line[{24, 25}]}, {GrayLevel[0.5], Line[{25, 26}]}, {GrayLevel[0.5], Line[{26, 27}]}, {GrayLevel[0.5], Line[{27, 28}]}, {GrayLevel[0.5], Line[{28, 29}]}, {GrayLevel[0.5], Line[{29, 30}]}, {GrayLevel[0.5], Line[{30, 31}]}, {GrayLevel[0.5], Line[{31, 32}]}, {GrayLevel[0.5], Line[{32, 33}]}, {GrayLevel[0.5], Line[{33, 34}]}, {GrayLevel[0.5], Line[{34, 35}]}, {GrayLevel[0.5], Line[{35, 36}]}, {GrayLevel[0.5], Line[{36, 37}]}, {GrayLevel[0.5], Line[{37, 38}]}, {GrayLevel[0.5], Line[{38, 39}]}, {GrayLevel[0.5], Line[{39, 40}]}, {GrayLevel[0.5], Line[{40, 41}]}, {GrayLevel[0.5], Line[{41, 42}]}, {GrayLevel[0.5], Line[{42, 43}]}, {GrayLevel[0.5], Line[{43, 44}]}, {GrayLevel[0.5], Line[{44, 45}]}, {GrayLevel[0.5], Line[{45, 46}]}, {GrayLevel[0.5], Line[{46, 47}]}, {GrayLevel[0.5], Line[{47, 48}]}, {GrayLevel[0.5], Line[{48, 49}]}, {GrayLevel[0.5], Line[{49, 50}]}, {GrayLevel[0.5], Line[{50, 51}]}, {GrayLevel[0.5], Line[{51, 52}]}, {GrayLevel[0.5], Line[{52, 53}]}, {GrayLevel[0.5], Line[{53, 54}]}, {GrayLevel[0.5], Line[{54, 55}]}, {GrayLevel[0.5], Line[{55, 56}]}, {GrayLevel[0.5], Line[{56, 57}]}, {GrayLevel[0.5], Line[{57, 58}]}, {GrayLevel[0.5], Line[{58, 59}]}, {GrayLevel[0.5], Line[{59, 60}]}, {GrayLevel[0.5], Line[{60, 61}]}, {GrayLevel[0.5], Line[{61, 62}]}, {GrayLevel[0.5], Line[{62, 63}]}, {GrayLevel[0.5], Line[{63, 64}]}, {GrayLevel[0.5], Line[{64, 65}]}, {GrayLevel[0.5], Line[{65, 66}]}, {GrayLevel[0.5], Line[{66, 67}]}, {GrayLevel[0.5], Line[{67, 68}]}, {GrayLevel[0.5], Line[{68, 69}]}, {GrayLevel[0.5], Line[{69, 70}]}, {GrayLevel[0.5], Line[{70, 71}]}, {GrayLevel[0.5], Line[{71, 72}]}, {GrayLevel[0.5], Line[{72, 73}]}, {GrayLevel[0.5], Line[{73, 74}]}, {GrayLevel[0.5], Line[{74, 75}]}, {GrayLevel[0.5], Line[{75, 76}]}, {GrayLevel[0.5], Line[{76, 77}]}, {GrayLevel[0.5], Line[{77, 78}]}, {GrayLevel[0.5], Line[{78, 79}]}, {GrayLevel[0.5], Line[{79, 80}]}, {GrayLevel[0.5], Line[{80, 81}]}, {GrayLevel[0.5], Line[{81, 82}]}, {GrayLevel[0.5], Line[{82, 83}]}, {GrayLevel[0.5], Line[{83, 84}]}, {GrayLevel[0.5], Line[{84, 85}]}}, (-3) $CellContext`x2 == -9]}}, {{RGBColor[0, 0, 1], Inset[Style[3, FontWeight -> Bold], {-0.07279411764705787, 3.}]}, {RGBColor[0, 0, 1], Inset[Style[3, FontWeight -> Bold], {5.022794117647058, 3.}]}}}, VertexColors -> {0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5}], GraphicsComplex[CompressedData["*) (*1:eJxd1LFKw1AYBeCkjY2NjSapL+CQwaHPkMy61UG4boKDiPgKPkPBvkHcCg7F*) (*xcEuLk5aRUEQFERBLLRYLBZFvTf4h3MOhMuBj0ObJl3a3GtulSzLWv+79Jln*) (*Mc3PYXLV16lLb+jsFt0caR19RL4Tod+P0Dcj9IMQfS9E3wrRN8hb5PsBff4A*) (*fRqgj4q+oXQWaL/omck8+jufvq+PPvbRqxrt19Bnc7Tv0b6HPvbQqyrtV9Fn*) (*s7Tv0r6LPnbRqwrtV9BnM7Tv0L6DPnbQqzLtl9FnJdq3ad9GH9volUX7RV9Z*) (*1flJ/vtxV+c7QV/03H+Rn5Kfkv8kPyE/If9Bfkx+TP6d/Ij8iPxQ+s62zoD8*) (*m/SDts4r+Rfyz+SfyD+Sf5C+bF7Ye+m35n2+IX9N/lJ6x/zgZ+R70g/N/ehK*) (*b5v91in6I+nK3J8T6Wvm+TknfyGd/z/Z/wLp/VWt*) (*"], {{{}, {}, {Hue[0.67, 0.6, 0.6], Tooltip[{{GrayLevel[0.5], Line[{1, 2}]}, {GrayLevel[0.5], Line[{2, 3}]}, {GrayLevel[0.5], Line[{3, 4}]}, {GrayLevel[0.5], Line[{4, 5}]}, {GrayLevel[0.5], Line[{5, 6}]}, {GrayLevel[0.5], Line[{6, 7}]}, {GrayLevel[0.5], Line[{7, 8}]}, {GrayLevel[0.5], Line[{8, 9}]}, {GrayLevel[0.5], Line[{9, 10}]}, {GrayLevel[0.5], Line[{10, 11}]}, {GrayLevel[0.5], Line[{11, 12}]}, {GrayLevel[0.5], Line[{12, 13}]}, {GrayLevel[0.5], Line[{13, 14}]}, {GrayLevel[0.5], Line[{14, 15}]}, {GrayLevel[0.5], Line[{15, 16}]}, {GrayLevel[0.5], Line[{16, 17}]}, {GrayLevel[0.5], Line[{17, 18}]}, {GrayLevel[0.5], Line[{18, 19}]}, {GrayLevel[0.5], Line[{19, 20}]}, {GrayLevel[0.5], Line[{20, 21}]}, {GrayLevel[0.5], Line[{21, 22}]}, {GrayLevel[0.5], Line[{22, 23}]}, {GrayLevel[0.5], Line[{23, 24}]}, {GrayLevel[0.5], Line[{24, 25}]}, {GrayLevel[0.5], Line[{25, 26}]}, {GrayLevel[0.5], Line[{26, 27}]}, {GrayLevel[0.5], Line[{27, 28}]}, {GrayLevel[0.5], Line[{28, 29}]}, {GrayLevel[0.5], Line[{29, 30}]}, {GrayLevel[0.5], Line[{30, 31}]}, {GrayLevel[0.5], Line[{31, 32}]}, {GrayLevel[0.5], Line[{32, 33}]}, {GrayLevel[0.5], Line[{33, 34}]}, {GrayLevel[0.5], Line[{34, 35}]}, {GrayLevel[0.5], Line[{35, 36}]}, {GrayLevel[0.5], Line[{36, 37}]}, {GrayLevel[0.5], Line[{37, 38}]}, {GrayLevel[0.5], Line[{38, 39}]}, {GrayLevel[0.5], Line[{39, 40}]}, {GrayLevel[0.5], Line[{40, 41}]}, {GrayLevel[0.5], Line[{41, 42}]}, {GrayLevel[0.5], Line[{42, 43}]}, {GrayLevel[0.5], Line[{43, 44}]}, {GrayLevel[0.5], Line[{44, 45}]}, {GrayLevel[0.5], Line[{45, 46}]}, {GrayLevel[0.5], Line[{46, 47}]}, {GrayLevel[0.5], Line[{47, 48}]}, {GrayLevel[0.5], Line[{48, 49}]}, {GrayLevel[0.5], Line[{49, 50}]}, {GrayLevel[0.5], Line[{50, 51}]}, {GrayLevel[0.5], Line[{51, 52}]}, {GrayLevel[0.5], Line[{52, 53}]}, {GrayLevel[0.5], Line[{53, 54}]}, {GrayLevel[0.5], Line[{54, 55}]}, {GrayLevel[0.5], Line[{55, 56}]}, {GrayLevel[0.5], Line[{56, 57}]}, {GrayLevel[0.5], Line[{57, 58}]}, {GrayLevel[0.5], Line[{58, 59}]}, {GrayLevel[0.5], Line[{59, 60}]}, {GrayLevel[0.5], Line[{60, 61}]}, {GrayLevel[0.5], Line[{61, 62}]}, {GrayLevel[0.5], Line[{62, 63}]}, {GrayLevel[0.5], Line[{63, 64}]}, {GrayLevel[0.5], Line[{64, 65}]}, {GrayLevel[0.5], Line[{65, 66}]}, {GrayLevel[0.5], Line[{66, 67}]}, {GrayLevel[0.5], Line[{67, 68}]}, {GrayLevel[0.5], Line[{68, 69}]}, {GrayLevel[0.5], Line[{69, 70}]}, {GrayLevel[0.5], Line[{70, 71}]}, {GrayLevel[0.5], Line[{71, 72}]}, {GrayLevel[0.5], Line[{72, 73}]}, {GrayLevel[0.5], Line[{73, 74}]}, {GrayLevel[0.5], Line[{74, 75}]}, {GrayLevel[0.5], Line[{75, 76}]}, {GrayLevel[0.5], Line[{76, 77}]}, {GrayLevel[0.5], Line[{77, 78}]}, {GrayLevel[0.5], Line[{78, 79}]}, {GrayLevel[0.5], Line[{79, 80}]}, {GrayLevel[0.5], Line[{80, 81}]}, {GrayLevel[0.5], Line[{81, 82}]}, {GrayLevel[0.5], Line[{82, 83}]}, {GrayLevel[0.5], Line[{83, 84}]}, {GrayLevel[0.5], Line[{84, 85}]}}, -$CellContext`x2 == -1]}}, {{RGBColor[0, 0, 1], Inset[Style[4, FontWeight -> Bold], {5.055882352941176, 1.}]}, {RGBColor[0, 0, 1], Inset[Style[4, FontWeight -> Bold], {-0.10588235294117626`, 1.}]}}}, VertexColors -> {0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5}], GraphicsComplex[CompressedData["*) (*1:eJxdVV1Ik2EU/vbt38/Ztk/owi4qJv1cSBBlUHJOF5HzQqluohBqUDe5DJKg*) (*i2p1UUhd6LILL5SKwiIDtbLywgkW/VCQEWFqVGQ/1Gxtus2cbn3veY/fhS+b*) (*z9553u8853nOebcq1Lj7kKooSoPxFijXm5jEOGwIRgPB6KvYWH94vD8ch5rl*) (*9U0jO57FCJZPwdWACOiL0b/H/8C3noznyKeeWKptZZ079Qfkc3r4eQm41twx*) (*2txxK3ZP39qrb03A2XaxbsNhWn85fogxCZXlz/ec2/YUjD/GpyQcEGn6X4PM*) (*lwQnJXoL2zv9k+u2p/jcO4mRFBw3WNY3fYCKgZ/XL1ZMQ+pJbfxJ7QR0jgoi*) (*0xz/mXEGGvfmzj9o/Awr3g8enCqbgWpRfuALEARngOiUT4Kb8qb53HfOl4bx*) (*yu5Tw2M/ICtUyKThEfH8BUeJdwZau2wna1p+g0HWYJyBRZ0lZuGCkGd0CkJr*) (*ThivLOwkoRPwVpQxkjV1kflm4fTwmJExCTGD7fvBWeaZArIt8A+onL3TIFhV*) (*dv8z65Q4BxtppYFsaJ+DXXEhUAbWivShnMlL5stBiSjLY/B8IJTKsT5ZoLTR*) (*HAiVB37OQhPpPm/mkbgARHd4DkKC7sEF/n4xLg9kY+c81JHAefOczJ8HQzxD*) (*wTxUkREFIBs78kC0awscr+AirqfGUPCKaMdSBSeocRQkPzYrKPvYgpeoQfhc*) (*xCJxSEEbJVSR2k+1oOw7FcnOsAVln6soefE5xWo+p4/63IofiKcFP9F8WPGY*) (*MCbB8RGbRFTxFS0bynlQsYWEtSPpfEPFxySwHYn2R5XzOcy81F5dDpTzYsWT*) (*1EhOrDJY9F5e5OXk+qxI9NucSPRfWrGBCnNynVZZ9zkXyvlinorb5C3nxY1b*) (*RIM9tEmegSIk2n7mWV+Ecv7sfK6I67UjCF66hlfJSDsGBY8aDcnem3a8L3js*) (*15Dun7t8Xik266ZxqyzGMppXB66mPizGViLgwADNqwepraKsU8TD9Tukvi88*) (*kvekQ8b5S1DOtxPnC2KVYF7AGadEY1/gveSxzNTVwXsr7+NEbBnSrZjhePQy*) (*upDGIOxldLH+XkYXfqV7wsvo4nw+0wcap30+RjceonvDx8g+DXH8kBtpLEd8*) (*jG6p65SP0Y3PhFzVfkb2KeI3fRuke8bPWIQ7yW+dUUO6FjbpjBrXqTNq+FyM*) (*aYPOqOEdupd0Rk3qG9MZ+bxSavq+9Pdoafx/YBp+IA==*) (*"], {{{}, {}, {Hue[0.67, 0.6, 0.6], Tooltip[{{GrayLevel[0.5], Line[{1, 2}]}, {GrayLevel[0.5], Line[{2, 3}]}, {GrayLevel[0.5], Line[{3, 4}]}, {GrayLevel[0.5], Line[{4, 5}]}, {GrayLevel[0.5], Line[{5, 6}]}, {GrayLevel[0.5], Line[{6, 7}]}, {GrayLevel[0.5], Line[{7, 8}]}, {GrayLevel[0.5], Line[{8, 9}]}, {GrayLevel[0.5], Line[{9, 10}]}, {GrayLevel[0.5], Line[{10, 11}]}, {GrayLevel[0.5], Line[{11, 12}]}, {GrayLevel[0.5], Line[{12, 13}]}, {GrayLevel[0.5], Line[{13, 14}]}, {GrayLevel[0.5], Line[{14, 15}]}, {GrayLevel[0.5], Line[{15, 16}]}, {GrayLevel[0.5], Line[{16, 17}]}, {GrayLevel[0.5], Line[{17, 18}]}, {GrayLevel[0.5], Line[{18, 19}]}, {GrayLevel[0.5], Line[{19, 20}]}, {GrayLevel[0.5], Line[{20, 21}]}, {GrayLevel[0.5], Line[{21, 22}]}, {GrayLevel[0.5], Line[{22, 23}]}, {GrayLevel[0.5], Line[{23, 24}]}, {GrayLevel[0.5], Line[{24, 25}]}, {GrayLevel[0.5], Line[{25, 26}]}, {GrayLevel[0.5], Line[{26, 27}]}, {GrayLevel[0.5], Line[{27, 28}]}, {GrayLevel[0.5], Line[{28, 29}]}, {GrayLevel[0.5], Line[{29, 30}]}, {GrayLevel[0.5], Line[{30, 31}]}, {GrayLevel[0.5], Line[{31, 32}]}, {GrayLevel[0.5], Line[{32, 33}]}, {GrayLevel[0.5], Line[{33, 34}]}, {GrayLevel[0.5], Line[{34, 35}]}, {GrayLevel[0.5], Line[{35, 36}]}, {GrayLevel[0.5], Line[{36, 37}]}, {GrayLevel[0.5], Line[{37, 38}]}, {GrayLevel[0.5], Line[{38, 39}]}, {GrayLevel[0.5], Line[{39, 40}]}, {GrayLevel[0.5], Line[{40, 41}]}, {GrayLevel[0.5], Line[{41, 42}]}, {GrayLevel[0.5], Line[{42, 43}]}, {GrayLevel[0.5], Line[{43, 44}]}, {GrayLevel[0.5], Line[{44, 45}]}, {GrayLevel[0.5], Line[{45, 46}]}, {GrayLevel[0.5], Line[{46, 47}]}, {GrayLevel[0.5], Line[{47, 48}]}, {GrayLevel[0.5], Line[{48, 49}]}, {GrayLevel[0.5], Line[{49, 50}]}, {GrayLevel[0.5], Line[{50, 51}]}, {GrayLevel[0.5], Line[{51, 52}]}, {GrayLevel[0.5], Line[{52, 53}]}, {GrayLevel[0.5], Line[{53, 54}]}, {GrayLevel[0.5], Line[{54, 55}]}, {GrayLevel[0.5], Line[{55, 56}]}, {GrayLevel[0.5], Line[{56, 57}]}, {GrayLevel[0.5], Line[{57, 58}]}, {GrayLevel[0.5], Line[{58, 59}]}, {GrayLevel[0.5], Line[{59, 60}]}, {GrayLevel[0.5], Line[{60, 61}]}, {GrayLevel[0.5], Line[{61, 62}]}, {GrayLevel[0.5], Line[{62, 63}]}, {GrayLevel[0.5], Line[{63, 64}]}, {GrayLevel[0.5], Line[{64, 65}]}, {GrayLevel[0.5], Line[{65, 66}]}, {GrayLevel[0.5], Line[{66, 67}]}, {GrayLevel[0.5], Line[{67, 68}]}, {GrayLevel[0.5], Line[{68, 69}]}, {GrayLevel[0.5], Line[{69, 70}]}, {GrayLevel[0.5], Line[{70, 71}]}, {GrayLevel[0.5], Line[{71, 72}]}, {GrayLevel[0.5], Line[{72, 73}]}, {GrayLevel[0.5], Line[{73, 74}]}, {GrayLevel[0.5], Line[{74, 75}]}, {GrayLevel[0.5], Line[{75, 76}]}, {GrayLevel[0.5], Line[{76, 77}]}, {GrayLevel[0.5], Line[{77, 78}]}, {GrayLevel[0.5], Line[{78, 79}]}, {GrayLevel[0.5], Line[{79, 80}]}, {GrayLevel[0.5], Line[{80, 81}]}, {GrayLevel[0.5], Line[{81, 82}]}, {GrayLevel[0.5], Line[{82, 83}]}, {GrayLevel[0.5], Line[{83, 84}]}, {GrayLevel[0.5], Line[{84, 85}]}, {GrayLevel[0.5], Line[{85, 86}]}, {GrayLevel[0.5], Line[{86, 87}]}, {GrayLevel[0.5], Line[{87, 88}]}, {GrayLevel[0.5], Line[{88, 89}]}, {GrayLevel[0.5], Line[{89, 90}]}, {GrayLevel[0.5], Line[{90, 91}]}, {GrayLevel[0.5], Line[{91, 92}]}, {GrayLevel[0.5], Line[{92, 93}]}, {GrayLevel[0.5], Line[{93, 94}]}, {GrayLevel[0.5], Line[{94, 95}]}, {GrayLevel[0.5], Line[{95, 96}]}, {GrayLevel[0.5], Line[{96, 97}]}, {GrayLevel[0.5], Line[{97, 98}]}, {GrayLevel[0.5], Line[{98, 99}]}, {GrayLevel[0.5], Line[{99, 100}]}, {GrayLevel[0.5], Line[{100, 101}]}, {GrayLevel[0.5], Line[{101, 102}]}, {GrayLevel[0.5], Line[{102, 103}]}, {GrayLevel[0.5], Line[{103, 104}]}, {GrayLevel[0.5], Line[{104, 105}]}, {GrayLevel[0.5], Line[{105, 106}]}}, $CellContext`x1 - 2 $CellContext`x2 == -2]}}, {{RGBColor[0, 0, 1], Inset[Style[5, FontWeight -> Bold], {-0.12321428571428603`, 0.888392857142857}]}, {RGBColor[0, 0, 1], Inset[Style[5, FontWeight -> Bold], {5.037499999999998, 3.468749999999999}]}}}, VertexColors -> {0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5}], GraphicsComplex[CompressedData["*) (*1:eJxd1U1IVFEYBuCb/ZdmtLQoJhKDFkKg24IWaStrWeAiaqVTqyAhKBAUAqFE*) (*Wmu4i8yFDggDGmThVDALS5mJysjIyCit6Vene96XeV+4F+WB4/n7vvOdY+r8*) (*5TMXq6Io6o5/g/zykxVb+w/FP0/l3lzH4QdDOcl+Nm7M9qw8kRjW+lg2D8y1*) (*tZ+Ylhz3SL7KpIuZ9EN5N56lLjcp2S8rT2HiCTkfz9pVMy7Zb1ROhGmLI3If*) (*4hiWWbTfkRx3S5ZGQlxnj1Vk+6Ccx77vSWxjcVyy35REWvptHGS8kWl5G3/I*) (*SY57JjFtU16yfVaieeCFxLaKLyXPqSA57q3sxUYsz9GulePv+oLEMnPvZAs6*) (*LsoOLGS5zgfJOD7KMr4lyTx+khz3WTLOZck4v0ikNfVVctw3ibLsshxn0e3G*) (*ikSYeYu07FmVnPe7ZB4t82g57w95KSybKUnm1XLen3I1TNv4S/Le/padyLvl*) (*uD8yzNr22rIu/0qeyz/JcWsS967Psn1dnsZAy/boeEWsu98i7mY7hHfAou4O*) (*bJBIU9pimSuWebV4x1JVEnU6bFEGzy23uVGyTizrxKJO6zdJrmtZN5Z1s1ky*) (*K3YmXJxrthAWum+5ry2yB3Vm+Y5YdJuyONZZW8K93irXE5YTMi6LYzqyTXbi*) (*YOxVBG657+0ylEF7t0WcyxZhNe6QXNfimS7Yk7gHOyXfIct1q2VDaD5qx8K5*) (*nrPvcU8s162RgyHwm/Yg7pFFXt9Y3qtdkuducW4zFu/YguW+ayXq+IJFmfZZ*) (*nnvCaLdM/v9P9v8P++EhaA==*) (*"], {{{}, {}, {Hue[0.67, 0.6, 0.6], Tooltip[{{GrayLevel[0.5], Line[{1, 2}]}, {GrayLevel[0.5], Line[{2, 3}]}, {GrayLevel[0.5], Line[{3, 4}]}, {GrayLevel[0.5], Line[{4, 5}]}, {GrayLevel[0.5], Line[{5, 6}]}, {GrayLevel[0.5], Line[{6, 7}]}, {GrayLevel[0.5], Line[{7, 8}]}, {GrayLevel[0.5], Line[{8, 9}]}, {GrayLevel[0.5], Line[{9, 10}]}, {GrayLevel[0.5], Line[{10, 11}]}, {GrayLevel[0.5], Line[{11, 12}]}, {GrayLevel[0.5], Line[{12, 13}]}, {GrayLevel[0.5], Line[{13, 14}]}, {GrayLevel[0.5], Line[{14, 15}]}, {GrayLevel[0.5], Line[{15, 16}]}, {GrayLevel[0.5], Line[{16, 17}]}, {GrayLevel[0.5], Line[{17, 18}]}, {GrayLevel[0.5], Line[{18, 19}]}, {GrayLevel[0.5], Line[{19, 20}]}, {GrayLevel[0.5], Line[{20, 21}]}, {GrayLevel[0.5], Line[{21, 22}]}, {GrayLevel[0.5], Line[{22, 23}]}, {GrayLevel[0.5], Line[{23, 24}]}, {GrayLevel[0.5], Line[{24, 25}]}, {GrayLevel[0.5], Line[{25, 26}]}, {GrayLevel[0.5], Line[{26, 27}]}, {GrayLevel[0.5], Line[{27, 28}]}, {GrayLevel[0.5], Line[{28, 29}]}, {GrayLevel[0.5], Line[{29, 30}]}, {GrayLevel[0.5], Line[{30, 31}]}, {GrayLevel[0.5], Line[{31, 32}]}, {GrayLevel[0.5], Line[{32, 33}]}, {GrayLevel[0.5], Line[{33, 34}]}, {GrayLevel[0.5], Line[{34, 35}]}, {GrayLevel[0.5], Line[{35, 36}]}, {GrayLevel[0.5], Line[{36, 37}]}, {GrayLevel[0.5], Line[{37, 38}]}, {GrayLevel[0.5], Line[{38, 39}]}, {GrayLevel[0.5], Line[{39, 40}]}, {GrayLevel[0.5], Line[{40, 41}]}, {GrayLevel[0.5], Line[{41, 42}]}, {GrayLevel[0.5], Line[{42, 43}]}, {GrayLevel[0.5], Line[{43, 44}]}, {GrayLevel[0.5], Line[{44, 45}]}, {GrayLevel[0.5], Line[{45, 46}]}, {GrayLevel[0.5], Line[{46, 47}]}, {GrayLevel[0.5], Line[{47, 48}]}, {GrayLevel[0.5], Line[{48, 49}]}, {GrayLevel[0.5], Line[{49, 50}]}, {GrayLevel[0.5], Line[{50, 51}]}, {GrayLevel[0.5], Line[{51, 52}]}, {GrayLevel[0.5], Line[{52, 53}]}, {GrayLevel[0.5], Line[{53, 54}]}, {GrayLevel[0.5], Line[{54, 55}]}, {GrayLevel[0.5], Line[{55, 56}]}, {GrayLevel[0.5], Line[{56, 57}]}, {GrayLevel[0.5], Line[{57, 58}]}, {GrayLevel[0.5], Line[{58, 59}]}, {GrayLevel[0.5], Line[{59, 60}]}, {GrayLevel[0.5], Line[{60, 61}]}, {GrayLevel[0.5], Line[{61, 62}]}, {GrayLevel[0.5], Line[{62, 63}]}, {GrayLevel[0.5], Line[{63, 64}]}, {GrayLevel[0.5], Line[{64, 65}]}, {GrayLevel[0.5], Line[{65, 66}]}, {GrayLevel[0.5], Line[{66, 67}]}, {GrayLevel[0.5], Line[{67, 68}]}, {GrayLevel[0.5], Line[{68, 69}]}, {GrayLevel[0.5], Line[{69, 70}]}, {GrayLevel[0.5], Line[{70, 71}]}, {GrayLevel[0.5], Line[{71, 72}]}, {GrayLevel[0.5], Line[{72, 73}]}, {GrayLevel[0.5], Line[{73, 74}]}, {GrayLevel[0.5], Line[{74, 75}]}, {GrayLevel[0.5], Line[{75, 76}]}, {GrayLevel[0.5], Line[{76, 77}]}, {GrayLevel[0.5], Line[{77, 78}]}, {GrayLevel[0.5], Line[{78, 79}]}, {GrayLevel[0.5], Line[{79, 80}]}, {GrayLevel[0.5], Line[{80, 81}]}, {GrayLevel[0.5], Line[{81, 82}]}, {GrayLevel[0.5], Line[{82, 83}]}, {GrayLevel[0.5], Line[{83, 84}]}, {GrayLevel[0.5], Line[{84, 85}]}, {GrayLevel[0.5], Line[{85, 86}]}, {GrayLevel[0.5], Line[{86, 87}]}, {GrayLevel[0.5], Line[{87, 88}]}, {GrayLevel[0.5], Line[{88, 89}]}, {GrayLevel[0.5], Line[{89, 90}]}, {GrayLevel[0.5], Line[{90, 91}]}, {GrayLevel[0.5], Line[{91, 92}]}, {GrayLevel[0.5], Line[{92, 93}]}, {GrayLevel[0.5], Line[{93, 94}]}, {GrayLevel[0.5], Line[{94, 95}]}, {GrayLevel[0.5], Line[{95, 96}]}, {GrayLevel[0.5], Line[{96, 97}]}, {GrayLevel[0.5], Line[{97, 98}]}, {GrayLevel[0.5], Line[{98, 99}]}, {GrayLevel[0.5], Line[{99, 100}]}, {GrayLevel[0.5], Line[{100, 101}]}, {GrayLevel[0.5], Line[{101, 102}]}, {GrayLevel[0.5], Line[{102, 103}]}, {GrayLevel[0.5], Line[{103, 104}]}, {GrayLevel[0.5], Line[{104, 105}]}, {GrayLevel[0.5], Line[{105, 106}]}, {GrayLevel[0.5], Line[{106, 107}]}, {GrayLevel[0.5], Line[{107, 108}]}, {GrayLevel[0.5], Line[{108, 109}]}, {GrayLevel[0.5], Line[{109, 110}]}, {GrayLevel[0.5], Line[{110, 111}]}, {GrayLevel[0.5], Line[{111, 112}]}, {GrayLevel[0.5], Line[{112, 113}]}, {GrayLevel[0.5], Line[{113, 114}]}, {GrayLevel[0.5], Line[{114, 115}]}, {GrayLevel[0.5], Line[{115, 116}]}, {GrayLevel[0.5], Line[{116, 117}]}, {GrayLevel[0.5], Line[{117, 118}]}, {GrayLevel[0.5], Line[{118, 119}]}, {GrayLevel[0.5], Line[{119, 120}]}, {GrayLevel[0.5], Line[{120, 121}]}, {GrayLevel[0.5], Line[{121, 122}]}, {GrayLevel[0.5], Line[{122, 123}]}, {GrayLevel[0.5], Line[{123, 124}]}, {GrayLevel[0.5], Line[{124, 125}]}, {GrayLevel[0.5], Line[{125, 126}]}, {GrayLevel[0.5], Line[{126, 127}]}, {GrayLevel[0.5], Line[{127, 128}]}, {GrayLevel[0.5], Line[{128, 129}]}}, $CellContext`x1 - $CellContext`x2 == 0]}}, {{RGBColor[0, 0, 1], Inset[Style[6, FontWeight -> Bold], {-0.11205357142857142`, -0.21205357142857142`}]}, {RGBColor[0, 0, 1], Inset[Style[6, FontWeight -> Bold], {4.06875, 3.96875}]}}}, VertexColors -> CompressedData["*) (*1:eJxTTMoPSmVkYGBoZoCBB/aj9Cg9So/SMBoAaK2VCw==*) (*"]], GraphicsComplex[CompressedData["*) (*1:eJxdVWuIlGUUntvO/dZSf1YwNlYUKoRgDCM4D1jZCoGJERSZhP7SqSAEl4KE*) (*YINIqG2h6EesZX+6GZVbgTGGCLkqbLil7pbNuo2ul213Z3f2Mjozzfuc882C*) (*HwzvfPOey3Oe85wznS+9sm13wOfzvdb8uFOfu6DncOHTY73lji+yeLj//NYd*) (*m04Wuvu6uvt2ZbHFnV0/Fzb1pEq57VmMNm97UkcLdEMWSX75tjC0Z92RQ7cz*) (*WPX1oWak52QNA2RwYTA/Npj/UtKlXDNyBs8499JR0bwZLDXcU5Cn3e+fpzHd*) (*4SKdFLo9n8b7Ln33kBDWcgp0zw2r/4EUQvwyIhccrqdSeNxd9/8hijsFxhn7*) (*U1YRYBJrjzh8o5Y/CeWiKC+6uJUE3i47AEX5yzlOJcAw3UUJOrOBBBYd3DfH*) (*hXxtSYDpzl+Wl12awTiepENJWO4ncbisx3pLhjdufF8R/lyOWT2Tdh8D6Whc*) (*k73O4IkYSHfpuuH17m8I3T+Ogu79U6I4o1bvfzLnypiOgHR2zlj8COr0n5EN*) (*JDSCD9zRM6v9ao+Y/6zaHw8r3gNlYZoPw2C5w2XDEzb+5uWXJh1DN9qMv3n5*) (*iXHajL95tUebxasI5XE4BOVtQbSOkPG5YHhDxteiMO5kEKxr/ZL2Jx8E43Qt*) (*C9v7aBDkrXfZ8AXh51lV/b0XAPtyqar92hAA4fXdUv2V/KAOcrctv9/y11R/*) (*G/2gzg7WLL53XxfSeNBn/arbvQ/eSb5PNFQXq32Wt2E4fCB/e+ri2tox1GSK*) (*/Nek5tp1r9/msCakK+/HYw7m3TVhun0rOD3cpym8W6r7zgBepWNVqM/DAWxz*) (*jrurivtswPyrLd4uspBlrbsniH4CXtK5HAtiJ8d4UXlbs9Inr29XXJ07FkT1*) (*EoLqoSIM076iA08XeQc3XJFTDvcbdt/Uld90VqGQ54Q0rQujSCBzovshjNdd*) (*omfLqqORMO6h36wscF4jtidmpW7vAXtv2LvP3r05aWejpoXtuD+KreR/SvaS*) (*xyhWU3A3ZT8LjJr/zdacPsABvSak4a0YPiLOSRnlXonhERJ7VVjO+pW9oPnj*) (*2MxB+Fe+4bqK4zLrnZDNxJGwPTouul8S5l9s7bXfSNQl4b57KImd9PtbfuBe*) (*TaLI/l6UCc5D0vxHWnu1Tr7PyYAj4J0UZpj/nNxH3aawnYa/K7//pCz+Wa3v*) (*wTTWkp8zqs9c2uKfsX6nUWG8E9q/U2m8QMNfhWnHPfvjrf8J34/O/nvV764M*) (*NjL+VzZ33v1nqq/vMuY/YHvMs39X996Ed+89Wdz5/3en/f+M97Dv*) (*"], {{{}, {}, {Hue[0.67, 0.6, 0.6], Tooltip[{{GrayLevel[0.5], Line[{1, 2}]}, {GrayLevel[0.5], Line[{2, 3}]}, {GrayLevel[0.5], Line[{3, 4}]}, {GrayLevel[0.5], Line[{4, 5}]}, {GrayLevel[0.5], Line[{5, 6}]}, {GrayLevel[0.5], Line[{6, 7}]}, {GrayLevel[0.5], Line[{7, 8}]}, {GrayLevel[0.5], Line[{8, 9}]}, {GrayLevel[0.5], Line[{9, 10}]}, {GrayLevel[0.5], Line[{10, 11}]}, {GrayLevel[0.5], Line[{11, 12}]}, {GrayLevel[0.5], Line[{12, 13}]}, {GrayLevel[0.5], Line[{13, 14}]}, {GrayLevel[0.5], Line[{14, 15}]}, {GrayLevel[0.5], Line[{15, 16}]}, {GrayLevel[0.5], Line[{16, 17}]}, {GrayLevel[0.5], Line[{17, 18}]}, {GrayLevel[0.5], Line[{18, 19}]}, {GrayLevel[0.5], Line[{19, 20}]}, {GrayLevel[0.5], Line[{20, 21}]}, {GrayLevel[0.5], Line[{21, 22}]}, {GrayLevel[0.5], Line[{22, 23}]}, {GrayLevel[0.5], Line[{23, 24}]}, {GrayLevel[0.5], Line[{24, 25}]}, {GrayLevel[0.5], Line[{25, 26}]}, {GrayLevel[0.5], Line[{26, 27}]}, {GrayLevel[0.5], Line[{27, 28}]}, {GrayLevel[0.5], Line[{28, 29}]}, {GrayLevel[0.5], Line[{29, 30}]}, {GrayLevel[0.5], Line[{30, 31}]}, {GrayLevel[0.5], Line[{31, 32}]}, {GrayLevel[0.5], Line[{32, 33}]}, {GrayLevel[0.5], Line[{33, 34}]}, {GrayLevel[0.5], Line[{34, 35}]}, {GrayLevel[0.5], Line[{35, 36}]}, {GrayLevel[0.5], Line[{36, 37}]}, {GrayLevel[0.5], Line[{37, 38}]}, {GrayLevel[0.5], Line[{38, 39}]}, {GrayLevel[0.5], Line[{39, 40}]}, {GrayLevel[0.5], Line[{40, 41}]}, {GrayLevel[0.5], Line[{41, 42}]}, {GrayLevel[0.5], Line[{42, 43}]}, {GrayLevel[0.5], Line[{43, 44}]}, {GrayLevel[0.5], Line[{44, 45}]}, {GrayLevel[0.5], Line[{45, 46}]}, {GrayLevel[0.5], Line[{46, 47}]}, {GrayLevel[0.5], Line[{47, 48}]}, {GrayLevel[0.5], Line[{48, 49}]}, {GrayLevel[0.5], Line[{49, 50}]}, {GrayLevel[0.5], Line[{50, 51}]}, {GrayLevel[0.5], Line[{51, 52}]}, {GrayLevel[0.5], Line[{52, 53}]}, {GrayLevel[0.5], Line[{53, 54}]}, {GrayLevel[0.5], Line[{54, 55}]}, {GrayLevel[0.5], Line[{55, 56}]}, {GrayLevel[0.5], Line[{56, 57}]}, {GrayLevel[0.5], Line[{57, 58}]}, {GrayLevel[0.5], Line[{58, 59}]}, {GrayLevel[0.5], Line[{59, 60}]}, {GrayLevel[0.5], Line[{60, 61}]}, {GrayLevel[0.5], Line[{61, 62}]}, {GrayLevel[0.5], Line[{62, 63}]}, {GrayLevel[0.5], Line[{63, 64}]}, {GrayLevel[0.5], Line[{64, 65}]}, {GrayLevel[0.5], Line[{65, 66}]}, {GrayLevel[0.5], Line[{66, 67}]}, {GrayLevel[0.5], Line[{67, 68}]}, {GrayLevel[0.5], Line[{68, 69}]}, {GrayLevel[0.5], Line[{69, 70}]}, {GrayLevel[0.5], Line[{70, 71}]}, {GrayLevel[0.5], Line[{71, 72}]}, {GrayLevel[0.5], Line[{72, 73}]}, {GrayLevel[0.5], Line[{73, 74}]}, {GrayLevel[0.5], Line[{74, 75}]}, {GrayLevel[0.5], Line[{75, 76}]}, {GrayLevel[0.5], Line[{76, 77}]}, {GrayLevel[0.5], Line[{77, 78}]}, {GrayLevel[0.5], Line[{78, 79}]}, {GrayLevel[0.5], Line[{79, 80}]}, {GrayLevel[0.5], Line[{80, 81}]}, {GrayLevel[0.5], Line[{81, 82}]}, {GrayLevel[0.5], Line[{82, 83}]}, {GrayLevel[0.5], Line[{83, 84}]}, {GrayLevel[0.5], Line[{84, 85}]}, {GrayLevel[0.5], Line[{85, 86}]}, {GrayLevel[0.5], Line[{86, 87}]}, {GrayLevel[0.5], Line[{87, 88}]}, {GrayLevel[0.5], Line[{88, 89}]}, {GrayLevel[0.5], Line[{89, 90}]}, {GrayLevel[0.5], Line[{90, 91}]}, {GrayLevel[0.5], Line[{91, 92}]}, {GrayLevel[0.5], Line[{92, 93}]}, {GrayLevel[0.5], Line[{93, 94}]}, {GrayLevel[0.5], Line[{94, 95}]}, {GrayLevel[0.5], Line[{95, 96}]}, {GrayLevel[0.5], Line[{96, 97}]}, {GrayLevel[0.5], Line[{97, 98}]}, {GrayLevel[0.5], Line[{98, 99}]}, {GrayLevel[0.5], Line[{99, 100}]}, {GrayLevel[0.5], Line[{100, 101}]}, {GrayLevel[0.5], Line[{101, 102}]}, {GrayLevel[0.5], Line[{102, 103}]}, {GrayLevel[0.5], Line[{103, 104}]}, {GrayLevel[0.5], Line[{104, 105}]}, {GrayLevel[0.5], Line[{105, 106}]}, {GrayLevel[0.5], Line[{106, 107}]}, {GrayLevel[0.5], Line[{107, 108}]}, {GrayLevel[0.5], Line[{108, 109}]}, {GrayLevel[0.5], Line[{109, 110}]}, {GrayLevel[0.5], Line[{110, 111}]}, {GrayLevel[0.5], Line[{111, 112}]}, {GrayLevel[0.5], Line[{112, 113}]}}, $CellContext`x1 + $CellContext`x2 == 4]}}, {{RGBColor[0, 0, 1], Inset[Style[7, FontWeight -> Bold], {4.2661184210526315`, -0.1661184210526316}]}, {RGBColor[0, 0, 1], Inset[Style[7, FontWeight -> Bold], {0.11398026315789479`, 3.9860197368421053`}]}}}, VertexColors -> {0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5}], GraphicsComplex[CompressedData["*) (*1:eJxdVX9o1HUY/t7v3+12F2Eoo9EyCgsiSYXF+xTGnGDCwjJtUYkhWeo0FBFy*) (*gkiDmHkMrBwUQwsciyVNktIblm5qybIp80c5Ozc7cnO7nXfbbXf2/bzv2z9+*) (*2bjt+3me93nf9/O871W+s7FurdOyrM32r/mUpy/56grznKKvW+qWPnHo12Si*) (*pmL7zX099Hhtoqo20ZvcV2X+OEuC/yXJ8BW/UeGvHXeO508k53V2O0Pz+vS8*) (*M9l+2DwX6KAJ15JI7q3O9D66t58aOU4bHf3gqv1zidq7Lq1a39FF5rQ6M6D8*) (*bhL+FbrHTw+J3p/K66MExxkkTrflAl3b9HNkwbVB5fcr/gbZhzbiMtmHNiJF*) (*AxzgOi3hAEP0U94UcJ2omP1jOQ0pf1D5Q/SGxkmkFjanFg5rvBQ1GLlNtxQ/*) (*TNK/NG05v+hW65Zhamua+8W6tjSVOP9hrSNNx1kvTbaYrfiv8m8rf4R8lbuj*) (*L/pHyPPQyrd2eUbpU65zVOseVf0xmr3hhxt7Zo8pf1z7NUYsf36c5v++OL5s*) (*/jjxvZ3JKH9c85ugZQzIKD+r/AzNWvPJM6dn3aV399dvvbx2gq4yMaf9ztIJ*) (*zj9HOwvpb778KKv8vPKzxPQ1k8T0/Xepg99PUUVPw5yehpziC4rPkZ/rLdB3*) (*xj6dOWKZWIE2sHCOvmV/TNMZY4/evPKLyp+kV8z11xXp7z2mI5PKKxLbtmqK*) (*ON18iTjdwpTyLQi/IP162UKcGzIt9XZY4DwS08TpVTog/p4RfqND+TP0oQkw*) (*4cBLHKCo506ITIlOc0OdeNO0o76k+i7ll+hhc/yZC9vM8dZ79JjRrXOD63jO*) (*gj0bdgPcWGLGscaSuI1uiD8t8H0870E/+8nC6/zeg+pmY1jFW16Ivxww1eyO*) (*emGby3aYAzIHXsg9O8D3l/Ji5Mi5p3+87ZD3ER8YFtC6Gn2q74SMrw8mu4rt*) (*TrxvCDmfxD3oBPv5Yz/Ezy7NJ6B8F1ZxgQHIvnFB/BqA+NUFLv+pILj8Crfq*) (*B5Wv/1shcDsPu4V3ICT3F/NA1kMIPP4LPIoPK98DHqdDYZzkBeHBFQ4QkTxq*) (*vWD/bIuA27vcq/oR5XtRbz4uRvCPWQ/NXrAdlz4A2UeKt8q0/z4cMAM4pwzr*) (*eEH4MGPWApWB98ROH0Z5QZQhywvCJ3xEle/HI8aPb0fB17fLL3k2RaXvR/2i*) (*fzYK2XcB1S/XfAPg+X+2HJ9zAgHp0+ZyqXskgFYGlMNsr6a5QeF3/88P4pTx*) (*szcGtseRILpMu1bHUMMGCeGOyX9jDLJPQ9qvmOYfAo/DQAxMPxdCl7m2J+No*) (*5f0a1nrjqhfGC7wg4niNDRsWncVxfM+6iv8qLvEvhvGeMUx7HCljx5th8P48*) (*Fsd69qXirQc1fgT3f9/dj/8PFwnp/A==*) (*"], {{{}, {}, {Hue[0.67, 0.6, 0.6], Tooltip[{{GrayLevel[0.5], Line[{1, 2}]}, {GrayLevel[0.5], Line[{2, 3}]}, {GrayLevel[0.5], Line[{3, 4}]}, {GrayLevel[0.5], Line[{4, 5}]}, {GrayLevel[0.5], Line[{5, 6}]}, {GrayLevel[0.5], Line[{6, 7}]}, {GrayLevel[0.5], Line[{7, 8}]}, {GrayLevel[0.5], Line[{8, 9}]}, {GrayLevel[0.5], Line[{9, 10}]}, {GrayLevel[0.5], Line[{10, 11}]}, {GrayLevel[0.5], Line[{11, 12}]}, {GrayLevel[0.5], Line[{12, 13}]}, {GrayLevel[0.5], Line[{13, 14}]}, {GrayLevel[0.5], Line[{14, 15}]}, {GrayLevel[0.5], Line[{15, 16}]}, {GrayLevel[0.5], Line[{16, 17}]}, {GrayLevel[0.5], Line[{17, 18}]}, {GrayLevel[0.5], Line[{18, 19}]}, {GrayLevel[0.5], Line[{19, 20}]}, {GrayLevel[0.5], Line[{20, 21}]}, {GrayLevel[0.5], Line[{21, 22}]}, {GrayLevel[0.5], Line[{22, 23}]}, {GrayLevel[0.5], Line[{23, 24}]}, {GrayLevel[0.5], Line[{24, 25}]}, {GrayLevel[0.5], Line[{25, 26}]}, {GrayLevel[0.5], Line[{26, 27}]}, {GrayLevel[0.5], Line[{27, 28}]}, {GrayLevel[0.5], Line[{28, 29}]}, {GrayLevel[0.5], Line[{29, 30}]}, {GrayLevel[0.5], Line[{30, 31}]}, {GrayLevel[0.5], Line[{31, 32}]}, {GrayLevel[0.5], Line[{32, 33}]}, {GrayLevel[0.5], Line[{33, 34}]}, {GrayLevel[0.5], Line[{34, 35}]}, {GrayLevel[0.5], Line[{35, 36}]}, {GrayLevel[0.5], Line[{36, 37}]}, {GrayLevel[0.5], Line[{37, 38}]}, {GrayLevel[0.5], Line[{38, 39}]}, {GrayLevel[0.5], Line[{39, 40}]}, {GrayLevel[0.5], Line[{40, 41}]}, {GrayLevel[0.5], Line[{41, 42}]}, {GrayLevel[0.5], Line[{42, 43}]}, {GrayLevel[0.5], Line[{43, 44}]}, {GrayLevel[0.5], Line[{44, 45}]}, {GrayLevel[0.5], Line[{45, 46}]}, {GrayLevel[0.5], Line[{46, 47}]}, {GrayLevel[0.5], Line[{47, 48}]}, {GrayLevel[0.5], Line[{48, 49}]}, {GrayLevel[0.5], Line[{49, 50}]}, {GrayLevel[0.5], Line[{50, 51}]}, {GrayLevel[0.5], Line[{51, 52}]}, {GrayLevel[0.5], Line[{52, 53}]}, {GrayLevel[0.5], Line[{53, 54}]}, {GrayLevel[0.5], Line[{54, 55}]}, {GrayLevel[0.5], Line[{55, 56}]}, {GrayLevel[0.5], Line[{56, 57}]}, {GrayLevel[0.5], Line[{57, 58}]}, {GrayLevel[0.5], Line[{58, 59}]}, {GrayLevel[0.5], Line[{59, 60}]}, {GrayLevel[0.5], Line[{60, 61}]}, {GrayLevel[0.5], Line[{61, 62}]}, {GrayLevel[0.5], Line[{62, 63}]}, {GrayLevel[0.5], Line[{63, 64}]}, {GrayLevel[0.5], Line[{64, 65}]}, {GrayLevel[0.5], Line[{65, 66}]}, {GrayLevel[0.5], Line[{66, 67}]}, {GrayLevel[0.5], Line[{67, 68}]}, {GrayLevel[0.5], Line[{68, 69}]}, {GrayLevel[0.5], Line[{69, 70}]}, {GrayLevel[0.5], Line[{70, 71}]}, {GrayLevel[0.5], Line[{71, 72}]}, {GrayLevel[0.5], Line[{72, 73}]}, {GrayLevel[0.5], Line[{73, 74}]}, {GrayLevel[0.5], Line[{74, 75}]}, {GrayLevel[0.5], Line[{75, 76}]}, {GrayLevel[0.5], Line[{76, 77}]}, {GrayLevel[0.5], Line[{77, 78}]}, {GrayLevel[0.5], Line[{78, 79}]}, {GrayLevel[0.5], Line[{79, 80}]}, {GrayLevel[0.5], Line[{80, 81}]}, {GrayLevel[0.5], Line[{81, 82}]}, {GrayLevel[0.5], Line[{82, 83}]}, {GrayLevel[0.5], Line[{83, 84}]}, {GrayLevel[0.5], Line[{84, 85}]}, {GrayLevel[0.5], Line[{85, 86}]}, {GrayLevel[0.5], Line[{86, 87}]}, {GrayLevel[0.5], Line[{87, 88}]}, {GrayLevel[0.5], Line[{88, 89}]}, {GrayLevel[0.5], Line[{89, 90}]}, {GrayLevel[0.5], Line[{90, 91}]}, {GrayLevel[0.5], Line[{91, 92}]}, {GrayLevel[0.5], Line[{92, 93}]}, {GrayLevel[0.5], Line[{93, 94}]}, {GrayLevel[0.5], Line[{94, 95}]}, {GrayLevel[0.5], Line[{95, 96}]}, {GrayLevel[0.5], Line[{96, 97}]}, {GrayLevel[0.5], Line[{97, 98}]}, {GrayLevel[0.5], Line[{98, 99}]}, {GrayLevel[0.5], Line[{99, 100}]}, {GrayLevel[0.5], Line[{100, 101}]}, {GrayLevel[0.5], Line[{101, 102}]}, {GrayLevel[0.5], Line[{102, 103}]}, {GrayLevel[0.5], Line[{103, 104}]}, {GrayLevel[0.5], Line[{104, 105}]}, {GrayLevel[0.5], Line[{105, 106}]}, {GrayLevel[0.5], Line[{106, 107}]}, {GrayLevel[0.5], Line[{107, 108}]}, {GrayLevel[0.5], Line[{108, 109}]}, {GrayLevel[0.5], Line[{109, 110}]}, {GrayLevel[0.5], Line[{110, 111}]}, {GrayLevel[0.5], Line[{111, 112}]}}, 2 $CellContext`x1 - 3 $CellContext`x2 == -1]}}, {{RGBColor[0, 0, 1], Inset[Style[8, FontWeight -> Bold], {-0.12379032258064518`, 0.18413978494623656`}]}, {RGBColor[0, 0, 1], Inset[Style[8, FontWeight -> Bold], {5.030803571428569, 3.620535714285713}]}}}, VertexColors -> {0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5}], GraphicsComplex[CompressedData["*) (*1:eJxdVU1oU2kUfflPmrz8utKVUpjOYhCECoJwLxTUFMQqgouBMMi4qp1Bd3Vj*) (*QVAQRI0FwYVtRUUUrchMRMgQURhsEaYwTKtm8AcatNraNm3TJmkT33fvycZA*) (*OaTvnO/ce+79XrYe/f3QMbdlWSecP4P62cSKcc62p7Ptr1NsIP0oxgczXf32*) (*fykeyZ8tb74bU95ESnE4Bl6K9xnhrzHwUuyoSp2HwR8Gn2PgpfiyEY5FwUux*) (*oxrvvRMFD/yBKJc6B6d6Gkke6+0YHVm3oUvyNuObtaFLckSINnRJ9GOzo8qf*) (*PZTkA6aPWxF+k+sr5n5K8h8O9P0cwbkJvjLVk+mqhlX3NIFzwmwe528neNJ5*) (*3L8/rP1mE7xXggqzLb4JdqpwmG3whd5qY2NfrsXZ2JdWQmxsi5W41jEX4vWm*) (*8/knzmsGh0LcI3XGNcfuEOekzjgfN7pckE27HQfiPG1srweRVxz1BvkHMYwB*) (*A9ww556OAQPaz54YMIA6MScrwA+M7FoU6Ne801GgX/3nbaAPvq3cfWzKynTZ*) (*QJ/mlbSBPjbjHrwaAXrhGwF62bQ9/iUM9Gj/xTDQw6aNJoWBHq7jex3fZb9u*) (*tgHdul+DbUA36sScBtx8zjTyKQR0scTcFwK6WMa6OwR0oc4Q0KXzvxQEWvyb*) (*Ee4MAi29P6UA35CGmqS+rdybJDHtCuienmxQaw6KDXplFuKCX/do8wae+/F8*) (*nbaYwj74eKc56HmdZC3HkHu2TmLf6dM8e2skub3z8obBQpWU52WB9ipJnee9*) (*mtemKuptzalKsu5vPZpHfpUey356tN9chaSPHR7t+1gF9Xqgr9BxIbhZ+rJX*) (*SH3d+t5JL5P0+cal9+XPJfi3cl8i2bvtLn0vZMoktnMW9mORJIczlu63fxH+*) (*FvQLqh9tktyfzLzq+5qk+q8ka3KkAd4s9BuaW3OGZE/u1Env/UfkVyPZt/YS*) (*/l8jfR9MQ18DTuu8q2uaU/kD+l/TuaXf6zxHVpHTW+hW4V/U+XRXSMbZ8Zo0*) (*9xXUP4lzl0nWZfxf6JeBE5rvtSXtc+olfJdI7/e45tdVpodywDPoFqkq/qP0*) (*l9zLBbpoaP33df7FBXolBdxTPi/QsvAval4d8/SjGhf0vHmqyfOhgupnqVsK*) (*eVIQ+/5Z+kUWMA/+LLkEnxb+F5/PqPvvgt6rGTol/bwAfwb8l/j+qdV/4fvf*) (*t+/53wBsIscl*) (*"], {{{}, {}, {Hue[0.67, 0.6, 0.6], Tooltip[{{GrayLevel[0.5], Line[{1, 2}]}, {GrayLevel[0.5], Line[{2, 3}]}, {GrayLevel[0.5], Line[{3, 4}]}, {GrayLevel[0.5], Line[{4, 5}]}, {GrayLevel[0.5], Line[{5, 6}]}, {GrayLevel[0.5], Line[{6, 7}]}, {GrayLevel[0.5], Line[{7, 8}]}, {GrayLevel[0.5], Line[{8, 9}]}, {GrayLevel[0.5], Line[{9, 10}]}, {GrayLevel[0.5], Line[{10, 11}]}, {GrayLevel[0.5], Line[{11, 12}]}, {GrayLevel[0.5], Line[{12, 13}]}, {GrayLevel[0.5], Line[{13, 14}]}, {GrayLevel[0.5], Line[{14, 15}]}, {GrayLevel[0.5], Line[{15, 16}]}, {GrayLevel[0.5], Line[{16, 17}]}, {GrayLevel[0.5], Line[{17, 18}]}, {GrayLevel[0.5], Line[{18, 19}]}, {GrayLevel[0.5], Line[{19, 20}]}, {GrayLevel[0.5], Line[{20, 21}]}, {GrayLevel[0.5], Line[{21, 22}]}, {GrayLevel[0.5], Line[{22, 23}]}, {GrayLevel[0.5], Line[{23, 24}]}, {GrayLevel[0.5], Line[{24, 25}]}, {GrayLevel[0.5], Line[{25, 26}]}, {GrayLevel[0.5], Line[{26, 27}]}, {GrayLevel[0.5], Line[{27, 28}]}, {GrayLevel[0.5], Line[{28, 29}]}, {GrayLevel[0.5], Line[{29, 30}]}, {GrayLevel[0.5], Line[{30, 31}]}, {GrayLevel[0.5], Line[{31, 32}]}, {GrayLevel[0.5], Line[{32, 33}]}, {GrayLevel[0.5], Line[{33, 34}]}, {GrayLevel[0.5], Line[{34, 35}]}, {GrayLevel[0.5], Line[{35, 36}]}, {GrayLevel[0.5], Line[{36, 37}]}, {GrayLevel[0.5], Line[{37, 38}]}, {GrayLevel[0.5], Line[{38, 39}]}, {GrayLevel[0.5], Line[{39, 40}]}, {GrayLevel[0.5], Line[{40, 41}]}, {GrayLevel[0.5], Line[{41, 42}]}, {GrayLevel[0.5], Line[{42, 43}]}, {GrayLevel[0.5], Line[{43, 44}]}, {GrayLevel[0.5], Line[{44, 45}]}, {GrayLevel[0.5], Line[{45, 46}]}, {GrayLevel[0.5], Line[{46, 47}]}, {GrayLevel[0.5], Line[{47, 48}]}, {GrayLevel[0.5], Line[{48, 49}]}, {GrayLevel[0.5], Line[{49, 50}]}, {GrayLevel[0.5], Line[{50, 51}]}, {GrayLevel[0.5], Line[{51, 52}]}, {GrayLevel[0.5], Line[{52, 53}]}, {GrayLevel[0.5], Line[{53, 54}]}, {GrayLevel[0.5], Line[{54, 55}]}, {GrayLevel[0.5], Line[{55, 56}]}, {GrayLevel[0.5], Line[{56, 57}]}, {GrayLevel[0.5], Line[{57, 58}]}, {GrayLevel[0.5], Line[{58, 59}]}, {GrayLevel[0.5], Line[{59, 60}]}, {GrayLevel[0.5], Line[{60, 61}]}, {GrayLevel[0.5], Line[{61, 62}]}, {GrayLevel[0.5], Line[{62, 63}]}, {GrayLevel[0.5], Line[{63, 64}]}, {GrayLevel[0.5], Line[{64, 65}]}, {GrayLevel[0.5], Line[{65, 66}]}, {GrayLevel[0.5], Line[{66, 67}]}, {GrayLevel[0.5], Line[{67, 68}]}, {GrayLevel[0.5], Line[{68, 69}]}, {GrayLevel[0.5], Line[{69, 70}]}, {GrayLevel[0.5], Line[{70, 71}]}, {GrayLevel[0.5], Line[{71, 72}]}, {GrayLevel[0.5], Line[{72, 73}]}, {GrayLevel[0.5], Line[{73, 74}]}, {GrayLevel[0.5], Line[{74, 75}]}, {GrayLevel[0.5], Line[{75, 76}]}, {GrayLevel[0.5], Line[{76, 77}]}, {GrayLevel[0.5], Line[{77, 78}]}, {GrayLevel[0.5], Line[{78, 79}]}, {GrayLevel[0.5], Line[{79, 80}]}, {GrayLevel[0.5], Line[{80, 81}]}, {GrayLevel[0.5], Line[{81, 82}]}, {GrayLevel[0.5], Line[{82, 83}]}, {GrayLevel[0.5], Line[{83, 84}]}, {GrayLevel[0.5], Line[{84, 85}]}, {GrayLevel[0.5], Line[{85, 86}]}, {GrayLevel[0.5], Line[{86, 87}]}, {GrayLevel[0.5], Line[{87, 88}]}, {GrayLevel[0.5], Line[{88, 89}]}, {GrayLevel[0.5], Line[{89, 90}]}, {GrayLevel[0.5], Line[{90, 91}]}, {GrayLevel[0.5], Line[{91, 92}]}, {GrayLevel[0.5], Line[{92, 93}]}, {GrayLevel[0.5], Line[{93, 94}]}, {GrayLevel[0.5], Line[{94, 95}]}, {GrayLevel[0.5], Line[{95, 96}]}, {GrayLevel[0.5], Line[{96, 97}]}, {GrayLevel[0.5], Line[{97, 98}]}, {GrayLevel[0.5], Line[{98, 99}]}, {GrayLevel[0.5], Line[{99, 100}]}, {GrayLevel[0.5], Line[{100, 101}]}, {GrayLevel[0.5], Line[{101, 102}]}, {GrayLevel[0.5], Line[{102, 103}]}, {GrayLevel[0.5], Line[{103, 104}]}, {GrayLevel[0.5], Line[{104, 105}]}, {GrayLevel[0.5], Line[{105, 106}]}, {GrayLevel[0.5], Line[{106, 107}]}, {GrayLevel[0.5], Line[{107, 108}]}, {GrayLevel[0.5], Line[{108, 109}]}, {GrayLevel[0.5], Line[{109, 110}]}, {GrayLevel[0.5], Line[{110, 111}]}}, 2 $CellContext`x1 - 2 $CellContext`x2 == 2]}}, {{RGBColor[0, 0, 1], Inset[Style[9, FontWeight -> Bold], {5.062053571428571, 3.962053571428571}]}, {RGBColor[0, 0, 1], Inset[Style[9, FontWeight -> Bold], {0.88125, -0.21875000000000003`}]}}}, VertexColors -> {0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5}], GraphicsComplex[CompressedData["*) (*1:eJxd1DlIA0EYBeD1vuKBjQFBUFJZSCobhQgWujaCohY2JiBYGC1s1MZKKxtT*) (*2Whl4dFYSEARRK0iqAFRJEEwiCYSFDXxiHfmDXk/2WHD8JHZnfnfzG69Z7xn*) (*ONcwjIH0T/W6Pbl0H9zL2G2fS18B+tBpUxfd5XOYPsc2XauGu1fp+KzyiCvj*) (*S7837Pdu0Gvrqu3QidbT9HVIL6jHm0c0luMO0tGp5fR1RuPx4Qsa09tDdB/a*) (*FY3lm9c0HheN0Lr+G7oTN9zSOo47y/pi9B/aPX2HCeL0CQp+sKz/0fL/E32A*) (*CZ8t48WIc/2FRjlTCUu9SRrLtb/S07hBPKYe73+jPSjwne5Hnh+0iQlSNLYz*) (*kbLk9UljOb4vS37fNIY7f2jd/9JJTCCOId+/7PU3G20Z6/MiRvmFOTTK9YpR*) (*3rkY5bTm0jgPK2K9P3k0ljch1vslxvD2fFrvnxj5VBfQSyhAXKPyioj1+1JI*) (*l6n8NsX6fSiif9XxnBFPYv/FL6qA7mJ6FAGJcZ7rSughbKA4pIY/iHsRaCl9*) (*rPLcFXfg/JTR+yqgeXELNkC8pfIctNFNODDiVRVnYzndgPMmRp4pMfJ0V2Tn*) (*GRAjT2dldp6LYt1X0dbvZ+b7+A8neqQb*) (*"], {{{}, {}, {Hue[0.67, 0.6, 0.6], Tooltip[{{GrayLevel[0.5], Line[{1, 2}]}, {GrayLevel[0.5], Line[{2, 3}]}, {GrayLevel[0.5], Line[{3, 4}]}, {GrayLevel[0.5], Line[{4, 5}]}, {GrayLevel[0.5], Line[{5, 6}]}, {GrayLevel[0.5], Line[{6, 7}]}, {GrayLevel[0.5], Line[{7, 8}]}, {GrayLevel[0.5], Line[{8, 9}]}, {GrayLevel[0.5], Line[{9, 10}]}, {GrayLevel[0.5], Line[{10, 11}]}, {GrayLevel[0.5], Line[{11, 12}]}, {GrayLevel[0.5], Line[{12, 13}]}, {GrayLevel[0.5], Line[{13, 14}]}, {GrayLevel[0.5], Line[{14, 15}]}, {GrayLevel[0.5], Line[{15, 16}]}, {GrayLevel[0.5], Line[{16, 17}]}, {GrayLevel[0.5], Line[{17, 18}]}, {GrayLevel[0.5], Line[{18, 19}]}, {GrayLevel[0.5], Line[{19, 20}]}, {GrayLevel[0.5], Line[{20, 21}]}, {GrayLevel[0.5], Line[{21, 22}]}, {GrayLevel[0.5], Line[{22, 23}]}, {GrayLevel[0.5], Line[{23, 24}]}, {GrayLevel[0.5], Line[{24, 25}]}, {GrayLevel[0.5], Line[{25, 26}]}, {GrayLevel[0.5], Line[{26, 27}]}, {GrayLevel[0.5], Line[{27, 28}]}, {GrayLevel[0.5], Line[{28, 29}]}, {GrayLevel[0.5], Line[{29, 30}]}, {GrayLevel[0.5], Line[{30, 31}]}, {GrayLevel[0.5], Line[{31, 32}]}, {GrayLevel[0.5], Line[{32, 33}]}, {GrayLevel[0.5], Line[{33, 34}]}, {GrayLevel[0.5], Line[{34, 35}]}, {GrayLevel[0.5], Line[{35, 36}]}, {GrayLevel[0.5], Line[{36, 37}]}, {GrayLevel[0.5], Line[{37, 38}]}, {GrayLevel[0.5], Line[{38, 39}]}, {GrayLevel[0.5], Line[{39, 40}]}, {GrayLevel[0.5], Line[{40, 41}]}, {GrayLevel[0.5], Line[{41, 42}]}, {GrayLevel[0.5], Line[{42, 43}]}, {GrayLevel[0.5], Line[{43, 44}]}, {GrayLevel[0.5], Line[{44, 45}]}, {GrayLevel[0.5], Line[{45, 46}]}, {GrayLevel[0.5], Line[{46, 47}]}, {GrayLevel[0.5], Line[{47, 48}]}, {GrayLevel[0.5], Line[{48, 49}]}, {GrayLevel[0.5], Line[{49, 50}]}, {GrayLevel[0.5], Line[{50, 51}]}, {GrayLevel[0.5], Line[{51, 52}]}, {GrayLevel[0.5], Line[{52, 53}]}, {GrayLevel[0.5], Line[{53, 54}]}, {GrayLevel[0.5], Line[{54, 55}]}, {GrayLevel[0.5], Line[{55, 56}]}, {GrayLevel[0.5], Line[{56, 57}]}, {GrayLevel[0.5], Line[{57, 58}]}, {GrayLevel[0.5], Line[{58, 59}]}, {GrayLevel[0.5], Line[{59, 60}]}, {GrayLevel[0.5], Line[{60, 61}]}, {GrayLevel[0.5], Line[{61, 62}]}, {GrayLevel[0.5], Line[{62, 63}]}, {GrayLevel[0.5], Line[{63, 64}]}, {GrayLevel[0.5], Line[{64, 65}]}, {GrayLevel[0.5], Line[{65, 66}]}, {GrayLevel[0.5], Line[{66, 67}]}, {GrayLevel[0.5], Line[{67, 68}]}, {GrayLevel[0.5], Line[{68, 69}]}, {GrayLevel[0.5], Line[{69, 70}]}, {GrayLevel[0.5], Line[{70, 71}]}, {GrayLevel[0.5], Line[{71, 72}]}, {GrayLevel[0.5], Line[{72, 73}]}, {GrayLevel[0.5], Line[{73, 74}]}, {GrayLevel[0.5], Line[{74, 75}]}, {GrayLevel[0.5], Line[{75, 76}]}, {GrayLevel[0.5], Line[{76, 77}]}, {GrayLevel[0.5], Line[{77, 78}]}, {GrayLevel[0.5], Line[{78, 79}]}, {GrayLevel[0.5], Line[{79, 80}]}, {GrayLevel[0.5], Line[{80, 81}]}, {GrayLevel[0.5], Line[{81, 82}]}, {GrayLevel[0.5], Line[{82, 83}]}, {GrayLevel[0.5], Line[{83, 84}]}, {GrayLevel[0.5], Line[{84, 85}]}}, 2 $CellContext`x1 == 2]}}, {{RGBColor[0, 0, 1], Inset[Style[10, FontWeight -> Bold], {1.1, 3.9}], Inset[Style[10, FontWeight -> Bold], {1.1, -0.2}]}, {RGBColor[0, 0, 1], Inset[Style[10, FontWeight -> Bold], {1.1, 3.9}], Inset[Style[10, FontWeight -> Bold], {1.1, -0.2}]}}}, VertexColors -> {0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5}], {Inset[24., {0, 1.}, {0, -1}], Inset[24., {1., 0}, {0, -1}], Inset[14., {1., 1.}, {0, -1}], Inset[14., {1., 1.5}, {0, -1}], Inset[20., {1., 3.}, {0, -1}], Inset[10.4, {1.5999999999999999`, 1.8}, {0, -1}], Inset[10., {1.666666666667, 1.666666666667}, {0, -1}], Inset[10., {1.75, 1.5}, {0, -1}], Inset[14., {2., 1.}, {0, -1}], Inset[10., {2., 1.666666666667}, {0, -1}], Inset[10., {2., 2.}, {0, -1}], Inset[18., {2., 3.}, {0, -1}], Inset[10.4, {2.2, 1.8}, {0, -1}], Inset[14., {2.5, 1.5}, {0, -1}], Inset[24., {3., 1.}, {0, -1}], Inset[18., {3., 3.}, {0, -1}], Inset[20., {4., 3.}, {0, -1}]}}, AspectRatio -> Automatic, Axes -> True, AxesStyle -> Directive[RGBColor[1, 0.5, 0], Opacity[0.5]], ImageSize -> {422., Automatic}, TicksStyle -> Directive[RGBColor[1, 0.5, 0], Opacity[0.5]]]*) (* ::Input:: *) (*X={{-2,-1},{-1,0},{0,-3},{0,-1},{1,-2},{1,-1},{1,1},{2,-3},{2,-2},{2,0}};*) (*z={-5,-2,-9,-1,-2,0,4,-1,2,2};*) (* ::Text:: *) (*Algorithm 3 starting from intersection of lines 4 and 10.*) (* ::Input:: *) (*rj=Alg3[X,z,{4,10}];*) (*a=Table[ToExpression["a"<>ToString[i]],{i,Length[X[[1]]]}];*) (*Print["NMinimize:",NMinimize[F[a,X,z],a]];*) (*Print["Algoritam 3:",{F[rj,X,z],rj}];*) (* ::Text:: *) (*Algorithm 3 starting from intersection of lines 5 and 10.*) (* ::Input:: *) (*rj=Alg3[X,z,{5,10}];*) (*a=Table[ToExpression["a"<>ToString[i]],{i,Length[X[[1]]]}];*) (*Print["NMinimize:",NMinimize[F[a,X,z],a]];*) (*Print["Algoritam 3:",{F[rj,X,z],rj}];*)