(************** Content-type: application/mathematica ************** CreatedBy='Mathematica 4.2' Mathematica-Compatible Notebook This notebook can be used with any Mathematica-compatible application, such as Mathematica, MathReader or Publicon. The data for the notebook starts with the line containing stars above. To get the notebook into a Mathematica-compatible application, do one of the following: * Save the data starting with the line of stars above into a file with a name ending in .nb, then open the file inside the application; * Copy the data starting with the line of stars above to the clipboard, then use the Paste menu command inside the application. Data for notebooks contains only printable 7-bit ASCII and can be sent directly in email or through ftp in text mode. Newlines can be CR, LF or CRLF (Unix, Macintosh or MS-DOS style). NOTE: If you modify the data for this notebook not in a Mathematica- compatible application, you must delete the line below containing the word CacheID, otherwise Mathematica-compatible applications may try to use invalid cache data. For more information on notebooks and Mathematica-compatible applications, contact Wolfram Research: web: http://www.wolfram.com email: info@wolfram.com phone: +1-217-398-0700 (U.S.) Notebook reader applications are available free of charge from Wolfram Research. *******************************************************************) (*CacheID: 232*) (*NotebookFileLineBreakTest NotebookFileLineBreakTest*) (*NotebookOptionsPosition[ 241747, 6769]*) (*NotebookOutlinePosition[ 242826, 6803]*) (* CellTagsIndexPosition[ 242782, 6799]*) (*WindowFrame->Normal*) Notebook[{ Cell[CellGroupData[{ Cell["\<\ Incremental/Iterative Solution Methods For Two-Bar Arch Modeled by \ the Total Lagrangian Formulation. \ \>", "Section", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, Evaluatable->False, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, Background->RGBColor[1, 1, 0]], Cell[TextData[ "This Notebook is a testbed of incremental /iterative (predictor/\ncorrector) \ solution methods, converted from ancient Fortran. The methods \nare applied \ to the 2-bar arch problem of Chapter 8. \nIt can used for the Exercises of \ Chapter 20. \n\nThe present implementation includes three predictors\n \ FE: Forward Euler\n MR: Midpoint Rule (a 2-step Runge Kutta)\n \ RK4: Classical 4-step Runge Kutta\n with three Increment Control \ Strategies\n LC Lambda Control, also called Load Control\n DC \ Displacement Control\n AC ArcLength Control\n The \ implementations below use a constant stepsize ell and a positive-work \ criterion to\n traverse limit points.\n \n To run the main programs in Cells \ 12ff, Cells 1 through 12 must be initialized. "], "Text", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, Evaluatable->False, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, FontSize->14, Background->RGBColor[0, 1, 0]], Cell[TextData[{ StyleBox["Cell 1. The modules ", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, Evaluatable->False, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, Background->RGBColor[0, 1, 1]], StyleBox["FormTanStiffnessOfTwoBarArch", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, Evaluatable->False, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, FontWeight->"Bold", Background->RGBColor[0, 1, 1]], StyleBox[" and ", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, Evaluatable->False, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, Background->RGBColor[0, 1, 1]], StyleBox["FormInternalForceOfTwoBarArch", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, Evaluatable->False, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, FontWeight->"Bold", Background->RGBColor[0, 1, 1]], StyleBox[ "\n below form return the internal force p and tangent stiffness K for the \ two-bar arch structure of Chapter 8,\n which is discretized by two Total \ Lagrangian (TL) bar elements.\n They are extracted directly from the \ Mathematica modules of that Chapter, with minor edits.\n", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, Evaluatable->False, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, Background->RGBColor[0, 1, 1]], StyleBox[" DetTanStiffness", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, Evaluatable->False, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, FontWeight->"Bold", Background->RGBColor[0, 1, 1]], StyleBox[" returns exactly that.\n", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, Evaluatable->False, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, Background->RGBColor[0, 1, 1]], StyleBox[" FormIncrementalLoadVector", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, Evaluatable->False, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, FontWeight->"Bold", Background->RGBColor[0, 1, 1]], StyleBox[ " simply returns its argument force as q. This trivial operation\n is \ implemented as a module to have a placeholder for more complicated cases.\n \n\ The inputs to ", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, Evaluatable->False, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, Background->RGBColor[0, 1, 1]], StyleBox["FormTanStiffnessOfTwoBarArch", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, Evaluatable->False, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, FontWeight->"Bold", Background->RGBColor[0, 1, 1]], StyleBox[" and ", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, Evaluatable->False, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, Background->RGBColor[0, 1, 1]], StyleBox["FormInternalForceOfTwoBarArch\n ", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, Evaluatable->False, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, FontWeight->"Bold", Background->RGBColor[0, 1, 1]], StyleBox[ "includes four structural properties collected in sprop:\n S \ arch span between supports\n H arch height \ under zero load\n Em bar elastic modulus\n A0 \ cross section area in reference state\n and the displacements \ u={uX,uY} of the crown from the reference state", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, Evaluatable->False, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, Background->RGBColor[0, 1, 1]] }], "Text", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, Evaluatable->False, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, Background->RGBColor[0, 1, 1]], Cell[CellGroupData[{ Cell[TextData[{ StyleBox[ "FormTanStiffnessOfTwoBarArch[{S_,H_,Em_,A0_},{uX_,uY_}]:=\n \ Module[{c,K},c=4*Em*A0/(4*H^2+S^2)^(3/2); \n \ K=c*{{(S^2+6*uX^2+4*H*uY+2*uY^2), \n 4*uX*(H+uY)},{4*uX*(H+uY), \n \ 2*(2*H^2+uX^2+6*H*uY+3*uY^2)}};\n Return[K]];\n \n\ FormInternalForceOfTwoBarArch[{S_,H_,Em_,A0_},{uX_,uY_}]:=\n \ Module[{c,p},c=4*Em*A0/(4*H^2+S^2)^(3/2); \n \ p=c*{uX*(S^2+2*uX^2+4*H*uY+2*uY^2),\n 2*(H+uY)*(uX^2+2*H*uY+uY^2)};\n \ Return[p]];\n \nFormIncrementalLoadVector[lambda_,force_,u_]:=Module[{},\n\ Return[force]];\n \nDetTanStiffness[sprop_,u_]:=Module[{K11,K12,K21,K22},\n\ {{K11,K12},{K21,K22}}=FormTanStiffnessOfTwoBarArch[sprop,u];\n \ Return[K11*K22-K12*K21]];\n ", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, InitializationCell->True, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, Background->RGBColor[1, 0.647715, 0.521981]], StyleBox[ "\nClearAll[uX,uY];\nK= FormTanStiffnessOfTwoBarArch[{S,H,Em,A0},{uX,uY}];\n\ p= FormInternalForceOfTwoBarArch[{S,H,Em,A0},{uX,uY}];\nKK={D[p,uX],D[p,uY]};\ \nPrint[\"Check K=grad(p): \",Simplify[K-KK]];\n\ Print[Simplify[Together[DetTanStiffness[{S,H,Em,A0},{0,uY}]]]];\n", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, InitializationCell->True, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, FontColor->RGBColor[1, 0, 0], Background->RGBColor[1, 0.647715, 0.521981]] }], "Input", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, InitializationCell->True, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, Background->RGBColor[1, 0.647715, 0.521981]], Cell[BoxData[ InterpretationBox[ \("Check K=grad(p): "\[InvisibleSpace]{{0, 0}, {0, 0}}\), SequenceForm[ "Check K=grad(p): ", {{0, 0}, {0, 0}}], Editable->False]], "Print"], Cell[BoxData[ \(\(32\ A0\^2\ Em\^2\ \((2\ H\^2 + 6\ H\ uY + 3\ uY\^2)\)\ \((S\^2 + 2\ uY\ \((2\ H + uY)\))\)\)\/\((4\ H\^2 + S\^2)\)\^3\)], "Print"] }, Open ]], Cell[TextData[{ StyleBox["Cell 2. ", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, Evaluatable->False, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, Background->RGBColor[0, 1, 1]], StyleBox["SolveTwoLinearEqs", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, Evaluatable->False, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, FontWeight->"Bold", Background->RGBColor[0, 1, 1]], StyleBox[ " is an ad-hoc module that solves a system of two linear equations: \n\ Ax=b. It returns x and a regularity/singularity indicator.", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, Evaluatable->False, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, Background->RGBColor[0, 1, 1]] }], "Text", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, Evaluatable->False, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, Background->RGBColor[0, 1, 1]], Cell[CellGroupData[{ Cell[TextData[{ StyleBox[ "SolveTwoLinearEqs[A_,b_]:=Module[{A11,A12,A21,A22,b1,b2,\n \ dnorm,Adet,x,eps=10.^(-15)},\n {{A11,A12},{A21,A22}}=A; \ Adet=A11*A22-A12*A21; {b1,b2}=b;\n dnorm=Max[Abs[A11*A22],Abs[A12*A21]];\n \ If [Abs[Adet]<=eps*dnorm, \n sing=True; Return[{Null,False}]];\n \ x={A22*b1-A12*b2,-A21*b1+A11*b2}/Adet; \n Return[{x,True}]];\n", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, InitializationCell->True, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, Background->RGBColor[1, 0.647715, 0.521981]], StyleBox[ " \nPrint[SolveTwoLinearEqs[{{2,6},{1,3}},{8,4}]];\n\ Print[SolveTwoLinearEqs[{{2,5},{1,3}},{7,4}]];", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, InitializationCell->True, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, FontColor->RGBColor[1, 0, 0], Background->RGBColor[1, 0.647715, 0.521981]] }], "Input", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, InitializationCell->True, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, Background->RGBColor[1, 0.647715, 0.521981]], Cell[BoxData[ \({Null, False}\)], "Print"], Cell[BoxData[ \({{1, 1}, True}\)], "Print"] }, Open ]], Cell[TextData[{ StyleBox["Cell 3. ", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, Evaluatable->False, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, Background->RGBColor[0, 1, 1]], StyleBox["IncVelocity ", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, Evaluatable->False, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, FontWeight->"Bold", Background->RGBColor[0, 1, 1]], StyleBox[ "computes of incremental velocity vector v at a given state (lambda,u). On \ entry\n K and q are evaluated. If K is numerically singular u is perturbed \ by a tiny random amount and\n the process repeated up to 10 times. If K is \ not numerically singular, the velocity v=Kinv.q is returned.\n \n ", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, Evaluatable->False, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, Background->RGBColor[0, 1, 1]], StyleBox["ResVelocity ", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, Evaluatable->False, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, FontWeight->"Bold", Background->RGBColor[0, 1, 1]], StyleBox["is similar to ", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, Evaluatable->False, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, Background->RGBColor[0, 1, 1]], StyleBox["IncVelocity", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, Evaluatable->False, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, FontWeight->"Bold", Background->RGBColor[0, 1, 1]], StyleBox[ " but it computes and returns the residual r and the\n residual velocity \ dr=-Kinv.r . It supports a correctrive process.", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, Evaluatable->False, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, Background->RGBColor[0, 1, 1]] }], "Text", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, Evaluatable->False, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, Background->RGBColor[0, 1, 1]], Cell[CellGroupData[{ Cell[TextData[{ StyleBox[ "IncVelocity[sprop_,force_,{lambda_,u_}]:= \n \ Module[{K,q,nudge,un=u,v,OK,ueps=10.^(-10)},\n nudge=0; \ SeedRandom[7654321];\n While [nudge<=10,\n \ K=FormTanStiffnessOfTwoBarArch[sprop,un]; K=N[K];\n ", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, InitializationCell->True, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, Background->RGBColor[1, 0.647715, 0.521981]], StyleBox[" (*Print[\"IncVelocity K=\",K];*)", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, InitializationCell->True, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, FontColor->RGBColor[1, 0, 0], Background->RGBColor[1, 0.647715, 0.521981]], StyleBox[ "\n q=FormIncrementalLoadVector[lambda,force,un]; q=N[q];\n \ {v,OK}=SolveTwoLinearEqs[K,q]; \n ", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, InitializationCell->True, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, Background->RGBColor[1, 0.647715, 0.521981]], StyleBox[" (*Print[\"IncVelocity v=\",v,\" OK= \",OK];*)\n ", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, InitializationCell->True, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, FontColor->RGBColor[1, 0, 0], Background->RGBColor[1, 0.647715, 0.521981]], StyleBox[ "If [OK, Break[]];\n nudge++; un+=Table[Random[]*ueps,{2}] ];\n If \ [OK, Return [{N[v],N[q],\" \"}]];\n Return [Null,Null,\"IncVelocity: Cant \ escape singularity\"]];\n\nResVelocity[sprop_,force_,{lambda_,u_}]:= \n \ Module[{K,r,p,nudge,un=u,dr,OK,ueps=10.^(-10)},\n nudge=0; \ SeedRandom[7654321];\n ", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, InitializationCell->True, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, Background->RGBColor[1, 0.647715, 0.521981]], StyleBox[" (*Print[\"ResVelocity, lambda=\",lambda,\" un=\",un];*)", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, InitializationCell->True, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, FontColor->RGBColor[1, 0, 0], Background->RGBColor[1, 0.647715, 0.521981]], StyleBox[ "\n While [nudge<=10,\n K=FormTanStiffnessOfTwoBarArch[sprop,un]; \ K=N[K];\n p=FormInternalForceOfTwoBarArch[sprop,un]; p=N[p];\n \ r=p-lambda*force; r=N[r];\n", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, InitializationCell->True, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, Background->RGBColor[1, 0.647715, 0.521981]], StyleBox[" (*Print[\"K=\",K,\" p=\",p,\" r=\",r];*)", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, InitializationCell->True, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, FontColor->RGBColor[1, 0, 0], Background->RGBColor[1, 0.647715, 0.521981]], StyleBox[" \n {dr,OK}=SolveTwoLinearEqs[K,-r]; \n ", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, InitializationCell->True, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, Background->RGBColor[1, 0.647715, 0.521981]], StyleBox[" (*Print[\"ResVelocity dr=\",dr,\" OK= \",OK];*)\n ", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, InitializationCell->True, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, FontColor->RGBColor[1, 0, 0], Background->RGBColor[1, 0.647715, 0.521981]], StyleBox[ "If [OK, Break[]];\n nudge++; un+=Table[Random[]*ueps,{2}] ];\n If \ [OK, Return [{N[dr],N[r],\" \"}]];\n Return [Null,Null,\"ResVelocity: Cant \ escape singularity\"]];\n ", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, InitializationCell->True, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, Background->RGBColor[1, 0.647715, 0.521981]], StyleBox[ "\nClearAll[uX,uY];\n{v,q,status}= IncVelocity[{2,3,1,1},{4,0},{1,{1,1}}];\n\ Print[\"v=\",v]; Print[\"q=\",q]; \nIf[status!=\" \", \ Print[\"status=\",status]];\n{dr,r,status}= \ ResVelocity[{2,3,1,1},{0,-0.33},{1,{0,-1}}];\nPrint[\"dr=\",dr]; \ Print[\"r=\",r]; \nIf[status!=\" \", Print[\"status=\",status]];", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, InitializationCell->True, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, FontColor->RGBColor[1, 0, 0], Background->RGBColor[1, 0.647715, 0.521981]] }], "Input", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, InitializationCell->True, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, Background->RGBColor[1, 0.647715, 0.521981]], Cell[BoxData[ InterpretationBox[ \("v="\[InvisibleSpace]{12.1626063852629973`, \(-2.4325212770525999`\)} \), SequenceForm[ "v=", {12.162606385262997, -2.4325212770525999}], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[\("q="\[InvisibleSpace]{4.`, 0}\), SequenceForm[ "q=", {4.0, 0}], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("dr="\[InvisibleSpace]{0, \(-0.14517209285188402`\)}\), SequenceForm[ "dr=", {0, -0.14517209285188401}], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[\("r="\[InvisibleSpace]{0, 0.0137722339831620743`}\), SequenceForm[ "r=", {0, 0.013772233983162074}], Editable->False]], "Print"] }, Open ]], Cell[TextData[{ StyleBox["Cell 4. ", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, Evaluatable->False, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, Background->RGBColor[0, 1, 1]], StyleBox["IncSolution ", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, Evaluatable->False, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, FontWeight->"Bold", Background->RGBColor[0, 1, 1]], StyleBox[ "advances the solution over one increment, without executing\nequilibrium \ corrections. The module arguments provide the following input data:\n \n \ method Keywords specifying the solution method: {integ,ics} \ where \n integ identifies the integrator: \ FE, MR or RK4\n ics identifies the \ increment control strategy: LC, DC or AC\n sprop \ Properties of arch structure: {H,S,Em,A0}. See Cell 1.\n force \ A list specifying the applied forces. In this problem, force is \ same as q.\n sol A list of entities at the last \ computed solution:\n n \ current increment number (n=0 for first step)\n \ lambdan the stage control parameter \n \ un ={uXn,uYn} the crown displacement components\n \ vn={vXn,vYn} incremental velocity over \ previous step\n qn={qXn,qYn} \ incremental force over previous step\n \ Kdetn the determinant of Kn\n \ elln same as ell for a fixed stepsize incrementation\n \ an cosine of angle \ between last 2 incremental directions\n \ kappan the limit point sensor\n \ kappa0 kappa at n=0\n \ rem a char string storing remarks about any remarkable \ happenings\n solpar A list of solution control \ parameters: \n nmax maximum \ number of incremental steps\n lambdamax \ bound on abs(lambda)\n umax \ {uXmax,uYmax} bounds on node displacements\n \ ell dimensionless length of increment\n \ ellmin not used here\n \ ellmax not used here\n \ acctol not used here\n \ lref a reference length for scaling u if ics is DC or \ AC\n\n If n>nmax, abs(lambda)>lambdamax, \ abs(uX)>uXmax or\n abs(uY)>uYmax, the \ solution is terminated. If a variable stepsize\n \ implementation is adopted, ellmin, elklmax and epsacc are used\n \ In corrective and perturbation-injection method, \ this list is\n expanded.\n \ \n The function returns {solnext,status} where\n \ solnext Updated values of sol at the next solution if status \ is blank\n status A character string status \ indicator. If nonblank, the solution process has been\n \ terminated for the reason stated in this string.\n \ \n In essence ", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, Evaluatable->False, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, Background->RGBColor[0, 1, 1]], StyleBox["IncSolution ", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, Evaluatable->False, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, FontWeight->"Bold", Background->RGBColor[0, 1, 1]], StyleBox[ "is a driver that checks whether certain bounds are exceeded (for example\n\ number of steps) and if satisfied simply calls the appropriate integrator. \ If termination is\n indicated by bound violation(s) or other reasons, it \ returns Null; else it reports the next solution.\n \n ", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, Evaluatable->False, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, Background->RGBColor[0, 1, 1]], StyleBox["IncIterSolution ", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, Evaluatable->False, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, FontWeight->"Bold", Background->RGBColor[0, 1, 1]], StyleBox[ "is a driver for an incremental-iterative solution method. After executing\ \n an increment step like", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, Evaluatable->False, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, Background->RGBColor[0, 1, 1]], StyleBox[" IncSolution", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, Evaluatable->False, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, FontWeight->"Bold", Background->RGBColor[0, 1, 1]], StyleBox[ ", it does a Conventional Newton iteration process via\n module ", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, Evaluatable->False, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, Background->RGBColor[0, 1, 1]], StyleBox["NewtonIteration", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, Evaluatable->False, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, FontWeight->"Bold", Background->RGBColor[0, 1, 1]], StyleBox[".", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, Evaluatable->False, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, Background->RGBColor[0, 1, 1]] }], "Text", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, Evaluatable->False, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, Background->RGBColor[0, 1, 1]], Cell[TextData[{ StyleBox[ "IncSolution[method_,sprop_,force_,sol_,solpar_]:=\n Module[ \ {integ,ics,n=sol[[1]],lambda=sol[[2]],u=sol[[3]],\n \ nmax=solpar[[1]],lambdamax=solpar[[2]],umax=solpar[[3]],\n \ solnext,status},\n If [n>nmax, Return[{Null,\n \"IncSolution: Max \ steps reached\"}]];\n If [Abs[lambda]>Abs[lambdamax], Return[{Null,\n \ \"IncSolution: lambda exceeds bound\"}]];\n If [Abs[u[[1]]]>Abs[umax[[1]]], \ Return[{Null,\n \"IncSolution: uX exceeds bound\"}]];\n If \ [Abs[u[[2]]]>Abs[umax[[2]]], Return[{Null,\n \"IncSolution: uY exceeds \ bound\"}]]; \n {integ,ics}=method;\n If [integ==\"FE\", \n \ {solnext,status}=FEstep [ics,sprop,force,sol,solpar]];\n If [integ==\"MR\", \ \n {solnext,status}=MRstep [ics,sprop,force,sol,solpar]];\n If \ [integ==\"RK4\", \n \ {solnext,status}=RK4step[ics,sprop,force,sol,solpar]];\n \ Return[{solnext,status}]];\n \n\ IncIterSolution[method_,sprop_,force_,sol_,solpar_]:=\n Module[ \ {integ,ics,n=sol[[1]],lambda=sol[[2]],u=sol[[3]],\n \ nmax=solpar[[1]],lambdamax=solpar[[2]],umax=solpar[[3]],\n \ solpred,solnext,status,sol0}, \n If [n>nmax, Return[{Null,\n \ \"IncSolution: Max steps reached\"}]];\n If [Abs[lambda]>Abs[lambdamax], \ Return[{Null,\n \"IncSolution: lambda exceeds bound\"}]];\n If \ [Abs[u[[1]]]>Abs[umax[[1]]], Return[{Null,\n \"IncSolution: uX exceeds \ bound\"}]];\n If [Abs[u[[2]]]>Abs[umax[[2]]], Return[{Null,\n \ \"IncSolution: uY exceeds bound\"}]]; \n {integ,ics}=method; sol0=sol;", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, InitializationCell->True, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, Background->RGBColor[1, 0.647715, 0.521981]], StyleBox[" ", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, InitializationCell->True, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, FontColor->RGBColor[1, 0, 0], Background->RGBColor[1, 0.647715, 0.521981]], StyleBox[ "\n If [integ==\"FE\", \n {solpred,status}=FEstep \ [ics,sprop,force,sol,solpar]];\n If [integ==\"MR\", \n \ {solpred,status}=MRstep [ics,sprop,force,sol,solpar]];\n If \ [integ==\"RK4\", \n \ {solpred,status}=RK4step[ics,sprop,force,sol,solpar]];\n \n \ {solnext,status}=NewtonIteration[ics,sprop,force,solpred,\n \ solpar,sol0];\n Return[{solnext,status}]];", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, InitializationCell->True, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, Background->RGBColor[1, 0.647715, 0.521981]] }], "Input", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, InitializationCell->True, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, Background->RGBColor[1, 0.647715, 0.521981]], Cell[TextData[{ StyleBox["Cell 5. ", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, Evaluatable->False, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, Background->RGBColor[0, 1, 1]], StyleBox["FEstep", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, Evaluatable->False, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, FontWeight->"Bold", Background->RGBColor[0, 1, 1]], StyleBox[ " is the Forward Euler (FE) driver module. It simply calls FELCstep, \ FEDCstep or\nFEACstep according to the increment control strategy specified \ in ics.", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, Evaluatable->False, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, Background->RGBColor[0, 1, 1]] }], "Text", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, Evaluatable->False, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, Background->RGBColor[0, 1, 1]], Cell["\<\ FEstep[ics_,sprop_,force_,sol_,solpar_]:=Module[ {solnext,status}, If [ics==\"LC\", {solnext,status}=FELCstep[ics,sprop,force,sol,solpar]]; If [ics==\"DC\", {solnext,status}=FEDCstep[ics,sprop,force,sol,solpar]]; If [ics==\"AC\", {solnext,status}=FEACstep[ics,sprop,force,sol,solpar]]; Return[{solnext,status}]];\ \>", "Input", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, InitializationCell->True, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, Background->RGBColor[1, 0.647715, 0.521981]], Cell[TextData[{ StyleBox["Cell 6. ", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, Evaluatable->False, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, Background->RGBColor[0, 1, 1]], StyleBox["FELCstep, FEDCstep ", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, Evaluatable->False, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, FontWeight->"Bold", Background->RGBColor[0, 1, 1]], StyleBox["and", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, Evaluatable->False, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, Background->RGBColor[0, 1, 1]], StyleBox[" FEACstep ", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, Evaluatable->False, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, FontWeight->"Bold", Background->RGBColor[0, 1, 1]], StyleBox[ " implement the Forward Euler (FE) integrator for Load\nControl, \ Displacement Control and ArcLength Control, respectively.", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, Evaluatable->False, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, Background->RGBColor[0, 1, 1]] }], "Text", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, Evaluatable->False, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, Background->RGBColor[0, 1, 1]], Cell["\<\ FELCstep[ics_,sprop_,force_,sol_,solpar_]:=Module[ {n,lambdan,un,vn,qn,Kdetn,elln,an,kappan,kappa0,rn, rem=\" \",v,q,qv,lambda,u,Kdet,a,kappa,nmax,lambdamax,umax, ell,ellmin,ellmax,acctol,lref,f,status}, {n,lambdan,un,vn,qn,Kdetn,elln,an,kappan,kappa0,rn}=sol; {nmax,lambdamax,umax,ell,ellmin,ellmax,acctol,lref}=solpar; {v,q,status}=IncVelocity[sprop,force,{lambdan,un}]; If [status!=\" \", rem=status; Return[{Null,status}]]; n++; qv=q.v; f=Sign[qv]; ell=elln; lambda=lambdan+ell/f; u=un+v*ell/f; Kdet=DetTanStiffness[sprop,u]; a=(v/Sqrt[v.v]).(vn/Sqrt[vn.vn]); kappa=((q.v)/(v.v))/kappa0; Return[{{n,lambda,u,v,q,Kdet,ell,a,kappa,kappa0,rem},\" \"}]]; FEDCstep[ics_,sprop_,force_,sol_,solpar_]:=Module[ {n,lambdan,un,vn,qn,Kdetn,elln,an,kappan,kappa0,rn, rem=\" \",v,q,qv,lambda,u,Kdet,a,kappa,nmax,lambdamax,umax, ell,ellmin,ellmax,acctol,lref,f,status}, {n,lambdan,un,vn,qn,Kdetn,elln,an,kappan,kappa0,rn}=sol; {nmax,lambdamax,umax,ell,ellmin,ellmax,acctol,lref}=solpar; {v,q,status}=IncVelocity[sprop,force,{lambdan,un}]; If [status!=\" \", rem=status; Return[{Null,status}]]; n++; qv=q.v; f=Sqrt[(v.v)/lref^2]*Sign[qv]; ell=elln; lambda=lambdan+ell/f; u=un+v*ell/f; Kdet=DetTanStiffness[sprop,u]; a=(v/Sqrt[v.v]).(vn/Sqrt[vn.vn]); kappa=((q.v)/(v.v))/kappa0; Return[{{n,lambda,u,v,q,Kdet,ell,a,kappa,kappa0,rem},\" \"}]]; FEACstep[ics_,sprop_,force_,sol_,solpar_]:=Module[ {n,lambdan,un,vn,qn,Kdetn,elln,an,kappan,kappa0,rn, rem=\" \",v,q,qv,lambda,u,Kdet,a,kappa,nmax,lambdamax,umax, ell,ellmin,ellmax,acctol,lref,f,status}, {n,lambdan,un,vn,qn,Kdetn,elln,an,kappan,kappa0,rn}=sol; {nmax,lambdamax,umax,ell,ellmin,ellmax,acctol,lref}=solpar; {v,q,status}=IncVelocity[sprop,force,{lambdan,un}]; If [status!=\" \", rem=status; Return[{Null,status}]]; n++; qv=q.v; f=Sqrt[1+(v.v)/lref^2]*Sign[qv]; ell=elln; lambda=lambdan+ell/f; u=un+v*ell/f; Kdet=DetTanStiffness[sprop,u]; a=(v/Sqrt[v.v]).(vn/Sqrt[vn.vn]); kappa=((q.v)/(v.v))/kappa0; Return[{{n,lambda,u,v,q,Kdet,ell,a,kappa,kappa0,rem},\" \"}]];\ \>", "Input", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, InitializationCell->True, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, Background->RGBColor[1, 0.647715, 0.521981]], Cell[TextData[{ StyleBox["Cell 7. ", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, Evaluatable->False, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, Background->RGBColor[0, 1, 1]], StyleBox[" MRstep", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, Evaluatable->False, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, FontWeight->"Bold", Background->RGBColor[0, 1, 1]], StyleBox[ " is the Midpoint Rule (MR) driver module. It simply calls MRLCstep, \ MRDCstep or\nMRACstep according to the increment control strategy specified \ in ics.", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, Evaluatable->False, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, Background->RGBColor[0, 1, 1]] }], "Text", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, Evaluatable->False, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, Background->RGBColor[0, 1, 1]], Cell["\<\ MRstep[ics_,sprop_,force_,sol_,solpar_]:=Module[ {solnext,status}, If [ics==\"LC\", {solnext,status}=MRLCstep[ics,sprop,force,sol,solpar]]; If [ics==\"DC\", {solnext,status}=MRDCstep[ics,sprop,force,sol,solpar]]; If [ics==\"AC\", {solnext,status}=MRACstep[ics,sprop,force,sol,solpar]]; Return[{solnext,status}]];\ \>", "Input", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, InitializationCell->True, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, Background->RGBColor[1, 0.647715, 0.521981]], Cell[TextData[{ StyleBox["Cell 8.", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, FontFamily->"Times", FontWeight->"Plain", Background->RGBColor[0, 1, 1]], StyleBox[" MRLCstep, MRDCstep ", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, FontFamily->"Times", Background->RGBColor[0, 1, 1]], StyleBox["and ", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, FontFamily->"Times", FontWeight->"Plain", Background->RGBColor[0, 1, 1]], StyleBox["MRACstep ", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, FontFamily->"Times", Background->RGBColor[0, 1, 1]], StyleBox[ " implement the Midpoint Rule (MR) integrator for Load\nControl, \ Displacement Control and ArcLength Control, respectively.", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, FontFamily->"Times", FontWeight->"Plain", Background->RGBColor[0, 1, 1]] }], "Input", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, Background->RGBColor[0, 1, 1]], Cell["\<\ MRLCstep[ics_,sprop_,force_,sol_,solpar_]:=Module[ {n,lambdan,un,vn,qn,Kdetn,elln,an,kappan,kappa0,rn, rem=\" \",v,q,qv,lambda,u,Kdet,a,kappa,nmax,lambdamax,umax, ell,ellmin,ellmax,acctol,lref,f,status}, {n,lambdan,un,vn,qn,Kdetn,elln,an,kappan,kappa0,rn}=sol; {nmax,lambdamax,umax,ell,ellmin,ellmax,acctol,lref}=solpar; {v,q,status}=IncVelocity[sprop,force,{lambdan,un}]; If [status!=\" \", Return[{Null,status}]]; qv=q.v; f=Sign[qv]; ell=elln; lambdaM=lambdan+ell/(2*f); uM=un+v*ell/(2*f); {v,q,status}=IncVelocity[sprop,force,{lambdaM,uM}]; If [status!=\" \", rem=status; Return[{Null,status}]]; n++; qv=q.v; f=Sign[qv]; ell=elln; lambda=lambdan+ell/f; u=un+v*ell/f; Kdet=DetTanStiffness[sprop,u]; a=(v/Sqrt[v.v]).(vn/Sqrt[vn.vn]); kappa=((q.v)/(v.v))/kappa0; Return[{{n,lambda,u,v,q,Kdet,ell,a,kappa,kappa0,rem},\" \"}]]; MRDCstep[ics_,sprop_,force_,sol_,solpar_]:=Module[ {n,lambdan,un,vn,qn,Kdetn,elln,an,kappan,kappa0,rn, rem=\" \",v,q,qv,lambda,u,Kdet,a,kappa,nmax,lambdamax,umax, ell,ellmin,ellmax,acctol,lref,f,status}, {n,lambdan,un,vn,qn,Kdetn,elln,an,kappan,kappa0,rem}=sol; {nmax,lambdamax,umax,ell,ellmin,ellmax,acctol,lref}=solpar; {v,q,status}=IncVelocity[sprop,force,{lambdan,un}]; If [status!=\" \", Return[{Null,status}]]; qv=q.v; f=Sqrt[(v.v)/lref^2]*Sign[qv]; ell=elln; lambdaM=lambdan+ell/(2*f); uM=un+v*ell/(2*f); {v,q,status}=IncVelocity[sprop,force,{lambdaM,uM}]; If [status!=\" \", rem=status; Return[{Null,status}]]; n++; qv=q.v; f=Sqrt[(v.v)/lref^2]*Sign[qv]; ell=elln; lambda=lambdan+ell/f; u=un+v*ell/f; Kdet=DetTanStiffness[sprop,u]; a=(v/Sqrt[v.v]).(vn/Sqrt[vn.vn]); kappa=((q.v)/(v.v))/kappa0; Return[{{n,lambda,u,v,q,Kdet,ell,a,kappa,kappa0,rem},\" \"}]]; MRACstep[ics_,sprop_,force_,sol_,solpar_]:=Module[ {n,lambdan,un,vn,qn,Kdetn,elln,an,kappan,kappa0,rn, rem=\" \",v,q,qv,lambda,u,Kdet,a,kappa,nmax,lambdamax,umax, ell,ellmin,ellmax,acctol,lref,f,status}, {n,lambdan,un,vn,qn,Kdetn,elln,an,kappan,kappa0,rn}=sol; {nmax,lambdamax,umax,ell,ellmin,ellmax,acctol,lref}=solpar; {v,q,status}=IncVelocity[sprop,force,{lambdan,un}]; If [status!=\" \", Return[{Null,status}]]; qv=q.v; f=Sqrt[1+(v.v)/lref^2]*Sign[qv]; ell=elln; lambdaM=lambdan+ell/(2*f); uM=un+v*ell/(2*f); {v,q,status}=IncVelocity[sprop,force,{lambdaM,uM}]; If [status!=\" \", rem=status; Return[{Null,status}]]; n++; qv=q.v; f=Sqrt[1+(v.v)/lref^2]*Sign[qv]; ell=elln; lambda=lambdan+ell/f; u=un+v*ell/f; Kdet=DetTanStiffness[sprop,u]; a=(v/Sqrt[v.v]).(vn/Sqrt[vn.vn]); kappa=((q.v)/(v.v))/kappa0; Return[{{n,lambda,u,v,q,Kdet,ell,a,kappa,kappa0,rem},\" \"}]]; \ \>", "Input", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, InitializationCell->True, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, Background->RGBColor[1, 0.647715, 0.521981]], Cell[TextData[{ StyleBox["Cell 9. ", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, Evaluatable->False, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, Background->RGBColor[1, 0, 0]], StyleBox[" RK4step", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, Evaluatable->False, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, FontWeight->"Bold", Background->RGBColor[1, 0, 0]], StyleBox[ " is the 4th order Runge Kutta (RK4) driver module. It simply calls \ RK4LCstep, \nRK4DCstep or RK4ACstep according to the increment control \ strategy specified in ics.", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, Evaluatable->False, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, Background->RGBColor[1, 0, 0]] }], "Text", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, Evaluatable->False, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, Background->RGBColor[1, 0, 0]], Cell["\<\ RK4step[ics_,sprop_,force_,sol_,solpar_]:=Module[ {solnext,status}, If [ics==\"LC\", {solnext,status}=RK4LCstep[ics,sprop,force,sol,solpar]]; If [ics==\"DC\", {solnext,status}=RK4DCstep[ics,sprop,force,sol,solpar]]; If [ics==\"AC\", {solnext,status}=RK4ACstep[ics,sprop,force,sol,solpar]]; If [status!=\" \", Return[{Null,status}]]; Return[{solnext,\" \"}]];\ \>", "Input", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, InitializationCell->True, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, Background->RGBColor[1, 0.647715, 0.521981]], Cell[TextData[{ StyleBox["Cell 10. ", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, FontFamily->"Times", FontWeight->"Plain", Background->RGBColor[1, 0, 0]], StyleBox["RK4LCstep, RK4DCstep ", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, FontFamily->"Times", Background->RGBColor[1, 0, 0]], StyleBox["and ", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, FontFamily->"Times", FontWeight->"Plain", Background->RGBColor[1, 0, 0]], StyleBox["RK4ACstep ", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, FontFamily->"Times", Background->RGBColor[1, 0, 0]], StyleBox[ " implement the classical Runge-Kutta (RK4) integrator\n for Load Control, \ Displacement Control and ArcLength Control, respectively. ", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, FontFamily->"Times", FontWeight->"Plain", Background->RGBColor[1, 0, 0]] }], "Input", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, Background->RGBColor[1, 0, 0]], Cell["\<\ RK4ACstep[ics_,sprop_,force_,sol_,solpar_]:=Module[ {n,lambdan,un,vn,qn,Kdetn,elln,an,kappan,kappa0,rn, rem=\" \",v,v1,v2,v3,v4,q,qv,lambda,u,Kdet,a,kappa,nmax, lambdamax,umax, ell,ellmin,ellmax,acctol,lref,f,status}, {n,lambdan,un,vn,qn,Kdetn,elln,an,kappan,kappa0,rn}=sol; {nmax,lambdamax,umax,ell,ellmin,ellmax,acctol,lref}=solpar; {v1,q,status}=IncVelocity[sprop,force,{lambdan,un}]; If [status!=\" \", Return[{Null,status}]]; qv=q.v1; f=Sqrt[(v1.v1)/lref^2]*Sign[qv]; ell=elln; lambdaM=lambdan+ell/(2*f); uM=un+v1*ell/(2*f); {v2,q,status}=IncVelocity[sprop,force,{lambdaM,uM}]; If [status!=\" \", Return[{Null,status}]]; qv=q.v2; f=Sqrt[(v2.v2)/lref^2]*Sign[qv]; ell=elln; lambdaM=lambdan+ell/(2*f); uM=un+v2*ell/(2*f); {v3,q,status}=IncVelocity[sprop,force,{lambdaM,uM}]; If [status!=\" \", Return[{Null,status}]]; qv=q.v3; f=Sqrt[(v3.v3)/lref^2]*Sign[qv]; ell=elln; lambdaM=lambdan+ell/(f); uM=un+v3*ell/(f); {v4,q,status}=IncVelocity[sprop,force,{lambdaM,uM}]; If [status!=\" \", rem=status; Return[{Null,status}]]; v=(v1+2*v2+2*v3+v4)/6; n++; qv=q.v; f=Sqrt[(v.v)/lref^2]*Sign[qv]; ell=elln; lambda=lambdan+ell/f; u=un+v*ell/f; Kdet=DetTanStiffness[sprop,u]; a=(v/Sqrt[v.v]).(vn/Sqrt[vn.vn]); kappa=((q.v)/(v.v))/kappa0; Return[{{n,lambda,u,v,q,Kdet,ell,a,kappa,kappa0,rem},\" \"}]]; \ \>", "Input", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, InitializationCell->True, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, Background->RGBColor[1, 0.647715, 0.521981]], Cell[TextData[ "Cell 11. Module to print computed solution table. Adapted from Fortran."], "Text", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, Evaluatable->False, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, Background->RGBColor[0, 1, 1]], Cell[CellGroupData[{ Cell["\<\ PrintSolutionTable[soltab_]:= Module[{numsteps, n,lambdan,un,vn,qn,Kdetn,elln,an,kappan,kappa0,rem,t}, numsteps=Length[soltab]; t=Table[\" \",{numsteps+1},{8}]; For [i=1, i<=numsteps, i++, {n,lambdan,un,vn,qn,Kdetn,elln,an,kappan,kappa0,rem}= Chop[soltab[[i]],.00000001]; t[[i+1,1]]=ToString[n]; t[[i+1,2]]=ToString[lambdan]; t[[i+1,3]]=ToString[un]; t[[i+1,4]]=ToString[vn]; t[[i+1,5]]=ToString[Kdetn]; t[[i+1,6]]=ToString[an]; t[[i+1,7]]=ToString[kappan]; t[[i+1,8]]=rem]; t[[1]] = {\"n\",\"lambda\",\"u={uX,uY}\",\"v={vX,vY}\", \"Kdet\", \"a\", \"kappa\", \"remarks\"}; Print[TableForm[t,TableAlignments->{Left}, TableDirections->{Column,Row},TableSpacing->{0,2}]]; ]; q={0,-1}; soltab= { {0,0.0,{ 0.0,0.},{0.42,0.0},q,45.,0.02,1.0,1.0,2.5,\"ref state\"}, {1,0.2,{-0.3,0.},{0.35,0.0},q,43.,0.02,1.0,0.9,2.5,\"step 1\"} }; PrintSolutionTable[soltab];\ \>", "Input", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, InitializationCell->True, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, Background->RGBColor[1, 0.572488, 0.43003]], Cell[BoxData[ TagBox[GridBox[{ {"\<\"n\"\>", "\<\"lambda\"\>", "\<\"u={uX,uY}\"\>", "\<\"v={vX,vY}\"\>", "\<\"Kdet\"\>", "\<\"a\"\>", "\<\"kappa\"\>", "\<\"remarks\"\>"}, {"\<\"0\"\>", "\<\"0\"\>", "\<\"{0, 0}\"\>", "\<\"{0.42, 0}\"\>", "\<\"45.\"\>", "\<\"1.\"\>", "\<\"1.\"\>", "\<\"ref state\"\>"}, {"\<\"1\"\>", "\<\"0.2\"\>", "\<\"{-0.3, 0}\"\>", "\<\"{0.35, 0}\"\>", "\<\"43.\"\>", "\<\"1.\"\>", "\<\"0.9\"\>", "\<\"step 1\"\>"} }, RowSpacings->0, ColumnSpacings->2, RowAlignments->Baseline, ColumnAlignments->{Left}], (TableForm[ #, TableAlignments -> {Left}, TableDirections -> {Column, Row}, TableSpacing -> {0, 2}]&)]], "Print"] }, Open ]], Cell[TextData[ "Cell 12. Driver of Newton iteration for corrective process. The following \ solution parameters are\nbuilt in: \n\n mNiter max \ number of Newton iterations, set to 8\n epsconv \ convergence tolerance for residual norm ratio, set to 0.00001\n \ facdiv divergence factor for residual norm ratio, set to 10000."], "Text", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, Evaluatable->False, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, Background->RGBColor[0, 1, 1]], Cell[TextData[{ StyleBox[ "NewtonIteration[ics_,sprop_,force_,sol_,solpar_,sol0_]:=Module[\n \ {n,lambdan,un,vn,qn,Kdetn,elln,an,kappan,kappa0,remn,\n rem=\" \ \",state,mNiter,epsconv,facdiv,r0,r0norm,rk,rknorm,\n \ statenext,converged,status,lambda0,u0},\n \ {n,lambdan,un,vn,qn,Kdetn,elln,an,kappan,kappa0,remn}=sol;\n mNiter=2; \ epsconv=0.00001; facdiv=10000.;\n \ p=FormInternalForceOfTwoBarArch[sprop,un]; p=N[p];\n", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, InitializationCell->True, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, Background->RGBColor[1, 0.572488, 0.43003]], StyleBox[ " (*Print[\"NewtonIteration entered\"];\n Print[\"sprop=\",sprop,\" \ un=\",un, \" p=\",p]; *)", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, InitializationCell->True, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, FontColor->RGBColor[1, 0, 0], Background->RGBColor[1, 0.572488, 0.43003]], StyleBox["\n rn=p-lambdan*force; rn=N[rn];\n rnnorm=Sqrt[rn.rn];", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, InitializationCell->True, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, Background->RGBColor[1, 0.572488, 0.43003]], StyleBox[" (*Print[\"rn=\",rn,\" rnnorm=\",rnnorm];*)", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, InitializationCell->True, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, FontColor->RGBColor[1, 0, 0], Background->RGBColor[1, 0.572488, 0.43003]], StyleBox[ "\n state={lambdan,un,rn}; lambda0=sol0[[2]]; u0=sol0[[3]];\n If \ [rnnorm==0, Return[sol,status]]; converged=False;\n For [k=0, k<=mNiter, \ k++, \n If [ics==\"LC\", \n \ {statenext,status}=LCNewtonCycle[sprop,force,state]];\n If [ics==\"DC\", \ \n {statenext,status}=DCNewtonCycle[sprop,force,state,\n \ {lambda0,u0,elln}]];\n If [ics==\"AC\", \n \ {statenext,status}=ACNewtonCycle[sprop,force,state,\n \ {lambda0,u0,v0,elln} ]];\n If [status!=\" \",Break[]];\n \ {lambda,u,r}=statenext; rnorm=Sqrt[r.r];\n ", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, InitializationCell->True, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, Background->RGBColor[1, 0.572488, 0.43003]], StyleBox[ "Print[\"lambda=\",lambda,\" u=\",u,\" r=\",r,\n \" \ rnorm=\",rnorm//InputForm];\n (*Print[\"norm ratio=\",rnorm/rnnorm,\" vs \ \",epsconv];*)", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, InitializationCell->True, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, FontColor->RGBColor[1, 0, 0], Background->RGBColor[1, 0.572488, 0.43003]], StyleBox[ "\n If [rnorm/rnnormfacdiv, status=\"Newton iteration diverges\";\n \ Return[Null,status]];\n state=statenext;\n ];\n If [!converged, \ status=\"Newton didnt converge in max iters\"];\n If [converged, \ Print[\"Newton converged after \",k,\" cycles\"]];\n {lambda,u,r}=state;\n \ solnext={n,lambda,u,vn,qn,Kdetn,elln,an,kappan,kappa0,remn};\n \ Return[{solnext,status}]];\n ", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, InitializationCell->True, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, Background->RGBColor[1, 0.572488, 0.43003]] }], "Input", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, InitializationCell->True, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, Background->RGBColor[1, 0.572488, 0.43003]], Cell[TextData[{ StyleBox[ "Cell 13. These modules perform one Newton iteration cycle under different \ increment \ncontrol constraints. They are called by ", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, Evaluatable->False, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, Background->RGBColor[0, 1, 1]], StyleBox["NewtonIteration.", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, Evaluatable->False, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, FontWeight->"Bold", Background->RGBColor[0, 1, 1]], StyleBox[ "\nOnly the Load Control (LC) and Displacement Control (DC)\nare \ implemented in modules ", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, Evaluatable->False, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, Background->RGBColor[0, 1, 1]], StyleBox["LCNewtonCycle", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, Evaluatable->False, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, FontWeight->"Bold", Background->RGBColor[0, 1, 1]], StyleBox[" and", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, Evaluatable->False, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, Background->RGBColor[0, 1, 1]], StyleBox[" DCNewtonCycle.", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, Evaluatable->False, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, FontWeight->"Bold", Background->RGBColor[0, 1, 1]], StyleBox["\nThe arguments of ", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, Evaluatable->False, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, Background->RGBColor[0, 1, 1]], StyleBox["LCNewtonCycle", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, Evaluatable->False, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, FontWeight->"Bold", Background->RGBColor[0, 1, 1]], StyleBox[ " are\n sprop structural properties\n \ force reference force\n {lambdak,uk,rk} \ lambda, u and r from the previous (kth) iteration cycle \n It returns \ {{lambda,u,r},\" \"} if no error detected: {lambda,u,r} are the updated \ values;\n otherwise {Null,Status} where status has an error message\n \n \ The arguments of", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, Evaluatable->False, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, Background->RGBColor[0, 1, 1]], StyleBox[" DCNewtonCycle", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, Evaluatable->False, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, FontWeight->"Bold", Background->RGBColor[0, 1, 1]], StyleBox[ " are\n sprop structural properties\n \ force reference force\n {lambdak,uk,rk} \ lambda, u, and r from the previous (kth) iteration cycle \n \ {lambda0,u0,ell} information needed to specify the increment control \ constraint c=0.\n {lambda0,u0} is the \ previous converged solution, ell the increment length\n It returns the same \ information as LCNewtonCycle", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, Evaluatable->False, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, Background->RGBColor[0, 1, 1]] }], "Text", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, Evaluatable->False, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, Background->RGBColor[0, 1, 1]], Cell[TextData[{ StyleBox[ "LCNewtonCycle[sprop_,force_,{lambdak_,uk_,rk_}]:=Module[\n \ {v,q,status,dr,r,lambda,u},\n", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, InitializationCell->True, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, Background->RGBColor[0.805646, 0.939986, 0.206806]], StyleBox[ " (*Print[\"Enter LCNewtoncycle\"];\n Print[\"sprop=\",sprop,\" force =\ \",force];\n Print[\"entry lambdak=\",lambdak,\" uk=\",uk,\" \ rk=\",rk];*)", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, InitializationCell->True, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, FontColor->RGBColor[1, 0, 0], Background->RGBColor[0.805646, 0.939986, 0.206806]], StyleBox["\n {v,q,status}= IncVelocity[sprop,force,{lambdak,uk}];\n ", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, InitializationCell->True, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, Background->RGBColor[0.805646, 0.939986, 0.206806]], StyleBox[" (*Print[\"v=\",v,\" q=\",q,\" status=\",status];*)", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, InitializationCell->True, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, FontColor->RGBColor[1, 0, 0], Background->RGBColor[0.805646, 0.939986, 0.206806]], StyleBox["\n {dr,r,status}=ResVelocity[sprop,force,{lambdak,uk}];\n ", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, InitializationCell->True, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, Background->RGBColor[0.805646, 0.939986, 0.206806]], StyleBox[" (*Print[\"dr=\",dr,\" r=\",r,\" status=\",status];*)", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, InitializationCell->True, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, FontColor->RGBColor[1, 0, 0], Background->RGBColor[0.805646, 0.939986, 0.206806]], StyleBox[ "\n If [status!=\" \",Return[Null,status]];\n lambda=lambdak; u=uk+dr;\n\ r=FormInternalForceOfTwoBarArch[sprop,u]-lambda*force; r=N[r];\n ", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, InitializationCell->True, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, Background->RGBColor[0.805646, 0.939986, 0.206806]], StyleBox[" (*Print[\"exit lambda=\",lambda,\" u=\",u,\" r=\",r];*)", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, InitializationCell->True, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, FontColor->RGBColor[1, 0, 0], Background->RGBColor[0.805646, 0.939986, 0.206806]], StyleBox[ "\n Return[{{lambda,u,r},\" \"}]];\n \n\ DCNewtonCycle[sprop_,force_,{lambdak_,uk_,rk_},\n \ {lambda0_,u0_,ell_}]:=Module[ {v,q,status,dr,r,lambda,u,\n eta,ck},\n", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, InitializationCell->True, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, Background->RGBColor[0.805646, 0.939986, 0.206806]], StyleBox[ " (*Print[\"Enter LCNewtoncycle\"];\n Print[\"sprop=\",sprop,\" force =\ \",force];\n Print[\"entry lambdak=\",lambdak,\" uk=\",uk,\" rk=\",rk];\n \ Print[\"entry lambda0=\",lambda0,\" u0=\",u0];*)", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, InitializationCell->True, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, FontColor->RGBColor[1, 0, 0], Background->RGBColor[0.805646, 0.939986, 0.206806]], StyleBox["\n {v,q,status}= IncVelocity[sprop,force,{lambdak,uk}];\n ", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, InitializationCell->True, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, Background->RGBColor[0.805646, 0.939986, 0.206806]], StyleBox[" (*Print[\"v=\",v,\" q=\",q,\" status=\",status];*)", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, InitializationCell->True, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, FontColor->RGBColor[1, 0, 0], Background->RGBColor[0.805646, 0.939986, 0.206806]], StyleBox["\n {dr,r,status}=ResVelocity[sprop,force,{lambdak,uk}];\n ", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, InitializationCell->True, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, Background->RGBColor[0.805646, 0.939986, 0.206806]], StyleBox[" (*Print[\"dr=\",dr,\" r=\",r,\" status=\",status];*)", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, InitializationCell->True, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, FontColor->RGBColor[1, 0, 0], Background->RGBColor[0.805646, 0.939986, 0.206806]], StyleBox[ "\n If [status!=\" \",Return[Null,status]];\n ck=(uk-u0).(uk-u0)-ell^2; \ ", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, InitializationCell->True, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, Background->RGBColor[0.805646, 0.939986, 0.206806]], StyleBox[ "Print[\"ck=\",ck];\n (*Print[\"uk-u0=\",uk-u0,\" \ (uk-u0).dr,v=\",(uk-u0).dr,\n \" (uk-u0).v=\",(uk-u0).v];*)", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, InitializationCell->True, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, FontColor->RGBColor[1, 0, 0], Background->RGBColor[0.805646, 0.939986, 0.206806]], StyleBox["\n eta=-(ck+(uk-u0).dr)/((uk-u0).v); ", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, InitializationCell->True, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, Background->RGBColor[0.805646, 0.939986, 0.206806]], StyleBox[" Print[\"eta=\",eta];", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, InitializationCell->True, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, FontColor->RGBColor[1, 0, 0], Background->RGBColor[0.805646, 0.939986, 0.206806]], StyleBox[ "\n lambda=lambdak+eta; u=uk+dr+eta*v;\n \ r=FormInternalForceOfTwoBarArch[sprop,u]-lambda*force; r=N[r];\n", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, InitializationCell->True, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, Background->RGBColor[0.805646, 0.939986, 0.206806]], StyleBox[" (*Print[\"exit lambda=\",lambda,\" u=\",u,\" r=\",r];*)", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, InitializationCell->True, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, FontColor->RGBColor[1, 0, 0], Background->RGBColor[0.805646, 0.939986, 0.206806]], StyleBox[ "\n Return[{{lambda,u,r},\" \"}]];\n \n\ ACNewtonCycle[sprop_,force_,{lambdak_,uk_,rk_},\n \ {lambda0_,u0_,v0_,ell_}]:=Module[ {v,q,status,dr,r,lambda,u,\n eta,ck,f0}, \ f0=Sqrt[1+v0.v0];\n", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, InitializationCell->True, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, Background->RGBColor[0.805646, 0.939986, 0.206806]], StyleBox[ " Print[\"Enter ACNewtoncycle\"];\n Print[\"sprop=\",sprop,\" force \ =\",force];\n Print[\"entry lambdak=\",lambdak,\" uk=\",uk,\" rk=\",rk];\n\ Print[\"entry lambda0=\",lambda0,\" u0=\",u0,\" f0=\",f0];", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, InitializationCell->True, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, FontColor->RGBColor[1, 0, 0], Background->RGBColor[0.805646, 0.939986, 0.206806]], StyleBox["\n {v,q,status}= IncVelocity[sprop,force,{lambdak,uk}];\n ", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, InitializationCell->True, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, Background->RGBColor[0.805646, 0.939986, 0.206806]], StyleBox[" Print[\"v=\",v,\" q=\",q,\" status=\",status];", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, InitializationCell->True, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, FontColor->RGBColor[1, 0, 0], Background->RGBColor[0.805646, 0.939986, 0.206806]], StyleBox["\n {dr,r,status}=ResVelocity[sprop,force,{lambdak,uk}];\n ", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, InitializationCell->True, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, Background->RGBColor[0.805646, 0.939986, 0.206806]], StyleBox[" Print[\"dr=\",dr,\" r=\",r,\" status=\",status];", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, InitializationCell->True, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, FontColor->RGBColor[1, 0, 0], Background->RGBColor[0.805646, 0.939986, 0.206806]], StyleBox[ "\n If [status!=\" \",Return[Null,status]];\n \ ck=(v0.(uk-u0)+(lambdak-lambda0))/f0-ell; ", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, InitializationCell->True, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, Background->RGBColor[0.805646, 0.939986, 0.206806]], StyleBox["Print[\"ck=\",ck];", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, InitializationCell->True, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, FontColor->RGBColor[1, 0, 0], Background->RGBColor[0.805646, 0.939986, 0.206806]], StyleBox["\n eta=-(ck+(v0/f0).dr)/((1+v0.v)/f0); ", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, InitializationCell->True, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, Background->RGBColor[0.805646, 0.939986, 0.206806]], StyleBox[" Print[\"eta=\",eta];", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, InitializationCell->True, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, FontColor->RGBColor[1, 0, 0], Background->RGBColor[0.805646, 0.939986, 0.206806]], StyleBox[ "\n lambda=lambdak+eta; u=uk+dr+eta*v;\n \ r=FormInternalForceOfTwoBarArch[sprop,u]-lambda*force; r=N[r];\n", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, InitializationCell->True, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, Background->RGBColor[0.805646, 0.939986, 0.206806]], StyleBox[" Print[\"exit lambda=\",lambda,\" u=\",u,\" r=\",r];", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, InitializationCell->True, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, FontColor->RGBColor[1, 0, 0], Background->RGBColor[0.805646, 0.939986, 0.206806]], StyleBox["\n Return[{{lambda,u,r},\" \"}]];", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, InitializationCell->True, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, Background->RGBColor[0.805646, 0.939986, 0.206806]] }], "Input", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, InitializationCell->True, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, Background->RGBColor[0.805646, 0.939986, 0.206806]], Cell[TextData[{ StyleBox[ "Cell 13. The main program to run the two-bar arch structure. S=2, H=1, \ Em=1, A0=1,\nq={0,-1}, ell=0.010. Method: Forward Euler (FE), with", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, Evaluatable->False, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, Background->RGBColor[0, 1, 1]], StyleBox[" Load Control ", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, Evaluatable->False, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, FontColor->RGBColor[1, 0, 0], Background->RGBColor[0, 1, 1]], StyleBox["(LC) and Conventional Newton", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, Evaluatable->False, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, Background->RGBColor[0, 1, 1]] }], "Text", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, Evaluatable->False, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, Background->RGBColor[0, 1, 1]], Cell[CellGroupData[{ Cell[TextData[ "ClearAll[S,H,Em,A0]; \n\n(* Define solution method via method={integ,ics} \ *)\n\nmethod={\"FE\",\"LC\"};\n\n(* Define structural properties: \ {S,H,Em,A0} (see Cell 1) *)\n\nS=2; H=1; Em=1; A0=1; sprop={S,H,Em,A0};\n\n(* \ Define applied force reference: fx=0, fy=-lambda *)\n\nforce={0,-1};\n\n(* \ Define solution parameters *) \n\nnmax=20; lambdamax=2; umax={2*S,2*H}; \ ell=0.010; acctol=0; lref=H; \n\ solpar={nmax,lambdamax,umax,ell,ell,ell,acctol,H};\n\n(* Define reference \ state solution *)\n\nlambda0=0.; u0={0.,0.}; \n\ {v0,q0,status}=IncVelocity[sprop,force,{lambda0,u0}];\nIf[status!=\" \", \ Print[\"status=\",status]];\nkappa0=(q0.v0)/(v0.v0); \ Kdet=DetTanStiffness[sprop,u0];\nsol={0,lambda0,u0,v0,q0,Kdet,ell,1,1,kappa0,\ \"C0 cfg\"};\n\n(* Prepare Table to save solutions and store initial one *)\n\ \nsolnull={0,0,0,0,0,0,0,0,0,0,\" \"};\nsoltab=Table[solnull,{nmax+1}];\n\ soltab[[1]]=sol; \n\n(* Now do nmax steps and record solutions in soltab *)\ \n\nn=0;\nWhile [status==\" \" && nTrue,PlotRange->All,\n \ AxesOrigin->{0,0},AxesLabel->{\"uY\",\"lambda\"},\n \ PlotLabel->\"Response: lambda vs uY\"];\n\ ListPlot[uYvskappa,PlotJoined->True,PlotRange->All,\n \ AxesOrigin->{0,0},AxesLabel->{\"uY\",\"kappa\"},\n PlotLabel->\"LP \ sensor kappa vs uY\"];"], "Input", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, Background->RGBColor[0, 1, 0]], Cell[TextData[ " -7\nlambda=0.01 u={0., \ -0.0144539} r={0., 1.01594 10 } rnorm=1.015936701900116*10^-7\n \ -14\nlambda=0.01 u={0., -0.014454} r={0., \ 2.3561 10 } rnorm=2.356101425071699*10^-14\nNewton converged after 1 cycles\ \n -7\nlambda=0.02 u={0., \ -0.029584} r={0., 1.27413 10 } rnorm=1.274128222129855*10^-7\n \ -14\nlambda=0.02 u={0., -0.0295842} r={0., \ 4.01311 10 } rnorm=4.013109289324746*10^-14\nNewton converged after 1 \ cycles\n -7\nlambda=0.03 \ u={0., -0.0454821} r={0., 1.62766 10 } rnorm=1.627655716873733*10^-7\n \ -14\nlambda=0.03 u={0., \ -0.0454823} r={0., 7.14186 10 } rnorm=7.141856550596515*10^-14\nNewton \ converged after 1 cycles\n -7\ \nlambda=0.04 u={0., -0.0622625} r={0., 2.12721 10 } \ rnorm=2.127209427146837*10^-7\n \ -13\nlambda=0.04 u={0., -0.0622629} r={0., 1.34191 10 } \ rnorm=1.341912692076619*10^-13\nNewton converged after 1 cycles\n \ -7\nlambda=0.05 u={0., -0.0800705} \ r={0., 2.85757 10 } rnorm=2.857566261255107*10^-7\n \ -13\nlambda=0.05 u={0., -0.0800711} r={0., 2.69181 10 \ } rnorm=2.691805112142731*10^-13\nNewton converged after 1 cycles\n \ -7\nlambda=0.06 u={0., -0.0990956} \ r={0., 3.97067 10 } rnorm=3.970674723532963*10^-7\n \ -13\nlambda=0.06 u={0., -0.0990964} r={0., 5.85393 10 \ } rnorm=5.853928453092295*10^-13\nNewton converged after 1 cycles\n \ -7\nlambda=0.07 u={0., -0.119592} \ r={0., 5.75742 10 } rnorm=5.757420976637562*10^-7\n \ -12\nlambda=0.07 u={0., -0.119594} r={0., 1.40976 10 \ } rnorm=1.409761196669024*10^-12\nNewton converged after 1 cycles\n \ -7\nlambda=0.08 u={0., -0.141918} \ r={0., 8.82291 10 } rnorm=8.82291380427391*10^-7\n \ -12\nlambda=0.08 u={0., -0.14192} r={0., 3.87816 10 } \ rnorm=3.878161680681558*10^-12\nNewton converged after 1 cycles\n \ -6\nlambda=0.09 u={0., -0.166597} r={0., \ 1.457 10 } rnorm=1.45699628448237*10^-6\n \ -11\nlambda=0.09 u={0., -0.166601} r={0., 1.27831 10 } \ rnorm=1.278305239438282*10^-11\nNewton converged after 1 cycles\n \ -6\nlambda=0.1 u={0., -0.194466} r={0., \ 2.67688 10 } rnorm=2.676875172219484*10^-6\n \ -11\nlambda=0.1 u={0., -0.194474} r={0., 5.4654 10 } \ rnorm=5.465397578952037*10^-11\nNewton converged after 1 cycles\n \ -6\nlambda=0.11 u={0., -0.227001} r={0., \ 5.79878 10 } rnorm=5.798782686930461*10^-6\n \ -10\nlambda=0.11 u={0., -0.227021} r={0., 3.51095 10 } \ rnorm=3.510948892726518*10^-10\nNewton converged after 1 cycles\nlambda=0.12 \ u={0., -0.267322} r={0., 0.0000167657} rnorm=0.00001676568659284006\n \ -9\nlambda=0.12 u={0., -0.267399} \ r={0., 4.68925 10 } rnorm=4.689254123890763*10^-9\nNewton converged after 1 \ cycles\nlambda=0.13 u={0., -0.324949} r={0., 0.0000923426} \ rnorm=0.0000923425802657474\n \ -7\nlambda=0.13 u={0., -0.32566} r={0., 3.62349 10 } \ rnorm=3.623493324556026*10^-7\n \ -12\nlambda=0.13 u={0., -0.325663} r={0., 5.66389 10 } \ rnorm=5.663886026852083*10^-12\nNewton converged after 2 cycles\nlambda=0.14 \ u={0., -0.575685} r={0., 0.0169917} rnorm=0.0169917169583524\nlambda=0.14 \ u={0., -0.471178} r={0., 0.00531894} rnorm=0.005318943076695671\nlambda=0.14 \ u={0., -0.377759} r={0., 0.00518327} rnorm=0.005183272960622098\nlambda=0.14 \ u={0., -0.468506} r={0., 0.00517086} rnorm=0.00517086240471476\nlambda=0.14 \ u={0., -0.37263} r={0., 0.00549365} rnorm=0.005493648493994085\nlambda=0.14 \ u={0., -0.458581} r={0., 0.00469145} rnorm=0.004691445952969803\nlambda=0.14 \ u={0., -0.34855} r={0., 0.00742351} rnorm=0.007423507500187859\nlambda=0.14 \ u={0., -0.425416} r={0., 0.00392191} rnorm=0.003921913521463483\nlambda=0.14 \ u={0., 0.735152} r={0., 1.37354} rnorm=1.373535016237653\nn lambda \ u={uX,uY} v={vX,vY} Kdet a kappa remarks\n0 0 \ {0, 0} {0, -1.41421} 0.5 1 1 C0 cfg\n1 0.01 \ {0, -0.0144539} {0, -1.41421} 0.472212 1. 1. \n2 0.02 {0, \ -0.029584} {0, -1.47783} 0.443635 1. 0.956952 \n3 0.03 {0, \ -0.0454821} {0, -1.54972} 0.414783 1. 0.912561 \n4 0.04 {0, \ -0.0622625} {0, -1.6318} 0.385599 1. 0.866657 \n5 0.05 {0, \ -0.0800705} {0, -1.7267} 0.356014 1. 0.819027 \n6 0.06 {0, \ -0.0990956} {0, -1.83806} 0.325933 1. 0.769405 \n7 0.07 {0, \ -0.119592} {0, -1.97119} 0.295232 1. 0.717443 \n8 0.08 {0, \ -0.141918} {0, -2.13409} 0.263738 1. 0.662676 \n9 0.09 {0, \ -0.166597} {0, -2.33964} 0.231197 1. 0.604458 \n10 0.1 {0, \ -0.194466} {0, -2.61002} 0.19722 1. 0.54184 \n11 0.11 {0, \ -0.227001} {0, -2.98781} 0.16115 1. 0.473327 \n12 0.12 {0, \ -0.267322} {0, -3.56861} 0.121733 1. 0.396292 \nNewton didnt converge \ in max iters"], "Print", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, Evaluatable->False, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, FontSize->10, Background->RGBColor[1, 1, 0]], Cell[CellGroupData[{ Cell[GraphicsData["PostScript", "\<\ %! %%Creator: Mathematica %%AspectRatio: .61803 MathPictureStart %% Graphics /Courier findfont 10 scalefont setfont % Scaling calculations 0.97619 2.92446 0.0147151 4.52772 [ [(-0.3)] .09885 .01472 0 2 Msboxa [(-0.25)] .24508 .01472 0 2 Msboxa [(-0.2)] .3913 .01472 0 2 Msboxa [(-0.15)] .53752 .01472 0 2 Msboxa [(-0.1)] .68374 .01472 0 2 Msboxa [(-0.05)] .82997 .01472 0 2 Msboxa [(uY)] 1.025 .01472 -1 0 Msboxa [(Response: lambda vs uY)] .5 .61803 0 -2 0 0 1 Mouter Mrotsboxa [(0.02)] .96369 .10527 1 0 Msboxa [(0.04)] .96369 .19582 1 0 Msboxa [(0.06)] .96369 .28638 1 0 Msboxa [(0.08)] .96369 .37693 1 0 Msboxa [(0.1)] .96369 .46749 1 0 Msboxa [(0.12)] .96369 .55804 1 0 Msboxa [(lambda)] .97619 .61803 0 -4 Msboxa [ -0.001 -0.001 0 0 ] [ 1.001 .61903 0 0 ] ] MathScale % Start of Graphics 1 setlinecap 1 setlinejoin newpath [ ] 0 setdash 0 g p p .002 w .09885 .01472 m .09885 .02097 L s P [(-0.3)] .09885 .01472 0 2 Mshowa p .002 w .24508 .01472 m .24508 .02097 L s P [(-0.25)] .24508 .01472 0 2 Mshowa p .002 w .3913 .01472 m .3913 .02097 L s P [(-0.2)] .3913 .01472 0 2 Mshowa p .002 w .53752 .01472 m .53752 .02097 L s P [(-0.15)] .53752 .01472 0 2 Mshowa p .002 w .68374 .01472 m .68374 .02097 L s P [(-0.1)] .68374 .01472 0 2 Mshowa p .002 w .82997 .01472 m .82997 .02097 L s P [(-0.05)] .82997 .01472 0 2 Mshowa p .001 w .1281 .01472 m .1281 .01847 L s P p .001 w .15734 .01472 m .15734 .01847 L s P p .001 w .18659 .01472 m .18659 .01847 L s P p .001 w .21583 .01472 m .21583 .01847 L s P p .001 w .27432 .01472 m .27432 .01847 L s P p .001 w .30356 .01472 m .30356 .01847 L s P p .001 w .33281 .01472 m .33281 .01847 L s P p .001 w .36205 .01472 m .36205 .01847 L s P p .001 w .42054 .01472 m .42054 .01847 L s P p .001 w .44979 .01472 m .44979 .01847 L s P p .001 w .47903 .01472 m .47903 .01847 L s P p .001 w .50828 .01472 m .50828 .01847 L s P p .001 w .56677 .01472 m .56677 .01847 L s P p .001 w .59601 .01472 m .59601 .01847 L s P p .001 w .62526 .01472 m .62526 .01847 L s P p .001 w .6545 .01472 m .6545 .01847 L s P p .001 w .71299 .01472 m .71299 .01847 L s P p .001 w .74223 .01472 m .74223 .01847 L s P p .001 w .77148 .01472 m .77148 .01847 L s P p .001 w .80072 .01472 m .80072 .01847 L s P p .001 w .85921 .01472 m .85921 .01847 L s P p .001 w .88846 .01472 m .88846 .01847 L s P p .001 w .9177 .01472 m .9177 .01847 L s P p .001 w .94695 .01472 m .94695 .01847 L s P p .001 w .06961 .01472 m .06961 .01847 L s P p .001 w .04036 .01472 m .04036 .01847 L s P p .001 w .01112 .01472 m .01112 .01847 L s P [(uY)] 1.025 .01472 -1 0 Mshowa p .002 w 0 .01472 m 1 .01472 L s P [(Response: lambda vs uY)] .5 .61803 0 -2 0 0 1 Mouter Mrotshowa p .002 w .97619 .10527 m .98244 .10527 L s P [(0.02)] .96369 .10527 1 0 Mshowa p .002 w .97619 .19582 m .98244 .19582 L s P [(0.04)] .96369 .19582 1 0 Mshowa p .002 w .97619 .28638 m .98244 .28638 L s P [(0.06)] .96369 .28638 1 0 Mshowa p .002 w .97619 .37693 m .98244 .37693 L s P [(0.08)] .96369 .37693 1 0 Mshowa p .002 w .97619 .46749 m .98244 .46749 L s P [(0.1)] .96369 .46749 1 0 Mshowa p .002 w .97619 .55804 m .98244 .55804 L s P [(0.12)] .96369 .55804 1 0 Mshowa p .001 w .97619 .03283 m .97994 .03283 L s P p .001 w .97619 .05094 m .97994 .05094 L s P p .001 w .97619 .06905 m .97994 .06905 L s P p .001 w .97619 .08716 m .97994 .08716 L s P p .001 w .97619 .12338 m .97994 .12338 L s P p .001 w .97619 .14149 m .97994 .14149 L s P p .001 w .97619 .1596 m .97994 .1596 L s P p .001 w .97619 .17771 m .97994 .17771 L s P p .001 w .97619 .21393 m .97994 .21393 L s P p .001 w .97619 .23205 m .97994 .23205 L s P p .001 w .97619 .25016 m .97994 .25016 L s P p .001 w .97619 .26827 m .97994 .26827 L s P p .001 w .97619 .30449 m .97994 .30449 L s P p .001 w .97619 .3226 m .97994 .3226 L s P p .001 w .97619 .34071 m .97994 .34071 L s P p .001 w .97619 .35882 m .97994 .35882 L s P p .001 w .97619 .39504 m .97994 .39504 L s P p .001 w .97619 .41315 m .97994 .41315 L s P p .001 w .97619 .43127 m .97994 .43127 L s P p .001 w .97619 .44938 m .97994 .44938 L s P p .001 w .97619 .4856 m .97994 .4856 L s P p .001 w .97619 .50371 m .97994 .50371 L s P p .001 w .97619 .52182 m .97994 .52182 L s P p .001 w .97619 .53993 m .97994 .53993 L s P p .001 w .97619 .57615 m .97994 .57615 L s P p .001 w .97619 .59426 m .97994 .59426 L s P p .001 w .97619 .61237 m .97994 .61237 L s P [(lambda)] .97619 .61803 0 -4 Mshowa p .002 w .97619 0 m .97619 .61803 L s P P .004 w .97619 .01472 m .93392 .05999 L .88967 .10527 L .84318 .15055 L .79411 .19582 L .74203 .2411 L .68639 .28638 L .62645 .33166 L .56116 .37693 L .48898 .42221 L .40748 .46749 L .31234 .51276 L .19442 .55804 L .02381 .60332 L s 0 0 m 1 0 L 1 .61803 L 0 .61803 L closepath clip newpath % End of Graphics MathPictureEnd \ \>"], "Graphics", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, Evaluatable->False, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, ImageSize->{282, 174}, ImageMargins->{{34, Inherited}, {Inherited, Inherited}}, ImageCacheValid->False, Background->RGBColor[1, 1, 0]], Cell[GraphicsData["PostScript", "\<\ %! %%Creator: Mathematica %%AspectRatio: .61803 MathPictureStart %% Graphics /Courier findfont 10 scalefont setfont % Scaling calculations 0.97619 2.92446 -0.243869 0.847187 [ [(-0.3)] .09885 -0.24387 0 2 Msboxa [(-0.25)] .24508 -0.24387 0 2 Msboxa [(-0.2)] .3913 -0.24387 0 2 Msboxa [(-0.15)] .53752 -0.24387 0 2 Msboxa [(-0.1)] .68374 -0.24387 0 2 Msboxa [(-0.05)] .82997 -0.24387 0 2 Msboxa [(uY)] 1.025 -0.24387 -1 0 Msboxa [(LP sensor kappa vs uY)] .5 .61803 0 -2 0 0 1 Mouter Mrotsboxa [(0.3)] .96369 .01029 1 0 Msboxa [(0.4)] .96369 .09501 1 0 Msboxa [(0.5)] .96369 .17973 1 0 Msboxa [(0.6)] .96369 .26444 1 0 Msboxa [(0.7)] .96369 .34916 1 0 Msboxa [(0.8)] .96369 .43388 1 0 Msboxa [(0.9)] .96369 .5186 1 0 Msboxa [(1)] .96369 .60332 1 0 Msboxa [(kappa)] .97619 .61803 0 -4 Msboxa [ -0.001 -0.24487 0 0 ] [ 1.001 .61903 0 0 ] ] MathScale % Start of Graphics 1 setlinecap 1 setlinejoin newpath [ ] 0 setdash 0 g p p .002 w .09885 -0.24387 m .09885 -0.23762 L s P [(-0.3)] .09885 -0.24387 0 2 Mshowa p .002 w .24508 -0.24387 m .24508 -0.23762 L s P [(-0.25)] .24508 -0.24387 0 2 Mshowa p .002 w .3913 -0.24387 m .3913 -0.23762 L s P [(-0.2)] .3913 -0.24387 0 2 Mshowa p .002 w .53752 -0.24387 m .53752 -0.23762 L s P [(-0.15)] .53752 -0.24387 0 2 Mshowa p .002 w .68374 -0.24387 m .68374 -0.23762 L s P [(-0.1)] .68374 -0.24387 0 2 Mshowa p .002 w .82997 -0.24387 m .82997 -0.23762 L s P [(-0.05)] .82997 -0.24387 0 2 Mshowa p .001 w .1281 -0.24387 m .1281 -0.24012 L s P p .001 w .15734 -0.24387 m .15734 -0.24012 L s P p .001 w .18659 -0.24387 m .18659 -0.24012 L s P p .001 w .21583 -0.24387 m .21583 -0.24012 L s P p .001 w .27432 -0.24387 m .27432 -0.24012 L s P p .001 w .30356 -0.24387 m .30356 -0.24012 L s P p .001 w .33281 -0.24387 m .33281 -0.24012 L s P p .001 w .36205 -0.24387 m .36205 -0.24012 L s P p .001 w .42054 -0.24387 m .42054 -0.24012 L s P p .001 w .44979 -0.24387 m .44979 -0.24012 L s P p .001 w .47903 -0.24387 m .47903 -0.24012 L s P p .001 w .50828 -0.24387 m .50828 -0.24012 L s P p .001 w .56677 -0.24387 m .56677 -0.24012 L s P p .001 w .59601 -0.24387 m .59601 -0.24012 L s P p .001 w .62526 -0.24387 m .62526 -0.24012 L s P p .001 w .6545 -0.24387 m .6545 -0.24012 L s P p .001 w .71299 -0.24387 m .71299 -0.24012 L s P p .001 w .74223 -0.24387 m .74223 -0.24012 L s P p .001 w .77148 -0.24387 m .77148 -0.24012 L s P p .001 w .80072 -0.24387 m .80072 -0.24012 L s P p .001 w .85921 -0.24387 m .85921 -0.24012 L s P p .001 w .88846 -0.24387 m .88846 -0.24012 L s P p .001 w .9177 -0.24387 m .9177 -0.24012 L s P p .001 w .94695 -0.24387 m .94695 -0.24012 L s P p .001 w .06961 -0.24387 m .06961 -0.24012 L s P p .001 w .04036 -0.24387 m .04036 -0.24012 L s P p .001 w .01112 -0.24387 m .01112 -0.24012 L s P [(uY)] 1.025 -0.24387 -1 0 Mshowa p .002 w 0 -0.24387 m 1 -0.24387 L s P [(LP sensor kappa vs uY)] .5 .61803 0 -2 0 0 1 Mouter Mrotshowa p .002 w .97619 .01029 m .98244 .01029 L s P [(0.3)] .96369 .01029 1 0 Mshowa p .002 w .97619 .09501 m .98244 .09501 L s P [(0.4)] .96369 .09501 1 0 Mshowa p .002 w .97619 .17973 m .98244 .17973 L s P [(0.5)] .96369 .17973 1 0 Mshowa p .002 w .97619 .26444 m .98244 .26444 L s P [(0.6)] .96369 .26444 1 0 Mshowa p .002 w .97619 .34916 m .98244 .34916 L s P [(0.7)] .96369 .34916 1 0 Mshowa p .002 w .97619 .43388 m .98244 .43388 L s P [(0.8)] .96369 .43388 1 0 Mshowa p .002 w .97619 .5186 m .98244 .5186 L s P [(0.9)] .96369 .5186 1 0 Mshowa p .002 w .97619 .60332 m .98244 .60332 L s P [(1)] .96369 .60332 1 0 Mshowa p .001 w .97619 .02723 m .97994 .02723 L s P p .001 w .97619 .04418 m .97994 .04418 L s P p .001 w .97619 .06112 m .97994 .06112 L s P p .001 w .97619 .07806 m .97994 .07806 L s P p .001 w .97619 .11195 m .97994 .11195 L s P p .001 w .97619 .12889 m .97994 .12889 L s P p .001 w .97619 .14584 m .97994 .14584 L s P p .001 w .97619 .16278 m .97994 .16278 L s P p .001 w .97619 .19667 m .97994 .19667 L s P p .001 w .97619 .21361 m .97994 .21361 L s P p .001 w .97619 .23056 m .97994 .23056 L s P p .001 w .97619 .2475 m .97994 .2475 L s P p .001 w .97619 .28139 m .97994 .28139 L s P p .001 w .97619 .29833 m .97994 .29833 L s P p .001 w .97619 .31528 m .97994 .31528 L s P p .001 w .97619 .33222 m .97994 .33222 L s P p .001 w .97619 .36611 m .97994 .36611 L s P p .001 w .97619 .38305 m .97994 .38305 L s P p .001 w .97619 .39999 m .97994 .39999 L s P p .001 w .97619 .41694 m .97994 .41694 L s P p .001 w .97619 .45083 m .97994 .45083 L s P p .001 w .97619 .46777 m .97994 .46777 L s P p .001 w .97619 .48471 m .97994 .48471 L s P p .001 w .97619 .50166 m .97994 .50166 L s P p .001 w .97619 .53554 m .97994 .53554 L s P p .001 w .97619 .55249 m .97994 .55249 L s P p .001 w .97619 .56943 m .97994 .56943 L s P p .001 w .97619 .58638 m .97994 .58638 L s P [(kappa)] .97619 .61803 0 -4 Mshowa p .002 w .97619 0 m .97619 .61803 L s P P .004 w .97619 .60332 m .93392 .60332 L .88967 .56685 L .84318 .52924 L .79411 .49035 L .74203 .45 L .68639 .40796 L .62645 .36394 L .56116 .31754 L .48898 .26822 L .40748 .21517 L .31234 .15713 L .19442 .09187 L .02381 .01472 L s 0 0 m 1 0 L 1 .61803 L 0 .61803 L closepath clip newpath % End of Graphics MathPictureEnd \ \>"], "Graphics", CellFrame->True, CellMargins->{{14, 26}, {Inherited, Inherited}}, Evaluatable->False, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, ImageSize->{282, 174}, ImageMargins->{{34, Inherited}, {Inherited, Inherited}}, ImageCacheValid->False, Background->RGBColor[1, 1, 0]] }, Open ]] }, Open ]], Cell[TextData[{ StyleBox["\n", CellMargins->{{14, 26}, {Inherited, Inherited}}, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, Background->RGBColor[0, 1, 1]], StyleBox[ "Cell 14. The main program to run the two-bar arch structure. S=2, H=1, \ Em=1, A0=1,\nq={0,-1}, ell=0.100. Method: Forward Euler (FE), with", CellMargins->{{14, 26}, {Inherited, Inherited}}, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, FontFamily->"Times", FontWeight->"Plain", Background->RGBColor[0, 1, 1]], StyleBox[" DisplacementControl ", CellMargins->{{14, 26}, {Inherited, Inherited}}, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, FontFamily->"Times", FontWeight->"Plain", FontColor->RGBColor[1, 0, 0], Background->RGBColor[0, 1, 1]], StyleBox["(DC) and \nConventional Newton", CellMargins->{{14, 26}, {Inherited, Inherited}}, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, FontFamily->"Times", FontWeight->"Plain", Background->RGBColor[0, 1, 1]] }], "Input", CellMargins->{{14, 26}, {Inherited, Inherited}}, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, Background->RGBColor[0, 1, 1]], Cell[CellGroupData[{ Cell["\<\ ClearAll[S,H,Em,A0]; (* Define solution method via method={integ,ics} *) method={\"FE\",\"AC\"}; (* Define structural properties: {S,H,Em,A0} (see Cell 1) *) S=2; H=1; Em=1; A0=1; sprop={S,H,Em,A0}; (* Define applied force reference: fx=0, fy=-lambda *) force={0,-1}; (* Define solution parameters *) nmax=10; lambdamax=2; umax={2*S,2*H}; ell=0.100; acctol=0; lref=H; solpar={nmax,lambdamax,umax,ell,ell,ell,acctol,H}; (* Define reference state solution *) lambda0=0.; u0={0.,0.}; {v0,q0,status}=IncVelocity[sprop,force,{lambda0,u0}]; If[status!=\" \", Print[\"status=\",status]]; kappa0=(q0.v0)/(v0.v0); Kdet=DetTanStiffness[sprop,u0]; sol={0,lambda0,u0,v0,q0,Kdet,ell,1,1,kappa0,\"C0 cfg\"}; (* Prepare Table to save solutions and store initial one *) solnull={0,0,0,0,0,0,0,0,0,0,\" \"}; soltab=Table[solnull,{nmax+1}]; soltab[[1]]=sol; (* Now do nmax steps and record solutions in soltab *) n=0; While [status==\" \" && nTrue,PlotRange->All, AxesOrigin->{0,0},AxesLabel->{\"uY\",\"lambda\"}, PlotLabel->\"Response: lambda vs uY\"]; ListPlot[uYvskappa,PlotJoined->True,PlotRange->All, AxesOrigin->{0,0},AxesLabel->{\"uY\",\"kappa\"}, PlotLabel->\"LP sensor kappa vs uY\"];\ \>", "Input", CellMargins->{{14, 26}, {Inherited, Inherited}}, CellLabelMargins->{{7, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, Background->RGBColor[0, 1, 0]], Cell[BoxData[ \("Enter ACNewtoncycle"\)], "Print"], Cell[BoxData[ InterpretationBox[ \("sprop="\[InvisibleSpace]{2, 1, 1, 1} \[InvisibleSpace]" force ="\[InvisibleSpace]{0, \(-1\)}\), SequenceForm[ "sprop=", {2, 1, 1, 1}, " force =", {0, -1}], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("entry lambdak="\[InvisibleSpace]0.0577350269189625731` \[InvisibleSpace]" uk="\[InvisibleSpace]{0.`, \(-0.0816496580927726078`\)} \[InvisibleSpace]" rk="\[InvisibleSpace]{0.`, 0.00687861772213559152`}\), SequenceForm[ "entry lambdak=", 0.057735026918962574, " uk=", { 0.0, -0.081649658092772609}, " rk=", {0.0, 0.006878617722135591}], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("entry lambda0="\[InvisibleSpace]0.` \[InvisibleSpace]" u0="\[InvisibleSpace]{0.`, 0.`} \[InvisibleSpace]" f0="\[InvisibleSpace]1.73205080756887745`\), SequenceForm[ "entry lambda0=", 0.0, " u0=", {0.0, 0.0}, " f0=", 1.7320508075688774], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("v="\[InvisibleSpace]{0.`, \(-1.84852188262743624`\)} \[InvisibleSpace]" q="\[InvisibleSpace]{0, \(-1.`\)} \[InvisibleSpace]" status="\[InvisibleSpace]" "\), SequenceForm[ "v=", {0.0, -1.8485218826274363}, " q=", {0, -1.0}, " status=", " "], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("dr="\[InvisibleSpace]{0.`, \(-0.0127152753815965313`\)} \[InvisibleSpace]" r="\[InvisibleSpace]{0.`, 0.00687861772213559152`} \[InvisibleSpace]" status="\[InvisibleSpace]" "\), SequenceForm[ "dr=", {0.0, -0.01271527538159653}, " r=", {0.0, 0.006878617722135591}, " status=", " "], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("ck="\[InvisibleSpace]\(-1.38777878078144567`*^-17\)\), SequenceForm[ "ck=", -1.3877787807814457*^-17], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("eta="\[InvisibleSpace]\(-0.00497540020646834424`\)\), SequenceForm[ "eta=", -0.0049754002064683438], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("exit lambda="\[InvisibleSpace]0.0527596267124942386` \[InvisibleSpace]" u="\[InvisibleSpace]{0.`, \(-0.0851677973178833624`\)} \[InvisibleSpace]" r="\[InvisibleSpace]{0.`, 0.0000120408115251088854`}\), SequenceForm[ "exit lambda=", 0.052759626712494234, " u=", { 0.0, -0.085167797317883348}, " r=", {0.0, 1.2040811525108885*^-05}], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ RowBox[{ "\<\"lambda=\"\>", "\[InvisibleSpace]", "0.0527596267124942386`", "\[InvisibleSpace]", "\<\" u=\"\>", "\[InvisibleSpace]", \({0.`, \(-0.0851677973178833624`\)}\), "\[InvisibleSpace]", "\<\" r=\"\>", "\[InvisibleSpace]", \({0.`, 0.0000120408115251088854`}\), "\[InvisibleSpace]", "\<\" rnorm=\"\>", "\[InvisibleSpace]", InterpretationBox[ StyleBox["0.00001204081152510888", ShowStringCharacters->True, NumberMarks->True], InputForm[ 1.2040811525108887*^-05], Editable->True, AutoDelete->True]}], SequenceForm[ "lambda=", 0.052759626712494234, " u=", {0.0, -0.085167797317883348}, " r=", {0.0, 1.2040811525108885*^-05}, " rnorm=", InputForm[ 1.2040811525108887*^-05]], Editable->False]], "Print"], Cell[BoxData[ \("Enter ACNewtoncycle"\)], "Print"], Cell[BoxData[ InterpretationBox[ \("sprop="\[InvisibleSpace]{2, 1, 1, 1} \[InvisibleSpace]" force ="\[InvisibleSpace]{0, \(-1\)}\), SequenceForm[ "sprop=", {2, 1, 1, 1}, " force =", {0, -1}], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("entry lambdak="\[InvisibleSpace]0.0527596267124942386` \[InvisibleSpace]" uk="\[InvisibleSpace]{0.`, \(-0.0851677973178833624`\)} \[InvisibleSpace]" rk="\[InvisibleSpace]{0.`, 0.0000120408115251088854`}\), SequenceForm[ "entry lambdak=", 0.052759626712494234, " uk=", { 0.0, -0.085167797317883348}, " rk=", {0.0, 1.2040811525108885*^-05}], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("entry lambda0="\[InvisibleSpace]0.` \[InvisibleSpace]" u0="\[InvisibleSpace]{0.`, 0.`} \[InvisibleSpace]" f0="\[InvisibleSpace]1.73205080756887745`\), SequenceForm[ "entry lambda0=", 0.0, " u0=", {0.0, 0.0}, " f0=", 1.7320508075688774], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("v="\[InvisibleSpace]{0.`, \(-1.8721958403986477`\)} \[InvisibleSpace]" q="\[InvisibleSpace]{0, \(-1.`\)} \[InvisibleSpace]" status="\[InvisibleSpace]" "\), SequenceForm[ "v=", {0.0, -1.8721958403986476}, " q=", {0, -1.0}, " status=", " "], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("dr="\[InvisibleSpace]{0.`, \(-0.0000225427572523329544`\)} \[InvisibleSpace]" r="\[InvisibleSpace]{0.`, 0.0000120408115251088854`} \[InvisibleSpace]" status="\[InvisibleSpace]" "\), SequenceForm[ "dr=", {0.0, -2.2542757252332954*^-05}, " r=", {0.0, 1.2040811525108885*^-05}, " status=", " "], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[\("ck="\[InvisibleSpace]2.77555756156289135`*^-17\), SequenceForm[ "ck=", 2.7755575615628914*^-17], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("eta="\[InvisibleSpace]\(-8.73986521151696216`*^-6\)\), SequenceForm[ "eta=", -8.7398652115169629*^-06], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("exit lambda="\[InvisibleSpace]0.0527508868472827163` \[InvisibleSpace]" u="\[InvisibleSpace]{0.`, \(-0.0851739773358410268`\)} \[InvisibleSpace]" r="\[InvisibleSpace]{0.`, 3.70592168064121096`*^-11}\), SequenceForm[ "exit lambda=", 0.05275088684728272, " u=", { 0.0, -0.085173977335841033}, " r=", {0.0, 3.705921680641211*^-11}], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ RowBox[{ "\<\"lambda=\"\>", "\[InvisibleSpace]", "0.0527508868472827163`", "\[InvisibleSpace]", "\<\" u=\"\>", "\[InvisibleSpace]", \({0.`, \(-0.0851739773358410268`\)}\), "\[InvisibleSpace]", "\<\" r=\"\>", "\[InvisibleSpace]", \({0.`, 3.70592168064121096`*^-11}\), "\[InvisibleSpace]", "\<\" rnorm=\"\>", "\[InvisibleSpace]", InterpretationBox[ StyleBox["3.70592168064121*^-11", ShowStringCharacters->True, NumberMarks->True], InputForm[ 3.705921680641211*^-11], Editable->True, AutoDelete->True]}], SequenceForm[ "lambda=", 0.05275088684728272, " u=", {0.0, -0.085173977335841033}, " r=", {0.0, 3.705921680641211*^-11}, " rnorm=", InputForm[ 3.705921680641211*^-11]], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("Newton converged after "\[InvisibleSpace]1 \[InvisibleSpace]" cycles"\), SequenceForm[ "Newton converged after ", 1, " cycles"], Editable->False]], "Print"], Cell[BoxData[ \("Enter ACNewtoncycle"\)], "Print"], Cell[BoxData[ InterpretationBox[ \("sprop="\[InvisibleSpace]{2, 1, 1, 1} \[InvisibleSpace]" force ="\[InvisibleSpace]{0, \(-1\)}\), SequenceForm[ "sprop=", {2, 1, 1, 1}, " force =", {0, -1}], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("entry lambdak="\[InvisibleSpace]0.099873298277002398` \[InvisibleSpace]" uk="\[InvisibleSpace]{0.`, \(-0.173373817246863559`\)} \[InvisibleSpace]" rk="\[InvisibleSpace]{0.`, 0.00731883781761639973`}\), SequenceForm[ "entry lambdak=", 0.099873298277002392, " uk=", { 0.0, -0.17337381724686357}, " rk=", {0.0, 0.0073188378176163998}], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("entry lambda0="\[InvisibleSpace]0.0527596267124942386` \[InvisibleSpace]" u0="\[InvisibleSpace]{0.`, \(-0.0851677973178833624`\)} \[InvisibleSpace]" f0="\[InvisibleSpace]1.73205080756887745`\), SequenceForm[ "entry lambda0=", 0.052759626712494234, " u0=", { 0.0, -0.085167797317883348}, " f0=", 1.7320508075688774], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("v="\[InvisibleSpace]{0.`, \(-2.69391320134604006`\)} \[InvisibleSpace]" q="\[InvisibleSpace]{0, \(-1.`\)} \[InvisibleSpace]" status="\[InvisibleSpace]" "\), SequenceForm[ "v=", {0.0, -2.6939132013460401}, " q=", {0, -1.0}, " status=", " "], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("dr="\[InvisibleSpace]{0.`, \(-0.0197163138153874619`\)} \[InvisibleSpace]" r="\[InvisibleSpace]{0.`, 0.00731883781761639973`} \[InvisibleSpace]" status="\[InvisibleSpace]" "\), SequenceForm[ "dr=", {0.0, -0.019716313815387462}, " r=", {0.0, 0.0073188378176163998}, " status=", " "], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("ck="\[InvisibleSpace]\(-0.000778995350464398139`\)\), SequenceForm[ "ck=", -0.00077899535046439816], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("eta="\[InvisibleSpace]\(-0.00551665187250587418`\)\), SequenceForm[ "eta=", -0.0055166518725058749], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("exit lambda="\[InvisibleSpace]0.0943566464044965158` \[InvisibleSpace]" u="\[InvisibleSpace]{0.`, \(-0.178228749755677107`\)}\[InvisibleSpace]" r="\[InvisibleSpace]{ 0.`, 0.0000206253224779295862`}\), SequenceForm[ "exit lambda=", 0.094356646404496522, " u=", { 0.0, -0.1782287497556771}, " r=", {0.0, 2.0625322477929586*^-05}], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ RowBox[{ "\<\"lambda=\"\>", "\[InvisibleSpace]", "0.0943566464044965158`", "\[InvisibleSpace]", "\<\" u=\"\>", "\[InvisibleSpace]", \({0.`, \(-0.178228749755677107`\)}\), "\[InvisibleSpace]", "\<\" r=\"\>", "\[InvisibleSpace]", \({0.`, 0.0000206253224779295862`}\), "\[InvisibleSpace]", "\<\" rnorm=\"\>", "\[InvisibleSpace]", InterpretationBox[ StyleBox["0.00002062532247792958", ShowStringCharacters->True, NumberMarks->True], InputForm[ 2.0625322477929586*^-05], Editable->True, AutoDelete->True]}], SequenceForm[ "lambda=", 0.094356646404496522, " u=", {0.0, -0.1782287497556771}, " r=", {0.0, 2.0625322477929586*^-05}, " rnorm=", InputForm[ 2.0625322477929586*^-05]], Editable->False]], "Print"], Cell[BoxData[ \("Enter ACNewtoncycle"\)], "Print"], Cell[BoxData[ InterpretationBox[ \("sprop="\[InvisibleSpace]{2, 1, 1, 1} \[InvisibleSpace]" force ="\[InvisibleSpace]{0, \(-1\)}\), SequenceForm[ "sprop=", {2, 1, 1, 1}, " force =", {0, -1}], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("entry lambdak="\[InvisibleSpace]0.0943566464044965158` \[InvisibleSpace]" uk="\[InvisibleSpace]{0.`, \(-0.178228749755677107`\)} \[InvisibleSpace]" rk="\[InvisibleSpace]{0.`, 0.0000206253224779295862`}\), SequenceForm[ "entry lambdak=", 0.094356646404496522, " uk=", { 0.0, -0.1782287497556771}, " rk=", {0.0, 2.0625322477929586*^-05}], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("entry lambda0="\[InvisibleSpace]0.0527596267124942386` \[InvisibleSpace]" u0="\[InvisibleSpace]{0.`, \(-0.0851677973178833624`\)} \[InvisibleSpace]" f0="\[InvisibleSpace]1.73205080756887745`\), SequenceForm[ "entry lambda0=", 0.052759626712494234, " u0=", { 0.0, -0.085167797317883348}, " f0=", 1.7320508075688774], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("v="\[InvisibleSpace]{0.`, \(-2.75695590145596147`\)} \[InvisibleSpace]" q="\[InvisibleSpace]{0, \(-1.`\)} \[InvisibleSpace]" status="\[InvisibleSpace]" "\), SequenceForm[ "v=", {0.0, -2.7569559014559615}, " q=", {0, -1.0}, " status=", " "], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("dr="\[InvisibleSpace]{0.`, \(-0.0000568631045249602617`\)} \[InvisibleSpace]" r="\[InvisibleSpace]{0.`, 0.0000206253224779295862`} \[InvisibleSpace]" status="\[InvisibleSpace]" "\), SequenceForm[ "dr=", {0.0, -5.6863104524960269*^-05}, " r=", {0.0, 2.0625322477929586*^-05}, " status=", " "], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("ck="\[InvisibleSpace]\(-1.38777878078144567`*^-17\)\), SequenceForm[ "ck=", -1.3877787807814457*^-17], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("eta="\[InvisibleSpace]\(-0.0000164151488395025513`\)\), SequenceForm[ "eta=", -1.6415148839502549*^-05], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("exit lambda="\[InvisibleSpace]0.0943402312556570166` \[InvisibleSpace]" u="\[InvisibleSpace]{0.`, \(-0.178240357018735728`\)}\[InvisibleSpace]" r="\[InvisibleSpace]{ 0.`, 1.17431564472525451`*^-10}\), SequenceForm[ "exit lambda=", 0.09434023125565702, " u=", { 0.0, -0.17824035701873572}, " r=", {0.0, 1.1743156447252545*^-10}], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ RowBox[{ "\<\"lambda=\"\>", "\[InvisibleSpace]", "0.0943402312556570166`", "\[InvisibleSpace]", "\<\" u=\"\>", "\[InvisibleSpace]", \({0.`, \(-0.178240357018735728`\)}\), "\[InvisibleSpace]", "\<\" r=\"\>", "\[InvisibleSpace]", \({0.`, 1.17431564472525451`*^-10}\), "\[InvisibleSpace]", "\<\" rnorm=\"\>", "\[InvisibleSpace]", InterpretationBox[ StyleBox["1.174315644725254*^-10", ShowStringCharacters->True, NumberMarks->True], InputForm[ 1.1743156447252545*^-10], Editable->True, AutoDelete->True]}], SequenceForm[ "lambda=", 0.09434023125565702, " u=", {0.0, -0.17824035701873572}, " r=", {0.0, 1.1743156447252545*^-10}, " rnorm=", InputForm[ 1.1743156447252545*^-10]], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("Newton converged after "\[InvisibleSpace]1 \[InvisibleSpace]" cycles"\), SequenceForm[ "Newton converged after ", 1, " cycles"], Editable->False]], "Print"], Cell[BoxData[ \("Enter ACNewtoncycle"\)], "Print"], Cell[BoxData[ InterpretationBox[ \("sprop="\[InvisibleSpace]{2, 1, 1, 1} \[InvisibleSpace]" force ="\[InvisibleSpace]{0, \(-1\)}\), SequenceForm[ "sprop=", {2, 1, 1, 1}, " force =", {0, -1}], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("entry lambdak="\[InvisibleSpace]0.128454765873795228` \[InvisibleSpace]" uk="\[InvisibleSpace]{0.`, \(-0.272235761455110658`\)} \[InvisibleSpace]" rk="\[InvisibleSpace]{0.`, 0.00742968745976240008`}\), SequenceForm[ "entry lambdak=", 0.12845476587379523, " uk=", { 0.0, -0.27223576145511064}, " rk=", {0.0, 0.0074296874597623991}], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("entry lambda0="\[InvisibleSpace]0.0943566464044965158` \[InvisibleSpace]" u0="\[InvisibleSpace]{0.`, \(-0.178228749755677107`\)} \[InvisibleSpace]" f0="\[InvisibleSpace]1.73205080756887745`\), SequenceForm[ "entry lambda0=", 0.094356646404496522, " u0=", { 0.0, -0.1782287497556771}, " f0=", 1.7320508075688774], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("v="\[InvisibleSpace]{0.`, \(-4.80271647575888138`\)} \[InvisibleSpace]" q="\[InvisibleSpace]{0, \(-1.`\)} \[InvisibleSpace]" status="\[InvisibleSpace]" "\), SequenceForm[ "v=", {0.0, -4.8027164757588814}, " q=", {0, -1.0}, " status=", " "], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("dr="\[InvisibleSpace]{0.`, \(-0.0356826823727400288`\)} \[InvisibleSpace]" r="\[InvisibleSpace]{0.`, 0.00742968745976240008`} \[InvisibleSpace]" status="\[InvisibleSpace]" "\), SequenceForm[ "dr=", {0.0, -0.035682682372740027}, " r=", {0.0, 0.0074296874597623991}, " status=", " "], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[\("ck="\[InvisibleSpace]\(-0.00355703790972014566`\)\), SequenceForm[ "ck=", -0.0035570379097201454], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("eta="\[InvisibleSpace]\(-0.00568552147221354786`\)\), SequenceForm[ "eta=", -0.0056855214722135484], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("exit lambda="\[InvisibleSpace]0.122769244401581678` \[InvisibleSpace]" u="\[InvisibleSpace]{0.`, \(-0.280612496179969728`\)}\[InvisibleSpace]" r="\[InvisibleSpace]{ 0.`, 0.0000539569031801911958`}\), SequenceForm[ "exit lambda=", 0.12276924440158168, " u=", { 0.0, -0.28061249617996975}, " r=", {0.0, 5.3956903180191196*^-05}], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ RowBox[{ "\<\"lambda=\"\>", "\[InvisibleSpace]", "0.122769244401581678`", "\[InvisibleSpace]", "\<\" u=\"\>", "\[InvisibleSpace]", \({0.`, \(-0.280612496179969728`\)}\), "\[InvisibleSpace]", "\<\" r=\"\>", "\[InvisibleSpace]", \({0.`, 0.0000539569031801911958`}\), "\[InvisibleSpace]", "\<\" rnorm=\"\>", "\[InvisibleSpace]", InterpretationBox[ StyleBox["0.00005395690318019119", ShowStringCharacters->True, NumberMarks->True], InputForm[ 5.3956903180191196*^-05], Editable->True, AutoDelete->True]}], SequenceForm[ "lambda=", 0.12276924440158168, " u=", {0.0, -0.28061249617996975}, " r=", {0.0, 5.3956903180191196*^-05}, " rnorm=", InputForm[ 5.3956903180191196*^-05]], Editable->False]], "Print"], Cell[BoxData[ \("Enter ACNewtoncycle"\)], "Print"], Cell[BoxData[ InterpretationBox[ \("sprop="\[InvisibleSpace]{2, 1, 1, 1} \[InvisibleSpace]" force ="\[InvisibleSpace]{0, \(-1\)}\), SequenceForm[ "sprop=", {2, 1, 1, 1}, " force =", {0, -1}], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("entry lambdak="\[InvisibleSpace]0.122769244401581678` \[InvisibleSpace]" uk="\[InvisibleSpace]{0.`, \(-0.280612496179969728`\)} \[InvisibleSpace]" rk="\[InvisibleSpace]{0.`, 0.0000539569031801911958`}\), SequenceForm[ "entry lambdak=", 0.12276924440158168, " uk=", { 0.0, -0.28061249617996975}, " rk=", {0.0, 5.3956903180191196*^-05}], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("entry lambda0="\[InvisibleSpace]0.0943566464044965158` \[InvisibleSpace]" u0="\[InvisibleSpace]{0.`, \(-0.178228749755677107`\)} \[InvisibleSpace]" f0="\[InvisibleSpace]1.73205080756887745`\), SequenceForm[ "entry lambda0=", 0.094356646404496522, " u0=", { 0.0, -0.1782287497556771}, " f0=", 1.7320508075688774], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("v="\[InvisibleSpace]{0.`, \(-5.11881423223649978`\)} \[InvisibleSpace]" q="\[InvisibleSpace]{0, \(-1.`\)} \[InvisibleSpace]" status="\[InvisibleSpace]" "\), SequenceForm[ "v=", {0.0, -5.1188142322364998}, " q=", {0, -1.0}, " status=", " "], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("dr="\[InvisibleSpace]{0.`, \(-0.000276195363926169523`\)} \[InvisibleSpace]" r="\[InvisibleSpace]{0.`, 0.0000539569031801911958`} \[InvisibleSpace]" status="\[InvisibleSpace]" "\), SequenceForm[ "dr=", {0.0, -0.00027619536392616955}, " r=", {0.0, 5.3956903180191196*^-05}, " status=", " "], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("ck="\[InvisibleSpace]\(-1.38777878078144567`*^-17\)\), SequenceForm[ "ck=", -1.3877787807814457*^-17], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("eta="\[InvisibleSpace]\(-0.000047408017254231991`\)\), SequenceForm[ "eta=", -4.740801725423199*^-05], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("exit lambda="\[InvisibleSpace]0.122721836384327454` \[InvisibleSpace]" u="\[InvisibleSpace]{0.`, \(-0.280646018710452782`\)}\[InvisibleSpace]" r="\[InvisibleSpace]{ 0.`, 8.57444434898546603`*^-10}\), SequenceForm[ "exit lambda=", 0.12272183638432745, " u=", { 0.0, -0.2806460187104528}, " r=", {0.0, 8.574444348985466*^-10}], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ RowBox[{ "\<\"lambda=\"\>", "\[InvisibleSpace]", "0.122721836384327454`", "\[InvisibleSpace]", "\<\" u=\"\>", "\[InvisibleSpace]", \({0.`, \(-0.280646018710452782`\)}\), "\[InvisibleSpace]", "\<\" r=\"\>", "\[InvisibleSpace]", \({0.`, 8.57444434898546603`*^-10}\), "\[InvisibleSpace]", "\<\" rnorm=\"\>", "\[InvisibleSpace]", InterpretationBox[ StyleBox["8.574444348985466*^-10", ShowStringCharacters->True, NumberMarks->True], InputForm[ 8.574444348985466*^-10], Editable->True, AutoDelete->True]}], SequenceForm[ "lambda=", 0.12272183638432745, " u=", {0.0, -0.2806460187104528}, " r=", {0.0, 8.574444348985466*^-10}, " rnorm=", InputForm[ 8.574444348985466*^-10]], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("Newton converged after "\[InvisibleSpace]1 \[InvisibleSpace]" cycles"\), SequenceForm[ "Newton converged after ", 1, " cycles"], Editable->False]], "Print"], Cell[BoxData[ \("Enter ACNewtoncycle"\)], "Print"], Cell[BoxData[ InterpretationBox[ \("sprop="\[InvisibleSpace]{2, 1, 1, 1} \[InvisibleSpace]" force ="\[InvisibleSpace]{0, \(-1\)}\), SequenceForm[ "sprop=", {2, 1, 1, 1}, " force =", {0, -1}], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("entry lambdak="\[InvisibleSpace]0.141942572901929731` \[InvisibleSpace]" uk="\[InvisibleSpace]{0.`, \(-0.378757202986897212`\)} \[InvisibleSpace]" rk="\[InvisibleSpace]{0.`, 0.00706947492184434267`}\), SequenceForm[ "entry lambdak=", 0.14194257290192974, " uk=", { 0.0, -0.37875720298689719}, " rk=", {0.0, 0.0070694749218443431}], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("entry lambda0="\[InvisibleSpace]0.122769244401581678` \[InvisibleSpace]" u0="\[InvisibleSpace]{0.`, \(-0.280612496179969728`\)} \[InvisibleSpace]" f0="\[InvisibleSpace]1.73205080756887745`\), SequenceForm[ "entry lambda0=", 0.12276924440158168, " u0=", { 0.0, -0.28061249617996975}, " f0=", 1.7320508075688774], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("v="\[InvisibleSpace]{0.`, \(-17.9209647121417559`\)} \[InvisibleSpace]" q="\[InvisibleSpace]{0, \(-1.`\)} \[InvisibleSpace]" status="\[InvisibleSpace]" "\), SequenceForm[ "v=", {0.0, -17.920964712141757}, " q=", {0, -1.0}, " status=", " "], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("dr="\[InvisibleSpace]{0.`, \(-0.126691810607743571`\)} \[InvisibleSpace]" r="\[InvisibleSpace]{0.`, 0.00706947492184434267`} \[InvisibleSpace]" status="\[InvisibleSpace]" "\), SequenceForm[ "dr=", {0.0, -0.12669181060774357}, " r=", {0.0, 0.0070694749218443431}, " status=", " "], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[\("ck="\[InvisibleSpace]\(-0.00879545608505255493`\)\), SequenceForm[ "ck=", -0.0087954560850525554], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("eta="\[InvisibleSpace]\(-0.00622284603736583541`\)\), SequenceForm[ "eta=", -0.0062228460373658359], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("exit lambda="\[InvisibleSpace]0.135719726864563905` \[InvisibleSpace]" u="\[InvisibleSpace]{0.`, \(-0.393929609349916489`\)}\[InvisibleSpace]" r="\[InvisibleSpace]{ 0.`, 0.000150451508418525392`}\), SequenceForm[ "exit lambda=", 0.13571972686456391, " u=", { 0.0, -0.3939296093499165}, " r=", {0.0, 0.00015045150841852539}], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ RowBox[{ "\<\"lambda=\"\>", "\[InvisibleSpace]", "0.135719726864563905`", "\[InvisibleSpace]", "\<\" u=\"\>", "\[InvisibleSpace]", \({0.`, \(-0.393929609349916489`\)}\), "\[InvisibleSpace]", "\<\" r=\"\>", "\[InvisibleSpace]", \({0.`, 0.000150451508418525392`}\), "\[InvisibleSpace]", "\<\" rnorm=\"\>", "\[InvisibleSpace]", InterpretationBox[ StyleBox["0.0001504515084185253", ShowStringCharacters->True, NumberMarks->True], InputForm[ 0.00015045150841852539], Editable->True, AutoDelete->True]}], SequenceForm[ "lambda=", 0.13571972686456391, " u=", {0.0, -0.3939296093499165}, " r=", {0.0, 0.00015045150841852539}, " rnorm=", InputForm[ 0.00015045150841852539]], Editable->False]], "Print"], Cell[BoxData[ \("Enter ACNewtoncycle"\)], "Print"], Cell[BoxData[ InterpretationBox[ \("sprop="\[InvisibleSpace]{2, 1, 1, 1} \[InvisibleSpace]" force ="\[InvisibleSpace]{0, \(-1\)}\), SequenceForm[ "sprop=", {2, 1, 1, 1}, " force =", {0, -1}], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("entry lambdak="\[InvisibleSpace]0.135719726864563905` \[InvisibleSpace]" uk="\[InvisibleSpace]{0.`, \(-0.393929609349916489`\)} \[InvisibleSpace]" rk="\[InvisibleSpace]{0.`, 0.000150451508418525392`}\), SequenceForm[ "entry lambdak=", 0.13571972686456391, " uk=", { 0.0, -0.3939296093499165}, " rk=", {0.0, 0.00015045150841852539}], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("entry lambda0="\[InvisibleSpace]0.122769244401581678` \[InvisibleSpace]" u0="\[InvisibleSpace]{0.`, \(-0.280612496179969728`\)} \[InvisibleSpace]" f0="\[InvisibleSpace]1.73205080756887745`\), SequenceForm[ "entry lambda0=", 0.12276924440158168, " u0=", { 0.0, -0.28061249617996975}, " f0=", 1.7320508075688774], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("v="\[InvisibleSpace]{0.`, \(-27.7394802634470849`\)} \[InvisibleSpace]" q="\[InvisibleSpace]{0, \(-1.`\)} \[InvisibleSpace]" status="\[InvisibleSpace]" "\), SequenceForm[ "v=", {0.0, -27.739480263447085}, " q=", {0, -1.0}, " status=", " "], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("dr="\[InvisibleSpace]{0.`, \(-0.00417344664838152823`\)} \[InvisibleSpace]" r="\[InvisibleSpace]{0.`, 0.000150451508418525392`} \[InvisibleSpace]" status="\[InvisibleSpace]" "\), SequenceForm[ "dr=", {0.0, -0.0041734466483815281}, " r=", {0.0, 0.00015045150841852539}, " status=", " "], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[\("ck="\[InvisibleSpace]2.77555756156289135`*^-17\), SequenceForm[ "ck=", 2.7755575615628914*^-17], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("eta="\[InvisibleSpace]\(-0.000146711682559081202`\)\), SequenceForm[ "eta=", -0.00014671168255908122], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("exit lambda="\[InvisibleSpace]0.135573015182004819` \[InvisibleSpace]" u="\[InvisibleSpace]{0.`, \(-0.394033350175533314`\)}\[InvisibleSpace]" r="\[InvisibleSpace]{ 0.`, 6.91789470153025831`*^-9}\), SequenceForm[ "exit lambda=", 0.13557301518200482, " u=", { 0.0, -0.39403335017553331}, " r=", {0.0, 6.9178947015302583*^-09}], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ RowBox[{ "\<\"lambda=\"\>", "\[InvisibleSpace]", "0.135573015182004819`", "\[InvisibleSpace]", "\<\" u=\"\>", "\[InvisibleSpace]", \({0.`, \(-0.394033350175533314`\)}\), "\[InvisibleSpace]", "\<\" r=\"\>", "\[InvisibleSpace]", \({0.`, 6.91789470153025831`*^-9}\), "\[InvisibleSpace]", "\<\" rnorm=\"\>", "\[InvisibleSpace]", InterpretationBox[ StyleBox["6.917894701530258*^-9", ShowStringCharacters->True, NumberMarks->True], InputForm[ 6.9178947015302583*^-09], Editable->True, AutoDelete->True]}], SequenceForm[ "lambda=", 0.13557301518200482, " u=", {0.0, -0.39403335017553331}, " r=", {0.0, 6.9178947015302583*^-09}, " rnorm=", InputForm[ 6.9178947015302583*^-09]], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("Newton converged after "\[InvisibleSpace]1 \[InvisibleSpace]" cycles"\), SequenceForm[ "Newton converged after ", 1, " cycles"], Editable->False]], "Print"], Cell[BoxData[ \("Enter ACNewtoncycle"\)], "Print"], Cell[BoxData[ InterpretationBox[ \("sprop="\[InvisibleSpace]{2, 1, 1, 1} \[InvisibleSpace]" force ="\[InvisibleSpace]{0, \(-1\)}\), SequenceForm[ "sprop=", {2, 1, 1, 1}, " force =", {0, -1}], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("entry lambdak="\[InvisibleSpace]0.13932235688017538` \[InvisibleSpace]" uk="\[InvisibleSpace]{0.`, \(-0.493864693564473622`\)} \[InvisibleSpace]" rk="\[InvisibleSpace]{0.`, 0.00621759013814204841`}\), SequenceForm[ "entry lambdak=", 0.13932235688017539, " uk=", { 0.0, -0.4938646935644736}, " rk=", {0.0, 0.0062175901381420484}], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("entry lambda0="\[InvisibleSpace]0.135719726864563905` \[InvisibleSpace]" u0="\[InvisibleSpace]{0.`, \(-0.393929609349916489`\)} \[InvisibleSpace]" f0="\[InvisibleSpace]1.73205080756887745`\), SequenceForm[ "entry lambda0=", 0.13571972686456391, " u0=", { 0.0, -0.3939296093499165}, " f0=", 1.7320508075688774], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("v="\[InvisibleSpace]{0.`, 12.2188224261473998`} \[InvisibleSpace]" q="\[InvisibleSpace]{0, \(-1.`\)} \[InvisibleSpace]" status="\[InvisibleSpace]" "\), SequenceForm[ "v=", {0.0, 12.218822426147399}, " q=", {0, -1.0}, " status=", " "], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("dr="\[InvisibleSpace]{0.`, 0.075971629816522963`} \[InvisibleSpace]" r="\[InvisibleSpace]{0.`, 0.00621759013814204841`} \[InvisibleSpace]" status="\[InvisibleSpace]" "\), SequenceForm[ "dr=", {0.0, 0.075971629816522973}, " r=", {0.0, 0.0062175901381420484}, " status=", " "], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[\("ck="\[InvisibleSpace]\(-0.0163233660147858579`\)\), SequenceForm[ "ck=", -0.016323366014785859], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("eta="\[InvisibleSpace]\(-0.00833616739560335595`\)\), SequenceForm[ "eta=", -0.0083361673956033559], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("exit lambda="\[InvisibleSpace]0.130986189484572035` \[InvisibleSpace]" u="\[InvisibleSpace]{0.`, \(-0.519751212869467682`\)}\[InvisibleSpace]" r="\[InvisibleSpace]{ 0.`, 0.000353608185588327739`}\), SequenceForm[ "exit lambda=", 0.13098618948457202, " u=", { 0.0, -0.51975121286946768}, " r=", {0.0, 0.00035360818558832774}], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ RowBox[{ "\<\"lambda=\"\>", "\[InvisibleSpace]", "0.130986189484572035`", "\[InvisibleSpace]", "\<\" u=\"\>", "\[InvisibleSpace]", \({0.`, \(-0.519751212869467682`\)}\), "\[InvisibleSpace]", "\<\" r=\"\>", "\[InvisibleSpace]", \({0.`, 0.000353608185588327739`}\), "\[InvisibleSpace]", "\<\" rnorm=\"\>", "\[InvisibleSpace]", InterpretationBox[ StyleBox["0.0003536081855883277", ShowStringCharacters->True, NumberMarks->True], InputForm[ 0.00035360818558832774], Editable->True, AutoDelete->True]}], SequenceForm[ "lambda=", 0.13098618948457202, " u=", {0.0, -0.51975121286946768}, " r=", {0.0, 0.00035360818558832774}, " rnorm=", InputForm[ 0.00035360818558832774]], Editable->False]], "Print"], Cell[BoxData[ \("Enter ACNewtoncycle"\)], "Print"], Cell[BoxData[ InterpretationBox[ \("sprop="\[InvisibleSpace]{2, 1, 1, 1} \[InvisibleSpace]" force ="\[InvisibleSpace]{0, \(-1\)}\), SequenceForm[ "sprop=", {2, 1, 1, 1}, " force =", {0, -1}], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("entry lambdak="\[InvisibleSpace]0.130986189484572035` \[InvisibleSpace]" uk="\[InvisibleSpace]{0.`, \(-0.519751212869467682`\)} \[InvisibleSpace]" rk="\[InvisibleSpace]{0.`, 0.000353608185588327739`}\), SequenceForm[ "entry lambdak=", 0.13098618948457202, " uk=", { 0.0, -0.51975121286946768}, " rk=", {0.0, 0.00035360818558832774}], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("entry lambda0="\[InvisibleSpace]0.135719726864563905` \[InvisibleSpace]" u0="\[InvisibleSpace]{0.`, \(-0.393929609349916489`\)} \[InvisibleSpace]" f0="\[InvisibleSpace]1.73205080756887745`\), SequenceForm[ "entry lambda0=", 0.13571972686456391, " u0=", { 0.0, -0.3939296093499165}, " f0=", 1.7320508075688774], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("v="\[InvisibleSpace]{0.`, 9.18072176259482475`} \[InvisibleSpace]" q="\[InvisibleSpace]{0, \(-1.`\)} \[InvisibleSpace]" status="\[InvisibleSpace]" "\), SequenceForm[ "v=", {0.0, 9.1807217625948248}, " q=", {0, -1.0}, " status=", " "], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("dr="\[InvisibleSpace]{0.`, 0.00324637836486243003`} \[InvisibleSpace]" r="\[InvisibleSpace]{0.`, 0.000353608185588327739`} \[InvisibleSpace]" status="\[InvisibleSpace]" "\), SequenceForm[ "dr=", {0.0, 0.00324637836486243}, " r=", {0.0, 0.00035360818558832774}, " status=", " "], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[\("ck="\[InvisibleSpace]0.`\), SequenceForm[ "ck=", 0.0], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("eta="\[InvisibleSpace]\(-0.000383116104754018138`\)\), SequenceForm[ "eta=", -0.00038311610475401809], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("exit lambda="\[InvisibleSpace]0.130603073379818002` \[InvisibleSpace]" u="\[InvisibleSpace]{0.`, \(-0.520022116865121031`\)}\[InvisibleSpace]" r="\[InvisibleSpace]{ 0.`, 3.73759027338937954`*^-8}\), SequenceForm[ "exit lambda=", 0.130603073379818, " u=", {0.0, -0.52002211686512101}, " r=", {0.0, 3.7375902733893795*^-08}], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ RowBox[{ "\<\"lambda=\"\>", "\[InvisibleSpace]", "0.130603073379818002`", "\[InvisibleSpace]", "\<\" u=\"\>", "\[InvisibleSpace]", \({0.`, \(-0.520022116865121031`\)}\), "\[InvisibleSpace]", "\<\" r=\"\>", "\[InvisibleSpace]", \({0.`, 3.73759027338937954`*^-8}\), "\[InvisibleSpace]", "\<\" rnorm=\"\>", "\[InvisibleSpace]", InterpretationBox[ StyleBox["3.737590273389379*^-8", ShowStringCharacters->True, NumberMarks->True], InputForm[ 3.7375902733893795*^-08], Editable->True, AutoDelete->True]}], SequenceForm[ "lambda=", 0.130603073379818, " u=", {0.0, -0.52002211686512101}, " r=", {0.0, 3.7375902733893795*^-08}, " rnorm=", InputForm[ 3.7375902733893795*^-08]], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("Newton converged after "\[InvisibleSpace]1 \[InvisibleSpace]" cycles"\), SequenceForm[ "Newton converged after ", 1, " cycles"], Editable->False]], "Print"], Cell[BoxData[ \("Enter ACNewtoncycle"\)], "Print"], Cell[BoxData[ InterpretationBox[ \("sprop="\[InvisibleSpace]{2, 1, 1, 1} \[InvisibleSpace]" force ="\[InvisibleSpace]{0, \(-1\)}\), SequenceForm[ "sprop=", {2, 1, 1, 1}, " force =", {0, -1}], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("entry lambdak="\[InvisibleSpace]0.120157846266779388` \[InvisibleSpace]" uk="\[InvisibleSpace]{0.`, \(-0.619163219101902662`\)} \[InvisibleSpace]" rk="\[InvisibleSpace]{0.`, 0.00504033599125684794`}\), SequenceForm[ "entry lambdak=", 0.12015784626677939, " uk=", { 0.0, -0.61916321910190264}, " rk=", {0.0, 0.0050403359912568474}], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("entry lambda0="\[InvisibleSpace]0.130986189484572035` \[InvisibleSpace]" u0="\[InvisibleSpace]{0.`, \(-0.519751212869467682`\)} \[InvisibleSpace]" f0="\[InvisibleSpace]1.73205080756887745`\), SequenceForm[ "entry lambda0=", 0.13098618948457202, " u0=", { 0.0, -0.51975121286946768}, " f0=", 1.7320508075688774], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("v="\[InvisibleSpace]{0.`, 5.00704018436224718`} \[InvisibleSpace]" q="\[InvisibleSpace]{0, \(-1.`\)} \[InvisibleSpace]" status="\[InvisibleSpace]" "\), SequenceForm[ "v=", {0.0, 5.0070401843622472}, " q=", {0, -1.0}, " status=", " "], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("dr="\[InvisibleSpace]{0.`, 0.0252371648509103607`} \[InvisibleSpace]" r="\[InvisibleSpace]{0.`, 0.00504033599125684794`} \[InvisibleSpace]" status="\[InvisibleSpace]" "\), SequenceForm[ "dr=", {0.0, 0.025237164850910359}, " r=", {0.0, 0.0050403359912568474}, " status=", " "], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[\("ck="\[InvisibleSpace]\(-0.0250821836797213126`\)\), SequenceForm[ "ck=", -0.025082183679721315], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[\("eta="\[InvisibleSpace]\(-0.0130133272842375746`\)\), SequenceForm[ "eta=", -0.013013327284237574], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("exit lambda="\[InvisibleSpace]0.107144518982541825` \[InvisibleSpace]" u="\[InvisibleSpace]{0.`, \(-0.659084306895427474`\)}\[InvisibleSpace]" r="\[InvisibleSpace]{ 0.`, 0.000621260147863053635`}\), SequenceForm[ "exit lambda=", 0.10714451898254182, " u=", { 0.0, -0.65908430689542752}, " r=", {0.0, 0.00062126014786305361}], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ RowBox[{ "\<\"lambda=\"\>", "\[InvisibleSpace]", "0.107144518982541825`", "\[InvisibleSpace]", "\<\" u=\"\>", "\[InvisibleSpace]", \({0.`, \(-0.659084306895427474`\)}\), "\[InvisibleSpace]", "\<\" r=\"\>", "\[InvisibleSpace]", \({0.`, 0.000621260147863053635`}\), "\[InvisibleSpace]", "\<\" rnorm=\"\>", "\[InvisibleSpace]", InterpretationBox[ StyleBox["0.0006212601478630536", ShowStringCharacters->True, NumberMarks->True], InputForm[ 0.00062126014786305361], Editable->True, AutoDelete->True]}], SequenceForm[ "lambda=", 0.10714451898254182, " u=", {0.0, -0.65908430689542752}, " r=", {0.0, 0.00062126014786305361}, " rnorm=", InputForm[ 0.00062126014786305361]], Editable->False]], "Print"], Cell[BoxData[ \("Enter ACNewtoncycle"\)], "Print"], Cell[BoxData[ InterpretationBox[ \("sprop="\[InvisibleSpace]{2, 1, 1, 1} \[InvisibleSpace]" force ="\[InvisibleSpace]{0, \(-1\)}\), SequenceForm[ "sprop=", {2, 1, 1, 1}, " force =", {0, -1}], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("entry lambdak="\[InvisibleSpace]0.107144518982541825` \[InvisibleSpace]" uk="\[InvisibleSpace]{0.`, \(-0.659084306895427474`\)} \[InvisibleSpace]" rk="\[InvisibleSpace]{0.`, 0.000621260147863053635`}\), SequenceForm[ "entry lambdak=", 0.10714451898254182, " uk=", { 0.0, -0.65908430689542752}, " rk=", {0.0, 0.00062126014786305361}], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("entry lambda0="\[InvisibleSpace]0.130986189484572035` \[InvisibleSpace]" u0="\[InvisibleSpace]{0.`, \(-0.519751212869467682`\)} \[InvisibleSpace]" f0="\[InvisibleSpace]1.73205080756887745`\), SequenceForm[ "entry lambda0=", 0.13098618948457202, " u0=", { 0.0, -0.51975121286946768}, " f0=", 1.7320508075688774], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("v="\[InvisibleSpace]{0.`, 4.34254436883598149`} \[InvisibleSpace]" q="\[InvisibleSpace]{0, \(-1.`\)} \[InvisibleSpace]" status="\[InvisibleSpace]" "\), SequenceForm[ "v=", {0.0, 4.3425443688359815}, " q=", {0, -1.0}, " status=", " "], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("dr="\[InvisibleSpace]{0.`, 0.00269784975668491266`} \[InvisibleSpace]" r="\[InvisibleSpace]{0.`, 0.000621260147863053635`} \[InvisibleSpace]" status="\[InvisibleSpace]" "\), SequenceForm[ "dr=", {0.0, 0.0026978497566849126}, " r=", {0.0, 0.00062126014786305361}, " status=", " "], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[\("ck="\[InvisibleSpace]6.93889390390722837`*^-17\), SequenceForm[ "ck=", 6.9388939039072284*^-17], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("eta="\[InvisibleSpace]\(-0.000742097668200976556`\)\), SequenceForm[ "eta=", -0.00074209766820097648], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("exit lambda="\[InvisibleSpace]0.106402421314340833` \[InvisibleSpace]" u="\[InvisibleSpace]{0.`, \(-0.659609049188915008`\)}\[InvisibleSpace]" r="\[InvisibleSpace]{ 0.`, 9.9515908333347447`*^-8}\), SequenceForm[ "exit lambda=", 0.10640242131434084, " u=", { 0.0, -0.65960904918891505}, " r=", {0.0, 9.9515908333347447*^-08}], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ RowBox[{ "\<\"lambda=\"\>", "\[InvisibleSpace]", "0.106402421314340833`", "\[InvisibleSpace]", "\<\" u=\"\>", "\[InvisibleSpace]", \({0.`, \(-0.659609049188915008`\)}\), "\[InvisibleSpace]", "\<\" r=\"\>", "\[InvisibleSpace]", \({0.`, 9.9515908333347447`*^-8}\), "\[InvisibleSpace]", "\<\" rnorm=\"\>", "\[InvisibleSpace]", InterpretationBox[ StyleBox["9.951590833334744*^-8", ShowStringCharacters->True, NumberMarks->True], InputForm[ 9.9515908333347447*^-08], Editable->True, AutoDelete->True]}], SequenceForm[ "lambda=", 0.10640242131434084, " u=", {0.0, -0.65960904918891505}, " r=", {0.0, 9.9515908333347447*^-08}, " rnorm=", InputForm[ 9.9515908333347447*^-08]], Editable->False]], "Print"], Cell[BoxData[ \("Enter ACNewtoncycle"\)], "Print"], Cell[BoxData[ InterpretationBox[ \("sprop="\[InvisibleSpace]{2, 1, 1, 1} \[InvisibleSpace]" force ="\[InvisibleSpace]{0, \(-1\)}\), SequenceForm[ "sprop=", {2, 1, 1, 1}, " force =", {0, -1}], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("entry lambdak="\[InvisibleSpace]0.106402421314340833` \[InvisibleSpace]" uk="\[InvisibleSpace]{0.`, \(-0.659609049188915008`\)} \[InvisibleSpace]" rk="\[InvisibleSpace]{0.`, 9.9515908333347447`*^-8}\), SequenceForm[ "entry lambdak=", 0.10640242131434084, " uk=", { 0.0, -0.65960904918891505}, " rk=", {0.0, 9.9515908333347447*^-08}], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("entry lambda0="\[InvisibleSpace]0.130986189484572035` \[InvisibleSpace]" u0="\[InvisibleSpace]{0.`, \(-0.519751212869467682`\)} \[InvisibleSpace]" f0="\[InvisibleSpace]1.73205080756887745`\), SequenceForm[ "entry lambda0=", 0.13098618948457202, " u0=", { 0.0, -0.51975121286946768}, " f0=", 1.7320508075688774], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("v="\[InvisibleSpace]{0.`, 4.33540534343149097`} \[InvisibleSpace]" q="\[InvisibleSpace]{0, \(-1.`\)} \[InvisibleSpace]" status="\[InvisibleSpace]" "\), SequenceForm[ "v=", {0.0, 4.335405343431491}, " q=", {0, -1.0}, " status=", " "], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("dr="\[InvisibleSpace]{0.`, 4.31441800744832892`*^-7} \[InvisibleSpace]" r="\[InvisibleSpace]{0.`, 9.9515908333347447`*^-8} \[InvisibleSpace]" status="\[InvisibleSpace]" "\), SequenceForm[ "dr=", {0.0, 4.3144180074483294*^-07}, " r=", {0.0, 9.9515908333347447*^-08}, " status=", " "], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("ck="\[InvisibleSpace]\(-1.38777878078144567`*^-17\)\), SequenceForm[ "ck=", -1.3877787807814457*^-17], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("eta="\[InvisibleSpace]\(-1.18910225649920842`*^-7\)\), SequenceForm[ "eta=", -1.1891022564992086*^-07], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("exit lambda="\[InvisibleSpace]0.106402302404115189` \[InvisibleSpace]" u="\[InvisibleSpace]{0.`, \(-0.659609133271142056`\)}\[InvisibleSpace]" r="\[InvisibleSpace]{ 0.`, 2.56739074444567449`*^-15}\), SequenceForm[ "exit lambda=", 0.10640230240411519, " u=", { 0.0, -0.65960913327114201}, " r=", {0.0, 2.5673907444456745*^-15}], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ RowBox[{ "\<\"lambda=\"\>", "\[InvisibleSpace]", "0.106402302404115189`", "\[InvisibleSpace]", "\<\" u=\"\>", "\[InvisibleSpace]", \({0.`, \(-0.659609133271142056`\)}\), "\[InvisibleSpace]", "\<\" r=\"\>", "\[InvisibleSpace]", \({0.`, 2.56739074444567449`*^-15}\), "\[InvisibleSpace]", "\<\" rnorm=\"\>", "\[InvisibleSpace]", InterpretationBox[ StyleBox["2.567390744445674*^-15", ShowStringCharacters->True, NumberMarks->True], InputForm[ 2.5673907444456745*^-15], Editable->True, AutoDelete->True]}], SequenceForm[ "lambda=", 0.10640230240411519, " u=", {0.0, -0.65960913327114201}, " r=", {0.0, 2.5673907444456745*^-15}, " rnorm=", InputForm[ 2.5673907444456745*^-15]], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("Newton converged after "\[InvisibleSpace]2 \[InvisibleSpace]" cycles"\), SequenceForm[ "Newton converged after ", 2, " cycles"], Editable->False]], "Print"], Cell[BoxData[ \("Enter ACNewtoncycle"\)], "Print"], Cell[BoxData[ InterpretationBox[ \("sprop="\[InvisibleSpace]{2, 1, 1, 1} \[InvisibleSpace]" force ="\[InvisibleSpace]{0, \(-1\)}\), SequenceForm[ "sprop=", {2, 1, 1, 1}, " force =", {0, -1}], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("entry lambdak="\[InvisibleSpace]0.083926674301550399` \[InvisibleSpace]" uk="\[InvisibleSpace]{0.`, \(-0.757050522885781163`\)} \[InvisibleSpace]" rk="\[InvisibleSpace]{0.`, 0.00310100400902273154`}\), SequenceForm[ "entry lambdak=", 0.083926674301550402, " uk=", { 0.0, -0.75705052288578112}, " rk=", {0.0, 0.0031010040090227314}], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("entry lambda0="\[InvisibleSpace]0.106402421314340833` \[InvisibleSpace]" u0="\[InvisibleSpace]{0.`, \(-0.659609049188915008`\)} \[InvisibleSpace]" f0="\[InvisibleSpace]1.73205080756887745`\), SequenceForm[ "entry lambda0=", 0.10640242131434084, " u0=", { 0.0, -0.65960904918891505}, " f0=", 1.7320508075688774], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("v="\[InvisibleSpace]{0.`, 3.43703428313873038`} \[InvisibleSpace]" q="\[InvisibleSpace]{0, \(-1.`\)} \[InvisibleSpace]" status="\[InvisibleSpace]" "\), SequenceForm[ "v=", {0.0, 3.4370342831387304}, " q=", {0, -1.0}, " status=", " "], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("dr="\[InvisibleSpace]{0.`, 0.0106582570911617713`} \[InvisibleSpace]" r="\[InvisibleSpace]{0.`, 0.00310100400902273154`} \[InvisibleSpace]" status="\[InvisibleSpace]" "\), SequenceForm[ "dr=", {0.0, 0.010658257091161772}, " r=", {0.0, 0.0031010040090227314}, " status=", " "], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[\("ck="\[InvisibleSpace]\(-0.0334157484740224131`\)\), SequenceForm[ "ck=", -0.033415748474022414], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[\("eta="\[InvisibleSpace]\(-0.0188957485579906752`\)\), SequenceForm[ "eta=", -0.018895748557990675], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("exit lambda="\[InvisibleSpace]0.0650309257435597309` \[InvisibleSpace]" u="\[InvisibleSpace]{0.`, \(-0.81133760139400266`\)}\[InvisibleSpace]" r="\[InvisibleSpace]{ 0.`, 0.000702860930741747136`}\), SequenceForm[ "exit lambda=", 0.065030925743559731, " u=", { 0.0, -0.81133760139400257}, " r=", {0.0, 0.00070286093074174716}], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ RowBox[{ "\<\"lambda=\"\>", "\[InvisibleSpace]", "0.0650309257435597309`", "\[InvisibleSpace]", "\<\" u=\"\>", "\[InvisibleSpace]", \({0.`, \(-0.81133760139400266`\)}\), "\[InvisibleSpace]", "\<\" r=\"\>", "\[InvisibleSpace]", \({0.`, 0.000702860930741747136`}\), "\[InvisibleSpace]", "\<\" rnorm=\"\>", "\[InvisibleSpace]", InterpretationBox[ StyleBox["0.0007028609307417471", ShowStringCharacters->True, NumberMarks->True], InputForm[ 0.00070286093074174716], Editable->True, AutoDelete->True]}], SequenceForm[ "lambda=", 0.065030925743559731, " u=", {0.0, -0.81133760139400257}, " r=", {0.0, 0.00070286093074174716}, " rnorm=", InputForm[ 0.00070286093074174716]], Editable->False]], "Print"], Cell[BoxData[ \("Enter ACNewtoncycle"\)], "Print"], Cell[BoxData[ InterpretationBox[ \("sprop="\[InvisibleSpace]{2, 1, 1, 1} \[InvisibleSpace]" force ="\[InvisibleSpace]{0, \(-1\)}\), SequenceForm[ "sprop=", {2, 1, 1, 1}, " force =", {0, -1}], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("entry lambdak="\[InvisibleSpace]0.0650309257435597309` \[InvisibleSpace]" uk="\[InvisibleSpace]{0.`, \(-0.81133760139400266`\)}\[InvisibleSpace]" rk="\[InvisibleSpace]{ 0.`, 0.000702860930741747136`}\), SequenceForm[ "entry lambdak=", 0.065030925743559731, " uk=", { 0.0, -0.81133760139400257}, " rk=", {0.0, 0.00070286093074174716}], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("entry lambda0="\[InvisibleSpace]0.106402421314340833` \[InvisibleSpace]" u0="\[InvisibleSpace]{0.`, \(-0.659609049188915008`\)} \[InvisibleSpace]" f0="\[InvisibleSpace]1.73205080756887745`\), SequenceForm[ "entry lambda0=", 0.10640242131434084, " u0=", { 0.0, -0.65960904918891505}, " f0=", 1.7320508075688774], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("v="\[InvisibleSpace]{0.`, 3.16655327262757735`} \[InvisibleSpace]" q="\[InvisibleSpace]{0, \(-1.`\)} \[InvisibleSpace]" status="\[InvisibleSpace]" "\), SequenceForm[ "v=", {0.0, 3.1665532726275774}, " q=", {0, -1.0}, " status=", " "], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("dr="\[InvisibleSpace]{0.`, 0.00222564658044234464`} \[InvisibleSpace]" r="\[InvisibleSpace]{0.`, 0.000702860930741747136`} \[InvisibleSpace]" status="\[InvisibleSpace]" "\), SequenceForm[ "dr=", {0.0, 0.0022256465804423443}, " r=", {0.0, 0.00070286093074174716}, " status=", " "], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[\("ck="\[InvisibleSpace]5.5511151231257827`*^-17\), SequenceForm[ "ck=", 5.5511151231257827*^-17], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("eta="\[InvisibleSpace]\(-0.000904937996491329776`\)\), SequenceForm[ "eta=", -0.00090493799649132985], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("exit lambda="\[InvisibleSpace]0.0641259877470684003` \[InvisibleSpace]" u="\[InvisibleSpace]{0.`, \(-0.811977489187875001`\)}\[InvisibleSpace]" r="\[InvisibleSpace]{ 0.`, 8.18423303849957761`*^-8}\), SequenceForm[ "exit lambda=", 0.064125987747068403, " u=", { 0.0, -0.81197748918787493}, " r=", {0.0, 8.1842330384995776*^-08}], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ RowBox[{ "\<\"lambda=\"\>", "\[InvisibleSpace]", "0.0641259877470684003`", "\[InvisibleSpace]", "\<\" u=\"\>", "\[InvisibleSpace]", \({0.`, \(-0.811977489187875001`\)}\), "\[InvisibleSpace]", "\<\" r=\"\>", "\[InvisibleSpace]", \({0.`, 8.18423303849957761`*^-8}\), "\[InvisibleSpace]", "\<\" rnorm=\"\>", "\[InvisibleSpace]", InterpretationBox[ StyleBox["8.184233038499577*^-8", ShowStringCharacters->True, NumberMarks->True], InputForm[ 8.1842330384995776*^-08], Editable->True, AutoDelete->True]}], SequenceForm[ "lambda=", 0.064125987747068403, " u=", {0.0, -0.81197748918787493}, " r=", {0.0, 8.1842330384995776*^-08}, " rnorm=", InputForm[ 8.1842330384995776*^-08]], Editable->False]], "Print"], Cell[BoxData[ \("Enter ACNewtoncycle"\)], "Print"], Cell[BoxData[ InterpretationBox[ \("sprop="\[InvisibleSpace]{2, 1, 1, 1} \[InvisibleSpace]" force ="\[InvisibleSpace]{0, \(-1\)}\), SequenceForm[ "sprop=", {2, 1, 1, 1}, " force =", {0, -1}], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("entry lambdak="\[InvisibleSpace]0.0641259877470684003` \[InvisibleSpace]" uk="\[InvisibleSpace]{0.`, \(-0.811977489187875001`\)} \[InvisibleSpace]" rk="\[InvisibleSpace]{0.`, 8.18423303849957761`*^-8}\), SequenceForm[ "entry lambdak=", 0.064125987747068403, " uk=", { 0.0, -0.81197748918787493}, " rk=", {0.0, 8.1842330384995776*^-08}], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("entry lambda0="\[InvisibleSpace]0.106402421314340833` \[InvisibleSpace]" u0="\[InvisibleSpace]{0.`, \(-0.659609049188915008`\)} \[InvisibleSpace]" f0="\[InvisibleSpace]1.73205080756887745`\), SequenceForm[ "entry lambda0=", 0.10640242131434084, " u0=", { 0.0, -0.65960904918891505}, " f0=", 1.7320508075688774], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("v="\[InvisibleSpace]{0.`, 3.1639918545828034`} \[InvisibleSpace]" q="\[InvisibleSpace]{0, \(-1.`\)} \[InvisibleSpace]" status="\[InvisibleSpace]" "\), SequenceForm[ "v=", {0.0, 3.1639918545828034}, " q=", {0, -1.0}, " status=", " "], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("dr="\[InvisibleSpace]{0.`, 2.58948466698201329`*^-7} \[InvisibleSpace]" r="\[InvisibleSpace]{0.`, 8.18423303849957761`*^-8} \[InvisibleSpace]" status="\[InvisibleSpace]" "\), SequenceForm[ "dr=", {0.0, 2.5894846669820134*^-07}, " r=", {0.0, 8.1842330384995776*^-08}, " status=", " "], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[\("ck="\[InvisibleSpace]5.5511151231257827`*^-17\), SequenceForm[ "ck=", 5.5511151231257827*^-17], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("eta="\[InvisibleSpace]\(-1.05397061276561432`*^-7\)\), SequenceForm[ "eta=", -1.0539706127656144*^-07], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("exit lambda="\[InvisibleSpace]0.0641258823500071351` \[InvisibleSpace]" u="\[InvisibleSpace]{0.`, \(-0.811977563714851591`\)}\[InvisibleSpace]" r="\[InvisibleSpace]{ 0.`, 1.11022302462515654`*^-15}\), SequenceForm[ "exit lambda=", 0.06412588235000713, " u=", { 0.0, -0.81197756371485164}, " r=", {0.0, 1.1102230246251565*^-15}], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ RowBox[{ "\<\"lambda=\"\>", "\[InvisibleSpace]", "0.0641258823500071351`", "\[InvisibleSpace]", "\<\" u=\"\>", "\[InvisibleSpace]", \({0.`, \(-0.811977563714851591`\)}\), "\[InvisibleSpace]", "\<\" r=\"\>", "\[InvisibleSpace]", \({0.`, 1.11022302462515654`*^-15}\), "\[InvisibleSpace]", "\<\" rnorm=\"\>", "\[InvisibleSpace]", InterpretationBox[ StyleBox["1.110223024625156*^-15", ShowStringCharacters->True, NumberMarks->True], InputForm[ 1.1102230246251565*^-15], Editable->True, AutoDelete->True]}], SequenceForm[ "lambda=", 0.06412588235000713, " u=", {0.0, -0.81197756371485164}, " r=", {0.0, 1.1102230246251565*^-15}, " rnorm=", InputForm[ 1.1102230246251565*^-15]], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("Newton converged after "\[InvisibleSpace]2 \[InvisibleSpace]" cycles"\), SequenceForm[ "Newton converged after ", 2, " cycles"], Editable->False]], "Print"], Cell[BoxData[ \("Enter ACNewtoncycle"\)], "Print"], Cell[BoxData[ InterpretationBox[ \("sprop="\[InvisibleSpace]{2, 1, 1, 1} \[InvisibleSpace]" force ="\[InvisibleSpace]{0, \(-1\)}\), SequenceForm[ "sprop=", {2, 1, 1, 1}, " force =", {0, -1}], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("entry lambdak="\[InvisibleSpace]0.0339897047021805454` \[InvisibleSpace]" uk="\[InvisibleSpace]{0.`, \(-0.90732844326930202`\)}\[InvisibleSpace]" rk="\[InvisibleSpace]{ 0.`, 0.00150674235795462259`}\), SequenceForm[ "entry lambdak=", 0.033989704702180545, " uk=", { 0.0, -0.90732844326930195}, " rk=", {0.0, 0.0015067423579546227}], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("entry lambda0="\[InvisibleSpace]0.0641259877470684003` \[InvisibleSpace]" u0="\[InvisibleSpace]{0.`, \(-0.811977489187875001`\)} \[InvisibleSpace]" f0="\[InvisibleSpace]1.73205080756887745`\), SequenceForm[ "entry lambda0=", 0.064125987747068403, " u0=", { 0.0, -0.81197748918787493}, " f0=", 1.7320508075688774], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("v="\[InvisibleSpace]{0.`, 2.90322599095988884`} \[InvisibleSpace]" q="\[InvisibleSpace]{0, \(-1.`\)} \[InvisibleSpace]" status="\[InvisibleSpace]" "\), SequenceForm[ "v=", {0.0, 2.9032259909598888}, " q=", {0, -1.0}, " status=", " "], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("dr="\[InvisibleSpace]{0.`, 0.00437441357529404939`} \[InvisibleSpace]" r="\[InvisibleSpace]{0.`, 0.00150674235795462259`} \[InvisibleSpace]" status="\[InvisibleSpace]" "\), SequenceForm[ "dr=", {0.0, 0.0043744135752940491}, " r=", {0.0, 0.0015067423579546227}, " status=", " "], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[\("ck="\[InvisibleSpace]\(-0.0395454631326589911`\)\), SequenceForm[ "ck=", -0.039545463132658991], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[\("eta="\[InvisibleSpace]\(-0.0240458334405596785`\)\), SequenceForm[ "eta=", -0.024045833440559678], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("exit lambda="\[InvisibleSpace]0.00994387126162086865` \[InvisibleSpace]" u="\[InvisibleSpace]{0.`, \(-0.97276451831293329`\)}\[InvisibleSpace]" r="\[InvisibleSpace]{ 0.`, 0.000321817030183000429`}\), SequenceForm[ "exit lambda=", 0.0099438712616208673, " u=", { 0.0, -0.9727645183129332}, " r=", {0.0, 0.00032181703018300048}], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ RowBox[{ "\<\"lambda=\"\>", "\[InvisibleSpace]", "0.00994387126162086865`", "\[InvisibleSpace]", "\<\" u=\"\>", "\[InvisibleSpace]", \({0.`, \(-0.97276451831293329`\)}\), "\[InvisibleSpace]", "\<\" r=\"\>", "\[InvisibleSpace]", \({0.`, 0.000321817030183000429`}\), "\[InvisibleSpace]", "\<\" rnorm=\"\>", "\[InvisibleSpace]", InterpretationBox[ StyleBox["0.0003218170301830004", ShowStringCharacters->True, NumberMarks->True], InputForm[ 0.00032181703018300048], Editable->True, AutoDelete->True]}], SequenceForm[ "lambda=", 0.0099438712616208673, " u=", {0.0, -0.9727645183129332}, " r=", {0.0, 0.00032181703018300048}, " rnorm=", InputForm[ 0.00032181703018300048]], Editable->False]], "Print"], Cell[BoxData[ \("Enter ACNewtoncycle"\)], "Print"], Cell[BoxData[ InterpretationBox[ \("sprop="\[InvisibleSpace]{2, 1, 1, 1} \[InvisibleSpace]" force ="\[InvisibleSpace]{0, \(-1\)}\), SequenceForm[ "sprop=", {2, 1, 1, 1}, " force =", {0, -1}], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("entry lambdak="\[InvisibleSpace]0.00994387126162086865` \[InvisibleSpace]" uk="\[InvisibleSpace]{0.`, \(-0.97276451831293329`\)}\[InvisibleSpace]" rk="\[InvisibleSpace]{ 0.`, 0.000321817030183000429`}\), SequenceForm[ "entry lambdak=", 0.0099438712616208673, " uk=", { 0.0, -0.9727645183129332}, " rk=", {0.0, 0.00032181703018300048}], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("entry lambda0="\[InvisibleSpace]0.0641259877470684003` \[InvisibleSpace]" u0="\[InvisibleSpace]{0.`, \(-0.811977489187875001`\)} \[InvisibleSpace]" f0="\[InvisibleSpace]1.73205080756887745`\), SequenceForm[ "entry lambda0=", 0.064125987747068403, " u0=", { 0.0, -0.81197748918787493}, " f0=", 1.7320508075688774], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("v="\[InvisibleSpace]{0.`, 2.83473530200041246`} \[InvisibleSpace]" q="\[InvisibleSpace]{0, \(-1.`\)} \[InvisibleSpace]" status="\[InvisibleSpace]" "\), SequenceForm[ "v=", {0.0, 2.8347353020004125}, " q=", {0, -1.0}, " status=", " "], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("dr="\[InvisibleSpace]{0.`, 0.000912266096244683843`} \[InvisibleSpace]" r="\[InvisibleSpace]{0.`, 0.000321817030183000429`} \[InvisibleSpace]" status="\[InvisibleSpace]" "\), SequenceForm[ "dr=", {0.0, 0.00091226609624468378}, " r=", {0.0, 0.00032181703018300048}, " status=", " "], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("ck="\[InvisibleSpace]\(-1.38777878078144567`*^-17\)\), SequenceForm[ "ck=", -1.3877787807814457*^-17], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("eta="\[InvisibleSpace]\(-0.00042877132324572269`\)\), SequenceForm[ "eta=", -0.00042877132324572277], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("exit lambda="\[InvisibleSpace]0.009515099938375144` \[InvisibleSpace]" u="\[InvisibleSpace]{0.`, \(-0.97306770542317853`\)}\[InvisibleSpace]" r="\[InvisibleSpace]{ 0.`, 2.64556391203762952`*^-9}\), SequenceForm[ "exit lambda=", 0.009515099938375144, " u=", { 0.0, -0.97306770542317855}, " r=", {0.0, 2.6455639120376295*^-09}], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ RowBox[{ "\<\"lambda=\"\>", "\[InvisibleSpace]", "0.009515099938375144`", "\[InvisibleSpace]", "\<\" u=\"\>", "\[InvisibleSpace]", \({0.`, \(-0.97306770542317853`\)}\), "\[InvisibleSpace]", "\<\" r=\"\>", "\[InvisibleSpace]", \({0.`, 2.64556391203762952`*^-9}\), "\[InvisibleSpace]", "\<\" rnorm=\"\>", "\[InvisibleSpace]", InterpretationBox[ StyleBox["2.645563912037629*^-9", ShowStringCharacters->True, NumberMarks->True], InputForm[ 2.6455639120376295*^-09], Editable->True, AutoDelete->True]}], SequenceForm[ "lambda=", 0.009515099938375144, " u=", {0.0, -0.97306770542317855}, " r=", {0.0, 2.6455639120376295*^-09}, " rnorm=", InputForm[ 2.6455639120376295*^-09]], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("Newton converged after "\[InvisibleSpace]1 \[InvisibleSpace]" cycles"\), SequenceForm[ "Newton converged after ", 1, " cycles"], Editable->False]], "Print"], Cell[BoxData[ \("Enter ACNewtoncycle"\)], "Print"], Cell[BoxData[ InterpretationBox[ \("sprop="\[InvisibleSpace]{2, 1, 1, 1} \[InvisibleSpace]" force ="\[InvisibleSpace]{0, \(-1\)}\), SequenceForm[ "sprop=", {2, 1, 1, 1}, " force =", {0, -1}], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("entry lambdak="\[InvisibleSpace]\(-0.0233235023842155975` \)\[InvisibleSpace]" uk="\[InvisibleSpace]{0.`, \(-1.06706871679162396`\)}\[InvisibleSpace]" rk="\[InvisibleSpace]{ 0.`, 0.000282206544607012368`}\), SequenceForm[ "entry lambdak=", -0.023323502384215596, " uk=", { 0.0, -1.067068716791624}, " rk=", {0.0, 0.00028220654460701239}], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("entry lambda0="\[InvisibleSpace]0.00994387126162086865` \[InvisibleSpace]" u0="\[InvisibleSpace]{0.`, \(-0.97276451831293329`\)} \[InvisibleSpace]" f0="\[InvisibleSpace]1.73205080756887745`\), SequenceForm[ "entry lambda0=", 0.0099438712616208673, " u0=", { 0.0, -0.9727645183129332}, " f0=", 1.7320508075688774], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("v="\[InvisibleSpace]{0.`, 2.86711784304812091`} \[InvisibleSpace]" q="\[InvisibleSpace]{0, \(-1.`\)} \[InvisibleSpace]" status="\[InvisibleSpace]" "\), SequenceForm[ "v=", {0.0, 2.867117843048121}, " q=", {0, -1.0}, " status=", " "], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("dr="\[InvisibleSpace]{0.`, 0.000809119419467720746`} \[InvisibleSpace]" r="\[InvisibleSpace]{0.`, 0.000282206544607012368`} \[InvisibleSpace]" status="\[InvisibleSpace]" "\), SequenceForm[ "dr=", {0.0, 0.00080911941946772072}, " r=", {0.0, 0.00028220654460701239}, " status=", " "], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[\("ck="\[InvisibleSpace]\(-0.0422078715046748698`\)\), SequenceForm[ "ck=", -0.042207871504674872], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[\("eta="\[InvisibleSpace]\(-0.0243068169900595476`\)\), SequenceForm[ "eta=", -0.02430681699005955], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("exit lambda="\[InvisibleSpace]\(-0.0476303193742751407` \)\[InvisibleSpace]" u="\[InvisibleSpace]{0.`, \(-1.13595010607206114`\)}\[InvisibleSpace]" r="\[InvisibleSpace]{ 0.`, \(-0.000453067960114697498`\)}\), SequenceForm[ "exit lambda=", -0.047630319374275143, " u=", { 0.0, -1.1359501060720611}, " r=", {0.0, -0.00045306796011469752}], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ RowBox[{ "\<\"lambda=\"\>", "\[InvisibleSpace]", \(-0.0476303193742751407`\), "\[InvisibleSpace]", "\<\" u=\"\>", "\[InvisibleSpace]", \({0.`, \(-1.13595010607206114`\)}\), "\[InvisibleSpace]", "\<\" r=\"\>", "\[InvisibleSpace]", \({0.`, \(-0.000453067960114697498`\)}\), "\[InvisibleSpace]", "\<\" rnorm=\"\>", "\[InvisibleSpace]", InterpretationBox[ StyleBox["0.0004530679601146974", ShowStringCharacters->True, NumberMarks->True], InputForm[ 0.00045306796011469752], Editable->True, AutoDelete->True]}], SequenceForm[ "lambda=", -0.047630319374275143, " u=", {0.0, -1.1359501060720611}, " r=", {0.0, -0.00045306796011469752}, " rnorm=", InputForm[ 0.00045306796011469752]], Editable->False]], "Print"], Cell[BoxData[ \("Enter ACNewtoncycle"\)], "Print"], Cell[BoxData[ InterpretationBox[ \("sprop="\[InvisibleSpace]{2, 1, 1, 1} \[InvisibleSpace]" force ="\[InvisibleSpace]{0, \(-1\)}\), SequenceForm[ "sprop=", {2, 1, 1, 1}, " force =", {0, -1}], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("entry lambdak="\[InvisibleSpace]\(-0.0476303193742751407` \)\[InvisibleSpace]" uk="\[InvisibleSpace]{0.`, \(-1.13595010607206114`\)}\[InvisibleSpace]" rk="\[InvisibleSpace]{ 0.`, \(-0.000453067960114697498`\)}\), SequenceForm[ "entry lambdak=", -0.047630319374275143, " uk=", { 0.0, -1.1359501060720611}, " rk=", {0.0, -0.00045306796011469752}], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("entry lambda0="\[InvisibleSpace]0.00994387126162086865` \[InvisibleSpace]" u0="\[InvisibleSpace]{0.`, \(-0.97276451831293329`\)} \[InvisibleSpace]" f0="\[InvisibleSpace]1.73205080756887745`\), SequenceForm[ "entry lambda0=", 0.0099438712616208673, " u0=", { 0.0, -0.9727645183129332}, " f0=", 1.7320508075688774], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("v="\[InvisibleSpace]{0.`, 2.99446193616124602`} \[InvisibleSpace]" q="\[InvisibleSpace]{0, \(-1.`\)} \[InvisibleSpace]" status="\[InvisibleSpace]" "\), SequenceForm[ "v=", {0.0, 2.9944619361612461}, " q=", {0, -1.0}, " status=", " "], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("dr="\[InvisibleSpace]{0.`, \(-0.0013566947610576836`\)} \[InvisibleSpace]" r="\[InvisibleSpace]{0.`, \(-0.000453067960114697498`\)} \[InvisibleSpace]" status="\[InvisibleSpace]" "\), SequenceForm[ "dr=", {0.0, -0.0013566947610576835}, " r=", { 0.0, -0.00045306796011469752}, " status=", " "], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("ck="\[InvisibleSpace]\(-4.16333634234433702`*^-17\)\), SequenceForm[ "ck=", -4.163336342344337*^-17], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[\("eta="\[InvisibleSpace]0.000593128162938316982`\), SequenceForm[ "eta=", 0.00059312816293831698], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("exit lambda="\[InvisibleSpace]\(-0.0470371912113368217` \)\[InvisibleSpace]" u="\[InvisibleSpace]{0.`, \(-1.13553070112593479`\)}\[InvisibleSpace]" r="\[InvisibleSpace]{ 0.`, \(-2.53382186690687127`*^-8\)}\), SequenceForm[ "exit lambda=", -0.047037191211336829, " u=", { 0.0, -1.1355307011259348}, " r=", {0.0, -2.5338218669068713*^-08}], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ RowBox[{ "\<\"lambda=\"\>", "\[InvisibleSpace]", \(-0.0470371912113368217`\), "\[InvisibleSpace]", "\<\" u=\"\>", "\[InvisibleSpace]", \({0.`, \(-1.13553070112593479`\)}\), "\[InvisibleSpace]", "\<\" r=\"\>", "\[InvisibleSpace]", \({0.`, \(-2.53382186690687127`*^-8\)}\), "\[InvisibleSpace]", "\<\" rnorm=\"\>", "\[InvisibleSpace]", InterpretationBox[ StyleBox["2.533821866906871*^-8", ShowStringCharacters->True, NumberMarks->True], InputForm[ 2.5338218669068713*^-08], Editable->True, AutoDelete->True]}], SequenceForm[ "lambda=", -0.047037191211336829, " u=", {0.0, -1.1355307011259348}, " r=", {0.0, -2.5338218669068713*^-08}, " rnorm=", InputForm[ 2.5338218669068713*^-08]], Editable->False]], "Print"], Cell[BoxData[ \("Enter ACNewtoncycle"\)], "Print"], Cell[BoxData[ InterpretationBox[ \("sprop="\[InvisibleSpace]{2, 1, 1, 1} \[InvisibleSpace]" force ="\[InvisibleSpace]{0, \(-1\)}\), SequenceForm[ "sprop=", {2, 1, 1, 1}, " force =", {0, -1}], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("entry lambdak="\[InvisibleSpace]\(-0.0470371912113368217` \)\[InvisibleSpace]" uk="\[InvisibleSpace]{0.`, \(-1.13553070112593479`\)}\[InvisibleSpace]" rk="\[InvisibleSpace]{ 0.`, \(-2.53382186690687127`*^-8\)}\), SequenceForm[ "entry lambdak=", -0.047037191211336829, " uk=", { 0.0, -1.1355307011259348}, " rk=", {0.0, -2.5338218669068713*^-08}], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("entry lambda0="\[InvisibleSpace]0.00994387126162086865` \[InvisibleSpace]" u0="\[InvisibleSpace]{0.`, \(-0.97276451831293329`\)} \[InvisibleSpace]" f0="\[InvisibleSpace]1.73205080756887745`\), SequenceForm[ "entry lambda0=", 0.0099438712616208673, " u0=", { 0.0, -0.9727645183129332}, " f0=", 1.7320508075688774], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("v="\[InvisibleSpace]{0.`, 2.99337943216542825`} \[InvisibleSpace]" q="\[InvisibleSpace]{0, \(-1.`\)} \[InvisibleSpace]" status="\[InvisibleSpace]" "\), SequenceForm[ "v=", {0.0, 2.9933794321654283}, " q=", {0, -1.0}, " status=", " "], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("dr="\[InvisibleSpace]{0.`, \(-7.58469026117003597`*^-8\)} \[InvisibleSpace]" r="\[InvisibleSpace]{0.`, \(-2.53382186690687127`*^-8\)} \[InvisibleSpace]" status="\[InvisibleSpace]" "\), SequenceForm[ "dr=", {0.0, -7.584690261170036*^-08}, " r=", { 0.0, -2.5338218669068713*^-08}, " status=", " "], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("ck="\[InvisibleSpace]\(-1.38777878078144567`*^-17\)\), SequenceForm[ "ck=", -1.3877787807814457*^-17], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[\("eta="\[InvisibleSpace]3.31749157574138475`*^-8\), SequenceForm[ "eta=", 3.3174915757413844*^-08], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("exit lambda="\[InvisibleSpace]\(-0.047037158036421065` \)\[InvisibleSpace]" u="\[InvisibleSpace]{0.`, \(-1.13553067766772675`\)}\[InvisibleSpace]" r="\[InvisibleSpace]{ 0.`, \(-1.2490009027033011`*^-16\)}\), SequenceForm[ "exit lambda=", -0.04703715803642107, " u=", { 0.0, -1.1355306776677268}, " r=", {0.0, -1.2490009027033011*^-16}], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ RowBox[{ "\<\"lambda=\"\>", "\[InvisibleSpace]", \(-0.047037158036421065`\), "\[InvisibleSpace]", "\<\" u=\"\>", "\[InvisibleSpace]", \({0.`, \(-1.13553067766772675`\)}\), "\[InvisibleSpace]", "\<\" r=\"\>", "\[InvisibleSpace]", \({0.`, \(-1.2490009027033011`*^-16\)}\), "\[InvisibleSpace]", "\<\" rnorm=\"\>", "\[InvisibleSpace]", InterpretationBox[ StyleBox["1.249000902703301*^-16", ShowStringCharacters->True, NumberMarks->True], InputForm[ 1.2490009027033011*^-16], Editable->True, AutoDelete->True]}], SequenceForm[ "lambda=", -0.04703715803642107, " u=", {0.0, -1.1355306776677268}, " r=", {0.0, -1.2490009027033011*^-16}, " rnorm=", InputForm[ 1.2490009027033011*^-16]], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("Newton converged after "\[InvisibleSpace]2 \[InvisibleSpace]" cycles"\), SequenceForm[ "Newton converged after ", 2, " cycles"], Editable->False]], "Print"], Cell[BoxData[ \("Enter ACNewtoncycle"\)], "Print"], Cell[BoxData[ InterpretationBox[ \("sprop="\[InvisibleSpace]{2, 1, 1, 1} \[InvisibleSpace]" force ="\[InvisibleSpace]{0, \(-1\)}\), SequenceForm[ "sprop=", {2, 1, 1, 1}, " force =", {0, -1}], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("entry lambdak="\[InvisibleSpace]\(-0.0787228940585015202` \)\[InvisibleSpace]" uk="\[InvisibleSpace]{0.`, \(-1.23037803232234321`\)}\[InvisibleSpace]" rk="\[InvisibleSpace]{ 0.`, \(-0.00159488958445384865`\)}\), SequenceForm[ "entry lambdak=", -0.078722894058501525, " uk=", { 0.0, -1.2303780323223432}, " rk=", {0.0, -0.0015948895844538485}], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("entry lambda0="\[InvisibleSpace]\(-0.0470371912113368217` \)\[InvisibleSpace]" u0="\[InvisibleSpace]{0.`, \(-1.13553070112593479`\)} \[InvisibleSpace]" f0="\[InvisibleSpace]1.73205080756887745`\), SequenceForm[ "entry lambda0=", -0.047037191211336829, " u0=", { 0.0, -1.1355307011259348}, " f0=", 1.7320508075688774], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("v="\[InvisibleSpace]{0.`, 3.36405984218851328`} \[InvisibleSpace]" q="\[InvisibleSpace]{0, \(-1.`\)} \[InvisibleSpace]" status="\[InvisibleSpace]" "\), SequenceForm[ "v=", {0.0, 3.3640598421885133}, " q=", {0, -1.0}, " status=", " "], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("dr="\[InvisibleSpace]{0.`, \(-0.00536530400378591743`\)} \[InvisibleSpace]" r="\[InvisibleSpace]{0.`, \(-0.00159488958445384865`\)} \[InvisibleSpace]" status="\[InvisibleSpace]" "\), SequenceForm[ "dr=", {0.0, -0.0053653040037859173}, " r=", { 0.0, -0.0015948895844538485}, " status=", " "], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[\("ck="\[InvisibleSpace]\(-0.0408512274362859084`\)\), SequenceForm[ "ck=", -0.040851227436285907], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[\("eta="\[InvisibleSpace]\(-0.0168113723740426905`\)\), SequenceForm[ "eta=", -0.016811372374042691], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("exit lambda="\[InvisibleSpace]\(-0.0955342664325442214` \)\[InvisibleSpace]" u="\[InvisibleSpace]{0.`, \(-1.29229779902172348`\)}\[InvisibleSpace]" r="\[InvisibleSpace]{ 0.`, \(-0.00102079767045698543`\)}\), SequenceForm[ "exit lambda=", -0.095534266432544215, " u=", { 0.0, -1.2922977990217235}, " r=", {0.0, -0.0010207976704569854}], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ RowBox[{ "\<\"lambda=\"\>", "\[InvisibleSpace]", \(-0.0955342664325442214`\), "\[InvisibleSpace]", "\<\" u=\"\>", "\[InvisibleSpace]", \({0.`, \(-1.29229779902172348`\)}\), "\[InvisibleSpace]", "\<\" r=\"\>", "\[InvisibleSpace]", \({0.`, \(-0.00102079767045698543`\)}\), "\[InvisibleSpace]", "\<\" rnorm=\"\>", "\[InvisibleSpace]", InterpretationBox[ StyleBox["0.001020797670456985", ShowStringCharacters->True, NumberMarks->True], InputForm[ 0.0010207976704569854], Editable->True, AutoDelete->True]}], SequenceForm[ "lambda=", -0.095534266432544215, " u=", {0.0, -1.2922977990217235}, " r=", {0.0, -0.0010207976704569854}, " rnorm=", InputForm[ 0.0010207976704569854]], Editable->False]], "Print"], Cell[BoxData[ \("Enter ACNewtoncycle"\)], "Print"], Cell[BoxData[ InterpretationBox[ \("sprop="\[InvisibleSpace]{2, 1, 1, 1} \[InvisibleSpace]" force ="\[InvisibleSpace]{0, \(-1\)}\), SequenceForm[ "sprop=", {2, 1, 1, 1}, " force =", {0, -1}], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("entry lambdak="\[InvisibleSpace]\(-0.0955342664325442214` \)\[InvisibleSpace]" uk="\[InvisibleSpace]{0.`, \(-1.29229779902172348`\)}\[InvisibleSpace]" rk="\[InvisibleSpace]{ 0.`, \(-0.00102079767045698543`\)}\), SequenceForm[ "entry lambdak=", -0.095534266432544215, " uk=", { 0.0, -1.2922977990217235}, " rk=", {0.0, -0.0010207976704569854}], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("entry lambda0="\[InvisibleSpace]\(-0.0470371912113368217` \)\[InvisibleSpace]" u0="\[InvisibleSpace]{0.`, \(-1.13553070112593479`\)} \[InvisibleSpace]" f0="\[InvisibleSpace]1.73205080756887745`\), SequenceForm[ "entry lambda0=", -0.047037191211336829, " u0=", { 0.0, -1.1355307011259348}, " f0=", 1.7320508075688774], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("v="\[InvisibleSpace]{0.`, 3.8032545490248534`} \[InvisibleSpace]" q="\[InvisibleSpace]{0, \(-1.`\)} \[InvisibleSpace]" status="\[InvisibleSpace]" "\), SequenceForm[ "v=", {0.0, 3.8032545490248535}, " q=", {0, -1.0}, " status=", " "], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("dr="\[InvisibleSpace]{0.`, \(-0.00388235338379950345`\)} \[InvisibleSpace]" r="\[InvisibleSpace]{0.`, \(-0.00102079767045698543`\)} \[InvisibleSpace]" status="\[InvisibleSpace]" "\), SequenceForm[ "dr=", {0.0, -0.0038823533837995033}, " r=", { 0.0, -0.0010207976704569854}, " status=", " "], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("ck="\[InvisibleSpace]\(-4.16333634234433702`*^-17\)\), SequenceForm[ "ck=", -4.163336342344337*^-17], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[\("eta="\[InvisibleSpace]0.00125393026267285234`\), SequenceForm[ "eta=", 0.0012539302626728523], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("exit lambda="\[InvisibleSpace]\(-0.0942803361698713615` \)\[InvisibleSpace]" u="\[InvisibleSpace]{0.`, \(-1.2914111364298526`\)}\[InvisibleSpace]" r="\[InvisibleSpace]{ 0.`, \(-2.43488931267465957`*^-7\)}\), SequenceForm[ "exit lambda=", -0.094280336169871359, " u=", { 0.0, -1.2914111364298526}, " r=", {0.0, -2.4348893126746596*^-07}], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ RowBox[{ "\<\"lambda=\"\>", "\[InvisibleSpace]", \(-0.0942803361698713615`\), "\[InvisibleSpace]", "\<\" u=\"\>", "\[InvisibleSpace]", \({0.`, \(-1.2914111364298526`\)}\), "\[InvisibleSpace]", "\<\" r=\"\>", "\[InvisibleSpace]", \({0.`, \(-2.43488931267465957`*^-7\)}\), "\[InvisibleSpace]", "\<\" rnorm=\"\>", "\[InvisibleSpace]", InterpretationBox[ StyleBox["2.434889312674659*^-7", ShowStringCharacters->True, NumberMarks->True], InputForm[ 2.4348893126746596*^-07], Editable->True, AutoDelete->True]}], SequenceForm[ "lambda=", -0.094280336169871359, " u=", {0.0, -1.2914111364298526}, " r=", {0.0, -2.4348893126746596*^-07}, " rnorm=", InputForm[ 2.4348893126746596*^-07]], Editable->False]], "Print"], Cell[BoxData[ \("Enter ACNewtoncycle"\)], "Print"], Cell[BoxData[ InterpretationBox[ \("sprop="\[InvisibleSpace]{2, 1, 1, 1} \[InvisibleSpace]" force ="\[InvisibleSpace]{0, \(-1\)}\), SequenceForm[ "sprop=", {2, 1, 1, 1}, " force =", {0, -1}], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("entry lambdak="\[InvisibleSpace]\(-0.0942803361698713615` \)\[InvisibleSpace]" uk="\[InvisibleSpace]{0.`, \(-1.2914111364298526`\)}\[InvisibleSpace]" rk="\[InvisibleSpace]{ 0.`, \(-2.43488931267465957`*^-7\)}\), SequenceForm[ "entry lambdak=", -0.094280336169871359, " uk=", { 0.0, -1.2914111364298526}, " rk=", {0.0, -2.4348893126746596*^-07}], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("entry lambda0="\[InvisibleSpace]\(-0.0470371912113368217` \)\[InvisibleSpace]" u0="\[InvisibleSpace]{0.`, \(-1.13553070112593479`\)} \[InvisibleSpace]" f0="\[InvisibleSpace]1.73205080756887745`\), SequenceForm[ "entry lambda0=", -0.047037191211336829, " u0=", { 0.0, -1.1355307011259348}, " f0=", 1.7320508075688774], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("v="\[InvisibleSpace]{0.`, 3.79533070339675049`} \[InvisibleSpace]" q="\[InvisibleSpace]{0, \(-1.`\)} \[InvisibleSpace]" status="\[InvisibleSpace]" "\), SequenceForm[ "v=", {0.0, 3.7953307033967505}, " q=", {0, -1.0}, " status=", " "], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("dr="\[InvisibleSpace]{0.`, \(-9.24121016776674508`*^-7\)} \[InvisibleSpace]" r="\[InvisibleSpace]{0.`, \(-2.43488931267465957`*^-7\)} \[InvisibleSpace]" status="\[InvisibleSpace]" "\), SequenceForm[ "dr=", {0.0, -9.2412101677667455*^-07}, " r=", { 0.0, -2.4348893126746596*^-07}, " status=", " "], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[\("ck="\[InvisibleSpace]8.32667268468867405`*^-17\), SequenceForm[ "ck=", 8.3266726846886741*^-17], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[\("eta="\[InvisibleSpace]2.99240288319919844`*^-7\), SequenceForm[ "eta=", 2.9924028831991985*^-07], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("exit lambda="\[InvisibleSpace]\(-0.0942800369295830265` \)\[InvisibleSpace]" u="\[InvisibleSpace]{0.`, \(-1.29141092483501541`\)}\[InvisibleSpace]" r="\[InvisibleSpace]{ 0.`, \(-1.38500322321988278`*^-14\)}\), SequenceForm[ "exit lambda=", -0.094280036929583033, " u=", { 0.0, -1.2914109248350154}, " r=", {0.0, -1.3850032232198828*^-14}], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ RowBox[{ "\<\"lambda=\"\>", "\[InvisibleSpace]", \(-0.0942800369295830265`\), "\[InvisibleSpace]", "\<\" u=\"\>", "\[InvisibleSpace]", \({0.`, \(-1.29141092483501541`\)}\), "\[InvisibleSpace]", "\<\" r=\"\>", "\[InvisibleSpace]", \({0.`, \(-1.38500322321988278`*^-14\)}\), "\[InvisibleSpace]", "\<\" rnorm=\"\>", "\[InvisibleSpace]", InterpretationBox[ StyleBox["1.385003223219882*^-14", ShowStringCharacters->True, NumberMarks->True], InputForm[ 1.3850032232198829*^-14], Editable->True, AutoDelete->True]}], SequenceForm[ "lambda=", -0.094280036929583033, " u=", {0.0, -1.2914109248350154}, " r=", {0.0, -1.3850032232198828*^-14}, " rnorm=", InputForm[ 1.3850032232198829*^-14]], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ \("Newton converged after "\[InvisibleSpace]2 \[InvisibleSpace]" cycles"\), SequenceForm[ "Newton converged after ", 2, " cycles"], Editable->False]], "Print"], Cell[BoxData[ TagBox[GridBox[{ {"\<\"n\"\>", "\<\"lambda\"\>", "\<\"u={uX,uY}\"\>", "\<\"v={vX,vY}\"\>", "\<\"Kdet\"\>", "\<\"a\"\>", "\<\"kappa\"\>", "\<\"remarks\"\>"}, {"\<\"0\"\>", "\<\"0\"\>", "\<\"{0, 0}\"\>", "\<\"{0, -1.41421}\"\>", "\<\"0.5\"\>", "\<\"1\"\>", "\<\"1\"\>", "\<\"C0 cfg\"\>"}, {"\<\"1\"\>", "\<\"0.0527596\"\>", "\<\"{0, -0.0851678}\"\>", "\<\"{0, -1.41421}\"\>", "\<\"0.352568\"\>", "\<\"1.\"\>", "\<\"1.\"\>", "\<\" \"\>"}, {"\<\"2\"\>", "\<\"0.0943566\"\>", "\<\"{0, -0.178229}\"\>", "\<\"{0, -1.8722}\"\>", "\<\"0.22092\"\>", "\<\"1.\"\>", "\<\"0.755377\"\>", "\<\" \"\>"}, {"\<\"3\"\>", "\<\"0.122769\"\>", "\<\"{0, -0.280612}\"\>", "\<\"{0, -2.75696}\"\>", "\<\"0.112605\"\>", "\<\"1.\"\>", "\<\"0.512962\"\>", "\<\" \"\>"}, {"\<\"4\"\>", "\<\"0.13572\"\>", "\<\"{0, -0.39393}\"\>", "\<\"{0, -5.11881}\"\>", "\<\"0.0273425\"\>", "\<\"1.\"\>", "\<\"0.276278\"\>", "\<\" \"\>"}, {"\<\"5\"\>", "\<\"0.130986\"\>", "\<\"{0, -0.519751}\"\>", "\<\"{0, -27.7395}\"\>", "\<\"-0.0363475\"\>", "\<\"1.\"\>", "\<\"0.050982\"\>", "\<\" \"\>"}, {"\<\"6\"\>", "\<\"0.106402\"\>", "\<\"{0, -0.659609}\"\>", "\<\"{0, 9.18072}\"\>", "\<\"-0.0808525\"\>", "\<\"-1.\"\>", "\<\"-0.154042\"\>", "\<\" \"\>"}, {"\<\"7\"\>", "\<\"0.064126\"\>", "\<\"{0, -0.811977}\"\>", "\<\"{0, 4.33541}\"\>", "\<\"-0.108937\"\>", "\<\"1.\"\>", "\<\"-0.326201\"\>", "\<\" \"\>"}, {"\<\"8\"\>", "\<\"0.00994387\"\>", "\<\"{0, -0.972765}\"\>", "\<\"{0, 3.16399}\"\>", "\<\"-0.122825\"\>", "\<\"1.\"\>", "\<\"-0.446971\"\>", "\<\" \"\>"}, {"\<\"9\"\>", "\<\"-0.0470372\"\>", "\<\"{0, -1.13553}\"\>", "\<\"{0, 2.83474}\"\>", "\<\"-0.123868\"\>", "\<\"1.\"\>", "\<\"-0.498887\"\>", "\<\" \"\>"} }, RowSpacings->0, ColumnSpacings->2, RowAlignments->Baseline, ColumnAlignments->{Left}], (TableForm[ #, TableAlignments -> {Left}, TableDirections -> {Column, Row}, TableSpacing -> {0, 2}]&)]], "Print"], Cell[BoxData[ \(General::"spell1" \( : \ \) "Possible spelling error: new symbol name \"\!\(uYvslambda\)\" is \ similar to existing symbol \"\!\(uXvslambda\)\"."\)], "Message"], Cell[GraphicsData["PostScript", "\<\ %! %%Creator: Mathematica %%AspectRatio: .61803 MathPictureStart /Mabs { Mgmatrix idtransform Mtmatrix dtransform } bind def /Mabsadd { Mabs 3 -1 roll add 3 1 roll add exch } bind def %% Graphics %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10 scalefont setfont % Scaling calculations 0.97619 0.737473 0.255992 2.55915 [ [.09122 .24349 -12 -9 ] [.09122 .24349 12 0 ] [.23872 .24349 -6 -9 ] [.23872 .24349 6 0 ] [.38621 .24349 -12 -9 ] [.38621 .24349 12 0 ] [.53371 .24349 -12 -9 ] [.53371 .24349 12 0 ] [.6812 .24349 -12 -9 ] [.6812 .24349 12 0 ] [.8287 .24349 -12 -9 ] [.8287 .24349 12 0 ] [1.025 .25599 0 -6 ] [1.025 .25599 16 6 ] [.96369 8e-05 -24 -4.5 ] [.96369 8e-05 0 4.5 ] [.96369 .12803 -30 -4.5 ] [.96369 .12803 0 4.5 ] [.96369 .38395 -24 -4.5 ] [.96369 .38395 0 4.5 ] [.96369 .51191 -18 -4.5 ] [.96369 .51191 0 4.5 ] [.97619 .64303 -20 0 ] [.97619 .64303 20 12 ] [.5 .63053 -68 0 ] [.5 .63053 68 12 ] [ 0 0 0 0 ] [ 1 .61803 0 0 ] ] MathScale % Start of Graphics 1 setlinecap 1 setlinejoin newpath 0 g .25 Mabswid [ ] 0 setdash .09122 .25599 m .09122 .26224 L s [(-1.2)] .09122 .24349 0 1 Mshowa .23872 .25599 m .23872 .26224 L s [(-1)] .23872 .24349 0 1 Mshowa .38621 .25599 m .38621 .26224 L s [(-0.8)] .38621 .24349 0 1 Mshowa .53371 .25599 m .53371 .26224 L s [(-0.6)] .53371 .24349 0 1 Mshowa .6812 .25599 m .6812 .26224 L s [(-0.4)] .6812 .24349 0 1 Mshowa .8287 .25599 m .8287 .26224 L s [(-0.2)] .8287 .24349 0 1 Mshowa .125 Mabswid .1281 .25599 m .1281 .25974 L s .16497 .25599 m .16497 .25974 L s .20184 .25599 m .20184 .25974 L s .27559 .25599 m .27559 .25974 L s .31246 .25599 m .31246 .25974 L s .34934 .25599 m .34934 .25974 L s .42309 .25599 m .42309 .25974 L s .45996 .25599 m .45996 .25974 L s .49683 .25599 m .49683 .25974 L s .57058 .25599 m .57058 .25974 L s .60745 .25599 m .60745 .25974 L s .64433 .25599 m .64433 .25974 L s .71807 .25599 m .71807 .25974 L s .75495 .25599 m .75495 .25974 L s .79182 .25599 m .79182 .25974 L s .86557 .25599 m .86557 .25974 L s .90244 .25599 m .90244 .25974 L s .93932 .25599 m .93932 .25974 L s .05435 .25599 m .05435 .25974 L s .01748 .25599 m .01748 .25974 L s .25 Mabswid 0 .25599 m 1 .25599 L s gsave 1.025 .25599 -61 -10 Mabsadd m 1 1 Mabs scale currentpoint translate 0 20 translate 1 -1 scale gsave 0.000000 0.000000 0.000000 setrgbcolor 1.000000 setlinewidth gsave newpath 61.000000 16.000000 moveto 558.000000 16.000000 lineto 558.000000 4.000000 lineto 61.000000 4.000000 lineto 61.000000 16.000000 lineto closepath clip newpath 63.000000 13.000000 moveto %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10.000000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000000 0.000000 0.000000 setrgbcolor 0.000000 0.000000 rmoveto 63.000000 13.000000 moveto %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10.000000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000000 0.000000 0.000000 setrgbcolor (uY) show 75.000000 13.000000 moveto %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10.000000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000000 0.000000 0.000000 setrgbcolor 0.000000 0.000000 rmoveto 1.000000 setlinewidth grestore grestore %%DocumentNeededResources: font Courier %%DocumentSuppliedResources: %%DocumentNeededFonts: Courier %%DocumentSuppliedFonts: %%DocumentFonts: font Courier grestore .97619 8e-05 m .98244 8e-05 L s [(-0.1)] .96369 8e-05 1 0 Mshowa .97619 .12803 m .98244 .12803 L s [(-0.05)] .96369 .12803 1 0 Mshowa .97619 .38395 m .98244 .38395 L s [(0.05)] .96369 .38395 1 0 Mshowa .97619 .51191 m .98244 .51191 L s [(0.1)] .96369 .51191 1 0 Mshowa .125 Mabswid .97619 .02567 m .97994 .02567 L s .97619 .05126 m .97994 .05126 L s .97619 .07685 m .97994 .07685 L s .97619 .10244 m .97994 .10244 L s .97619 .15363 m .97994 .15363 L s .97619 .17922 m .97994 .17922 L s .97619 .20481 m .97994 .20481 L s .97619 .2304 m .97994 .2304 L s .97619 .28158 m .97994 .28158 L s .97619 .30718 m .97994 .30718 L s .97619 .33277 m .97994 .33277 L s .97619 .35836 m .97994 .35836 L s .97619 .40954 m .97994 .40954 L s .97619 .43513 m .97994 .43513 L s .97619 .46072 m .97994 .46072 L s .97619 .48632 m .97994 .48632 L s .97619 .5375 m .97994 .5375 L s .97619 .56309 m .97994 .56309 L s .97619 .58868 m .97994 .58868 L s .97619 .61427 m .97994 .61427 L s .25 Mabswid .97619 0 m .97619 .61803 L s gsave .97619 .64303 -81 -4 Mabsadd m 1 1 Mabs scale currentpoint translate 0 20 translate 1 -1 scale gsave 0.000000 0.000000 0.000000 setrgbcolor 1.000000 setlinewidth gsave newpath 61.000000 16.000000 moveto 558.000000 16.000000 lineto 558.000000 4.000000 lineto 61.000000 4.000000 lineto 61.000000 16.000000 lineto closepath clip newpath 63.000000 13.000000 moveto %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10.000000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000000 0.000000 0.000000 setrgbcolor 0.000000 0.000000 rmoveto 63.000000 13.000000 moveto %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10.000000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000000 0.000000 0.000000 setrgbcolor (lambda) show 99.000000 13.000000 moveto %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10.000000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000000 0.000000 0.000000 setrgbcolor 0.000000 0.000000 rmoveto 1.000000 setlinewidth grestore grestore %%DocumentNeededResources: font Courier %%DocumentSuppliedResources: %%DocumentNeededFonts: Courier %%DocumentSuppliedFonts: %%DocumentFonts: font Courier grestore gsave .5 .63053 -129 -4 Mabsadd m 1 1 Mabs scale currentpoint translate 0 20 translate 1 -1 scale gsave 0.000000 0.000000 0.000000 setrgbcolor 1.000000 setlinewidth gsave newpath 61.000000 16.000000 moveto 558.000000 16.000000 lineto 558.000000 4.000000 lineto 61.000000 4.000000 lineto 61.000000 16.000000 lineto closepath clip newpath 63.000000 13.000000 moveto %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10.000000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000000 0.000000 0.000000 setrgbcolor 0.000000 0.000000 rmoveto 63.000000 13.000000 moveto %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10.000000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000000 0.000000 0.000000 setrgbcolor (Response) show 111.000000 13.000000 moveto %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10.000000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000000 0.000000 0.000000 setrgbcolor (:) show 123.000000 13.000000 moveto %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10.000000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000000 0.000000 0.000000 setrgbcolor (lambda) show 165.000000 13.000000 moveto (vs) show 183.000000 13.000000 moveto (uY) show 195.000000 13.000000 moveto %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10.000000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000000 0.000000 0.000000 setrgbcolor 0.000000 0.000000 rmoveto 1.000000 setlinewidth grestore grestore %%DocumentNeededResources: font Courier %%DocumentSuppliedResources: %%DocumentNeededFonts: Courier %%DocumentSuppliedFonts: %%DocumentFonts: font Courier grestore .5 Mabswid .97619 .25599 m .91338 .39101 L .84475 .49746 L .76925 .57018 L .68568 .60332 L .59289 .59121 L .48975 .52829 L .37738 .4201 L .2588 .28144 L .13877 .13562 L .02381 .01472 L s 0 0 m 1 0 L 1 .61803 L 0 .61803 L closepath clip newpath % End of Graphics MathPictureEnd \ \>"], "Graphics", Evaluatable->False, AspectRatioFixed->True, ImageSize->{282, 174}, ImageMargins->{{34, 0}, {0, 0}}, ImageRegion->{{0, 1}, {0, 1}}, Background->RGBColor[0, 1, 1], ImageCache->GraphicsData["Bitmap", "\<\ CF5dJ6E]HGAYHf4PAg9QL6QYHg400?oo00@00000003oo`00ool000000`00ool2000000@00?oo00<00000003oo`00ool0 :@00ool00=/00?oo1@0000001@00ool00000003oo`00ool0000000T00?oo00<00000003oo`00ool0 1000ool2000002<00?oo003Q003oo`0400000000ool00?oo000000T00?oo00<00000003oo`00ool0 1000ool00`0000000?oo003oo`0R003oo`00h@00ool0100000000?oo003oo`000009003oo`030000 0000ool00?oo00@00?oo00<00000003oo`00ool08P00ool00>800?oo0P000009003oo`8000001P00 ool00`0000000?oo003oo`0R003oo`003@00ool00`0000000?oo003oo`3U003oo`0300000000ool0 0?oo02800?oo000=003oo`0300000000ool00?oo0>D00?oo00<00000003oo`00ool08P00ool000h0 0?oo00<00000003oo`00ool0i000ool2000002<00?oo000?003oo`0300000000ool00?oo0><00?oo 00<00000003oo`00ool08P00ool001000?oo00<00000003oo`00ool0hP00ool00`0000000?oo003o o`0R003oo`004@00ool00`0000000?oo003oo`3Q003oo`0300000000ool00?oo02800?oo000B003o o`0300000000ool00?oo0>000?oo00<00000003oo`00ool08P00ool001<00?oo00<00000003oo`00 ool0g`00ool00`0000000?oo003oo`0R003oo`005000ool00`0000000?oo003oo`3N003oo`800000 8`00ool001D00?oo00<00000003oo`00ool0g@00ool00`0000000?oo003oo`0R003oo`005P00ool0 0`0000000?oo003oo`3L003oo`0300000000ool00?oo02800?oo000G003oo`0300000000ool00?oo 0=/00?oo00<00000003oo`00ool08P00ool001P00?oo00<00000003oo`00ool0fP00ool00`000000 0?oo003oo`0R003oo`006@00ool00`0000000?oo003oo`3I003oo`0300000000ool00?oo02800?oo 000J003oo`0300000000ool00?oo0=P00?oo00<00000003oo`00ool08P00ool001/00?oo00<00000 003oo`00ool0e`00ool2000002<00?oo000L003oo`0300000000ool00?oo0=H00?oo00<00000003o o`00ool08P00ool001d00?oo00<00000003oo`00ool0e@00ool00`0000000?oo003oo`0R003oo`00 7P00ool00`0000000?oo003oo`3D003oo`0300000000ool00?oo02800?oo000O003oo`0300000000 ool00?oo0=<00?oo00<00000003oo`00ool08P00ool002000?oo00<00000003oo`00ool0dP00ool0 0`0000000?oo003oo`0R003oo`008@00ool00`0000000?oo003oo`3A003oo`8000008`00ool00280 0?oo00<00000003oo`00ool0d000ool00`0000000?oo003oo`0R003oo`008`00ool00`0000000?oo 003oo`3?003oo`0300000000ool00?oo02800?oo000T003oo`0300000000ool00?oo0003oo`D0 000000D00?oo00000000ool00?oo00000008003oo`@000003P00ool500000005003oo`0000000?oo 003oo`0000002P00ool00`0000000?oo003oo`0H003oo`0300000000ool00?oo02800?oo000I003o o`0300000000ool00?oo00/00?oo00<00000003oo`00ool04@00ool00`0000000?oo003oo`05003o o`0300000000ool00?oo01X00?oo00@00000003oo`00ool000002@00ool2000001D00?oo00@00000 003oo`00ool000002000ool3000001D00?oo00@00000003oo`00ool000002000ool00`0000000?oo 0000000E003oo`0400000000ool00?oo000000/00?oo00<00000003oo`00ool05`00ool00`000000 0?oo003oo`0R003oo`006@00ool00`0000000?oo003oo`08003oo`0400000000ool00?oo000001@0 0?oo00<00000003oo`00ool01000ool00`0000000?oo003oo`0J003oo`0400000000ool00?oo0000 00P00?oo00@00000003oo`00ool000005000ool0100000000?oo003oo`000008003oo`0300000000 ool00?oo01D00?oo00@00000003oo`00ool000002000ool00`0000000?oo0000000E003oo`040000 0000ool00?oo000000P00?oo00@00000003oo`00ool000006@00ool2000002<00?oo000H003oo`80 00002`00ool2000001D00?oo00<00000003oo`00ool00`00ool2000001d00?oo0P00000:003oo`80 00005P00ool2000000X00?oo0`00000E003oo`8000002P00ool2000001H00?oo0P00000:003oo`80 00006P00ool00`0000000?oo003oo`0R003oo`00?@00ool00`0000000?oo003oo`2e003oo`030000 0000ool00?oo02800?oo000n003oo`0300000000ool00?oo0;@00?oo00<00000003oo`00ool08P00 ool003l00?oo00<00000003oo`00ool0/`00ool00`0000000?oo003oo`0<003oo`D000000P00ool3 000000`00?oo0010003oo`0300000000ool00?oo0;800?oo00<00000003oo`00ool03000ool01000 00000?oo003oo`000004003oo`0300000000ool00?oo00/00?oo0007003oooD000002000ool01000 00000?oo003oo`000004003oo`0300000000ool00?oo00/00?oo000<003oo`0300000000ool00?oo 00D00?oo00<00000003oo`00ool01P00ool00`0000000?oo003oo`06003oo`0300000000ool00?oo 00H00?oo00<00000003oo`00ool01P00ool00`0000000?oo003oo`06003oo`0300000000ool00?oo 00H00?oo00<00000003oo`00ool01P00ool00`0000000?oo003oo`06003oo`0300000000ool00?oo 00H00?oo00<00000003oo`00ool01P00ool00`0000000?oo003oo`06003oo`0300000000ool00?oo 00H00?oo00<00000003oo`00ool01P00ool00`0000000?oo003oo`06003oo`0300000000ool00?oo 00H00?oo00<00000003oo`00ool01P00ool00`0000000?oo003oo`06003oo`0300000000ool00?oo 00H00?oo00<00000003oo`00ool01P00ool00`0000000?oo003oo`06003oo`0300000000ool00?oo 00H00?oo00<00000003oo`00ool01P00ool00`0000000?oo003oo`06003oo`0300000000ool00?oo 00H00?oo00<00000003oo`00ool01P00ool00`0000000?oo003oo`0<003oo`0400000000ool00?oo 000000<00?oo0`00000<003oo`007@00ool00`0000000?oo003oo`0Q003oo`8000008P00ool00`00 00000?oo003oo`0Q003oo`0300000000ool00?oo02400?oo00<00000003oo`00ool08@00ool00`00 00000?oo003oo`0P003oo`8000003@00ool200000003003oo`000000000000<00?oo00<00000003o o`0000003000ool004<00?oo00<00000003oo`00ool0[P00ool2000001@00?oo0P0000000`00ool0 00000000000;003oo`00A000ool00`0000000?oo003oo`2/003oo`0300000000ool0000002@00?oo 0015003oo`0300000000ool00?oo0:/00?oo00<00000003oo`0000009000ool004H00?oo00<00000 003oo`00ool0Z@00ool0100000000?oo003oo`00000T003oo`00AP00ool00`0000000?oo003oo`2Y 003oo`0300000000ool00?oo008000008`00ool004L00?oo00<00000003oo`00ool0Z000ool01000 00000?oo003oo`00000T003oo`00B000ool00`0000000?oo003oo`2V003oo`0500000000ool00?oo 003oo`0000009000ool004T00?oo00<00000003oo`00ool0Y@00ool01@0000000?oo003oo`00ool0 000002@00?oo001:003oo`0300000000ool00?oo0:<00?oo00<00000003oo`00ool00P00ool00`00 00000?oo003oo`0R003oo`00B`00ool00`0000000?oo003oo`2R003oo`0300000000ool00?oo0080 0?oo00<00000003oo`00ool08P00ool004`00?oo00<00000003oo`00ool0X000ool00`0000000?oo 003oo`03003oo`8000008`00ool004`00?oo00<00000003oo`00ool0X000ool00`0000000?oo003o o`03003oo`0300000000ool00?oo02800?oo001=003oo`0300000000ool00?oo09h00?oo00<00000 003oo`00ool01000ool00`0000000?oo003oo`0R003oo`00CP00ool00`0000000?oo003oo`2M003o o`0300000000ool00?oo00@00?oo00<00000003oo`00ool08P00ool004l00?oo00<00000003oo`00 ool0V`00ool00`0000000?oo003oo`05003oo`0300000000ool00?oo02800?oo001@003oo`030000 0000ool00?oo09X00?oo00<00000003oo`00ool01@00ool00`0000000?oo003oo`0R003oo`00D@00 ool00`0000000?oo003oo`2I003oo`0300000000ool00?oo00D00?oo0P00000S003oo`00DP00ool0 0`0000000?oo003oo`2G003oo`0300000000ool00?oo00H00?oo00<00000003oo`00ool08P00ool0 05800?oo00<00000003oo`00ool0U`00ool00`0000000?oo003oo`06003oo`0300000000ool00?oo 02800?oo001C003oo`0300000000ool00?oo09D00?oo00<00000003oo`00ool01`00ool00`000000 0?oo003oo`0R003oo`00E000ool00`0000000?oo003oo`2D003oo`0300000000ool00?oo00L00?oo 00<00000003oo`00ool08P00ool005D00?oo00<00000003oo`00ool0TP00ool00`0000000?oo003o o`08003oo`0300000000ool00?oo02800?oo001F003oo`0300000000ool00?oo09400?oo00<00000 003oo`00ool02000ool2000002<00?oo001G003oo`0300000000ool00?oo08l00?oo00<00000003o o`00ool02@00ool00`0000000?oo003oo`0R003oo`00E`00ool00`0000000?oo003oo`2?003oo`03 00000000ool00?oo00T00?oo00<00000003oo`00ool08P00ool005P00?oo00<00000003oo`00ool0 S@00ool00`0000000?oo003oo`0:003oo`0300000000ool00?oo02800?oo001I003oo`0300000000 ool00?oo08000?oo0P000004003oo`8000001000ool2000000<00?oo0`000005003oo`0300000000 ool00?oo02800?oo001J003oo`0300000000ool00?oo07h00?oo00@00000003oo`00ool000000`00 ool2000000<00?oo0P0000000`00ool00000003oo`04003oo`0300000000ool00?oo00800?oo00<0 0000003oo`00ool08P00ool005/00?oo00<00000003oo`00ool0O@00ool0100000000?oo003oo`00 0008003oo`0400000000ool00?oo000000D00?oo00<00000003oo`00ool00P00ool2000002<00?oo 001L003oo`0300000000ool00?oo07`00?oo00@00000003oo`00ool000002000ool0100000000?oo 003oo`000002003oo`<000001@00ool00`0000000?oo003oo`0R003oo`00G@00ool00`0000000?oo 003oo`1k003oo`0400000000ool00?oo000000L00?oo0P000002003oo`0400000000ool00?oo0000 00L00?oo00<00000003oo`00ool08P00ool005d00?oo00<00000003oo`00ool0O000ool2000000P0 0?oo00@00000003oo`00000000000`00ool4000000@00?oo00<00000003oo`00ool08P00ool005h0 0?oo00<00000003oo`00ool0Q000ool00`0000000?oo003oo`0=003oo`0300000000ool00?oo0280 0?oo001O003oo`0300000000ool00?oo08800?oo00<00000003oo`00ool03P00ool00`0000000?oo 003oo`0R003oo`00H000ool00`0000000?oo003oo`21003oo`0300000000ool00?oo00h00?oo00<0 0000003oo`00ool08P00ool006400?oo00<00000003oo`00ool0O`00ool00`0000000?oo003oo`0? 003oo`8000008`00ool006800?oo00<00000003oo`00ool0O@00ool00`0000000?oo003oo`0@003o o`0300000000ool00?oo02800?oo001S003oo`0300000000ool00?oo07`00?oo00<00000003oo`00 ool04000ool00`0000000?oo003oo`0R003oo`00I000ool00`0000000?oo003oo`1j003oo`030000 0000ool00?oo01400?oo00<00000003oo`00ool08P00ool006D00?oo00<00000003oo`00ool0N000 ool00`0000000?oo003oo`0B003oo`0300000000ool00?oo02800?oo001V003oo`0300000000ool0 0?oo07L00?oo00<00000003oo`00ool04P00ool00`0000000?oo003oo`0R003oo`00I`00ool00`00 00000?oo003oo`1e003oo`0300000000ool00?oo01<00?oo0P00000S003oo`00J000ool00`000000 0?oo003oo`1c003oo`0300000000ool00?oo01@00?oo00<00000003oo`00ool08P00ool006T00?oo 00<00000003oo`00ool0LP00ool00`0000000?oo003oo`0D003oo`0300000000ool00?oo02800?oo 001Z003oo`0300000000ool00?oo07000?oo00<00000003oo`00ool05@00ool00`0000000?oo003o o`0R003oo`00J`00ool00`0000000?oo003oo`1_003oo`0300000000ool00?oo01D00?oo00<00000 003oo`00ool08P00ool006`00?oo00<00000003oo`00ool0K@00ool00`0000000?oo003oo`0F003o o`0300000000ool00?oo02800?oo001]003oo`0300000000ool00?oo06/00?oo00<00000003oo`00 ool05`00ool2000002<00?oo001^003oo`0300000000ool00?oo06X00?oo00<00000003oo`00ool0 5`00ool00`0000000?oo003oo`0R003oo`00K`00ool00`0000000?oo003oo`1X003oo`0300000000 ool00?oo01P00?oo00<00000003oo`00ool08P00ool007000?oo00<00000003oo`00ool0IP00ool0 0`0000000?oo003oo`0I003oo`0300000000ool00?oo02800?oo001a003oo`0300000000ool00?oo 06D00?oo00<00000003oo`00ool06@00ool00`0000000?oo003oo`0R003oo`00LP00ool00`000000 0?oo003oo`1S003oo`0300000000ool00?oo01X00?oo00<00000003oo`00ool08P00ool007<00?oo 00<00000003oo`00ool0H@00ool00`0000000?oo003oo`0K003oo`8000008`00ool007@00?oo00<0 0000003oo`00ool0H000ool00`0000000?oo003oo`0K003oo`0300000000ool00?oo02800?oo001e 003oo`0300000000ool00?oo05h00?oo00<00000003oo`00ool07000ool00`0000000?oo003oo`0R 003oo`00MP00ool00`0000000?oo003oo`1K003oo`8000007`00ool00`0000000?oo003oo`0R003o o`00M`00ool00`0000000?oo003oo`1I003oo`0300000000ool00?oo01l00?oo00<00000003oo`00 ool08P00ool007P00?oo00<00000003oo`00ool0E`00ool00`0000000?oo003oo`0=003oo`800000 1000ool2000000<00?oo0`000005003oo`0300000000ool00?oo02800?oo001i003oo`0300000000 ool00?oo05D00?oo00<00000003oo`00ool03@00ool0100000000?oo003oo`000003003oo`800000 1000ool00`0000000?oo003oo`04003oo`0300000000ool00?oo02800?oo001j003oo`0300000000 ool00?oo05<00?oo00<00000003oo`00ool03P00ool0100000000?oo003oo`000009003oo`030000 0000ool00?oo00@00?oo0P00000S003oo`00N`00ool00`0000000?oo003oo`1A003oo`0300000000 ool00?oo00l00?oo00@00000003oo`00ool000002@00ool00`0000000?oo003oo`04003oo`030000 0000ool00?oo02800?oo001l003oo`0300000000ool00?oo04l00?oo00<00000003oo`00ool04000 ool0100000000?oo003oo`000009003oo`0300000000ool00?oo00@00?oo00<00000003oo`00ool0 8P00ool007d00?oo00<00000003oo`00ool0C@00ool00`0000000?oo003oo`0B003oo`8000002@00 ool2000000H00?oo00<00000003oo`00ool08P00ool007h00?oo0P00001<003oo`0300000000ool0 0?oo02H00?oo00<00000003oo`00ool08P00ool008000?oo0P000019003oo`0300000000ool00?oo 02L00?oo00<00000003oo`00ool08P00ool008800?oo00<00000003oo`00ool0A@00ool00`000000 0?oo003oo`0X003oo`8000008`00ool008<00?oo0P000014003oo`0300000000ool00?oo02T00?oo 00<00000003oo`00ool08P00ool008D00?oo0P000011003oo`0300000000ool00?oo02X00?oo00<0 0000003oo`00ool08P00ool008L00?oo0P00000n003oo`0300000000ool00?oo02/00?oo00<00000 003oo`00ool08P00ool008T00?oo00<00000003oo`00ool0>P00ool00`0000000?oo003oo`0/003o o`0300000000ool00?oo02800?oo002:003oo`800000>@00ool00`0000000?oo003oo`0]003oo`03 00000000ool00?oo02800?oo002<003oo`800000=P00ool00`0000000?oo003oo`0^003oo`800000 8`00ool008h00?oo00<00000003oo`00ool0`00ool00`0000000?oo 003oo`0=003oo`0300000000ool00?oo00L00?oo00<00000003oo`00ool08`00ool200000003003o o`000000000002<00?oo00@00000003oo`00ool00?oo10000002003oo`0700000000ool00000003o o`0000000?oo00000003003oo`0300000000ool0000000<00?oo00<00000003oo`0000000`00000D 003oo`00O000ool3000000h00?oo0P000008003oo`800000C@00ool00`0000000?oo003oo`04003o o`0:00000000ool00?oo00000000ool00000003oo`0000000?oo000000<00?oo00<00000003oo`00 00000`00ool00`0000000?oo003oo`02003oo`0300000000ool00?oo01800?oo003V003oo`030000 0000ool00?oo00800?oo0`0000000`00ool000000000000400000003003oo`000000000000800000 0`00ool4000000800?oo0`00000D003oo`00iP00ool00`0000000?oo003oo`0=003oo`0300000000 ool00?oo00L00?oo00<00000003oo`00ool05`00ool00>@00?oo0`00000>003oo`8000002000ool2 000001T00?oo003o003ooa/00?oo003o003ooa/00?oo0000\ \>"], ImageRangeCache->{{{0, 281}, {173, 0}} -> {-1.36689, -0.107224, 0.00557084, 0.00160536}}], Cell[GraphicsData["PostScript", "\<\ %! %%Creator: Mathematica %%AspectRatio: .61803 MathPictureStart /Mabs { Mgmatrix idtransform Mtmatrix dtransform } bind def /Mabsadd { Mabs 3 -1 roll add 3 1 roll add exch } bind def %% Graphics %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10 scalefont setfont % Scaling calculations 0.97619 0.737473 0.210625 0.392694 [ [.09122 .19813 -12 -9 ] [.09122 .19813 12 0 ] [.23872 .19813 -6 -9 ] [.23872 .19813 6 0 ] [.38621 .19813 -12 -9 ] [.38621 .19813 12 0 ] [.53371 .19813 -12 -9 ] [.53371 .19813 12 0 ] [.6812 .19813 -12 -9 ] [.6812 .19813 12 0 ] [.8287 .19813 -12 -9 ] [.8287 .19813 12 0 ] [1.025 .21063 0 -6 ] [1.025 .21063 16 6 ] [.96369 .05355 -24 -4.5 ] [.96369 .05355 0 4.5 ] [.96369 .13209 -24 -4.5 ] [.96369 .13209 0 4.5 ] [.96369 .28916 -18 -4.5 ] [.96369 .28916 0 4.5 ] [.96369 .3677 -18 -4.5 ] [.96369 .3677 0 4.5 ] [.96369 .44624 -18 -4.5 ] [.96369 .44624 0 4.5 ] [.96369 .52478 -18 -4.5 ] [.96369 .52478 0 4.5 ] [.96369 .60332 -6 -4.5 ] [.96369 .60332 0 4.5 ] [.97619 .64303 -17 0 ] [.97619 .64303 17 12 ] [.5 .63053 -65 0 ] [.5 .63053 65 12 ] [ 0 0 0 0 ] [ 1 .61803 0 0 ] ] MathScale % Start of Graphics 1 setlinecap 1 setlinejoin newpath 0 g .25 Mabswid [ ] 0 setdash .09122 .21063 m .09122 .21688 L s [(-1.2)] .09122 .19813 0 1 Mshowa .23872 .21063 m .23872 .21688 L s [(-1)] .23872 .19813 0 1 Mshowa .38621 .21063 m .38621 .21688 L s [(-0.8)] .38621 .19813 0 1 Mshowa .53371 .21063 m .53371 .21688 L s [(-0.6)] .53371 .19813 0 1 Mshowa .6812 .21063 m .6812 .21688 L s [(-0.4)] .6812 .19813 0 1 Mshowa .8287 .21063 m .8287 .21688 L s [(-0.2)] .8287 .19813 0 1 Mshowa .125 Mabswid .1281 .21063 m .1281 .21438 L s .16497 .21063 m .16497 .21438 L s .20184 .21063 m .20184 .21438 L s .27559 .21063 m .27559 .21438 L s .31246 .21063 m .31246 .21438 L s .34934 .21063 m .34934 .21438 L s .42309 .21063 m .42309 .21438 L s .45996 .21063 m .45996 .21438 L s .49683 .21063 m .49683 .21438 L s .57058 .21063 m .57058 .21438 L s .60745 .21063 m .60745 .21438 L s .64433 .21063 m .64433 .21438 L s .71807 .21063 m .71807 .21438 L s .75495 .21063 m .75495 .21438 L s .79182 .21063 m .79182 .21438 L s .86557 .21063 m .86557 .21438 L s .90244 .21063 m .90244 .21438 L s .93932 .21063 m .93932 .21438 L s .05435 .21063 m .05435 .21438 L s .01748 .21063 m .01748 .21438 L s .25 Mabswid 0 .21063 m 1 .21063 L s gsave 1.025 .21063 -61 -10 Mabsadd m 1 1 Mabs scale currentpoint translate 0 20 translate 1 -1 scale gsave 0.000000 0.000000 0.000000 setrgbcolor 1.000000 setlinewidth gsave newpath 61.000000 16.000000 moveto 558.000000 16.000000 lineto 558.000000 4.000000 lineto 61.000000 4.000000 lineto 61.000000 16.000000 lineto closepath clip newpath 63.000000 13.000000 moveto %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10.000000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000000 0.000000 0.000000 setrgbcolor 0.000000 0.000000 rmoveto 63.000000 13.000000 moveto %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10.000000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000000 0.000000 0.000000 setrgbcolor (uY) show 75.000000 13.000000 moveto %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10.000000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000000 0.000000 0.000000 setrgbcolor 0.000000 0.000000 rmoveto 1.000000 setlinewidth grestore grestore %%DocumentNeededResources: font Courier %%DocumentSuppliedResources: %%DocumentNeededFonts: Courier %%DocumentSuppliedFonts: %%DocumentFonts: font Courier grestore .97619 .05355 m .98244 .05355 L s [(-0.4)] .96369 .05355 1 0 Mshowa .97619 .13209 m .98244 .13209 L s [(-0.2)] .96369 .13209 1 0 Mshowa .97619 .28916 m .98244 .28916 L s [(0.2)] .96369 .28916 1 0 Mshowa .97619 .3677 m .98244 .3677 L s [(0.4)] .96369 .3677 1 0 Mshowa .97619 .44624 m .98244 .44624 L s [(0.6)] .96369 .44624 1 0 Mshowa .97619 .52478 m .98244 .52478 L s [(0.8)] .96369 .52478 1 0 Mshowa .97619 .60332 m .98244 .60332 L s [(1)] .96369 .60332 1 0 Mshowa .125 Mabswid .97619 .07318 m .97994 .07318 L s .97619 .09282 m .97994 .09282 L s .97619 .11245 m .97994 .11245 L s .97619 .15172 m .97994 .15172 L s .97619 .17136 m .97994 .17136 L s .97619 .19099 m .97994 .19099 L s .97619 .23026 m .97994 .23026 L s .97619 .24989 m .97994 .24989 L s .97619 .26953 m .97994 .26953 L s .97619 .3088 m .97994 .3088 L s .97619 .32843 m .97994 .32843 L s .97619 .34807 m .97994 .34807 L s .97619 .38734 m .97994 .38734 L s .97619 .40697 m .97994 .40697 L s .97619 .42661 m .97994 .42661 L s .97619 .46588 m .97994 .46588 L s .97619 .48551 m .97994 .48551 L s .97619 .50515 m .97994 .50515 L s .97619 .54441 m .97994 .54441 L s .97619 .56405 m .97994 .56405 L s .97619 .58368 m .97994 .58368 L s .97619 .03391 m .97994 .03391 L s .97619 .01428 m .97994 .01428 L s .25 Mabswid .97619 0 m .97619 .61803 L s gsave .97619 .64303 -78 -4 Mabsadd m 1 1 Mabs scale currentpoint translate 0 20 translate 1 -1 scale gsave 0.000000 0.000000 0.000000 setrgbcolor 1.000000 setlinewidth gsave newpath 61.000000 16.000000 moveto 558.000000 16.000000 lineto 558.000000 4.000000 lineto 61.000000 4.000000 lineto 61.000000 16.000000 lineto closepath clip newpath 63.000000 13.000000 moveto %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10.000000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000000 0.000000 0.000000 setrgbcolor 0.000000 0.000000 rmoveto 63.000000 13.000000 moveto %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10.000000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000000 0.000000 0.000000 setrgbcolor (kappa) show 93.000000 13.000000 moveto %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10.000000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000000 0.000000 0.000000 setrgbcolor 0.000000 0.000000 rmoveto 1.000000 setlinewidth grestore grestore %%DocumentNeededResources: font Courier %%DocumentSuppliedResources: %%DocumentNeededFonts: Courier %%DocumentSuppliedFonts: %%DocumentFonts: font Courier grestore gsave .5 .63053 -126 -4 Mabsadd m 1 1 Mabs scale currentpoint translate 0 20 translate 1 -1 scale gsave 0.000000 0.000000 0.000000 setrgbcolor 1.000000 setlinewidth gsave newpath 61.000000 16.000000 moveto 558.000000 16.000000 lineto 558.000000 4.000000 lineto 61.000000 4.000000 lineto 61.000000 16.000000 lineto closepath clip newpath 63.000000 13.000000 moveto %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10.000000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000000 0.000000 0.000000 setrgbcolor 0.000000 0.000000 rmoveto 63.000000 13.000000 moveto %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10.000000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000000 0.000000 0.000000 setrgbcolor (LP) show 81.000000 13.000000 moveto (sensor) show 123.000000 13.000000 moveto (kappa) show 159.000000 13.000000 moveto (vs) show 177.000000 13.000000 moveto (uY) show 189.000000 13.000000 moveto %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10.000000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000000 0.000000 0.000000 setrgbcolor 0.000000 0.000000 rmoveto 1.000000 setlinewidth grestore grestore %%DocumentNeededResources: font Courier %%DocumentSuppliedResources: %%DocumentNeededFonts: Courier %%DocumentSuppliedFonts: %%DocumentFonts: font Courier grestore .5 Mabswid .97619 .60332 m .91338 .60332 L .84475 .50726 L .76925 .41206 L .68568 .31912 L .59289 .23065 L .48975 .15013 L .37738 .08253 L .2588 .0351 L .13877 .01472 L .02381 .0251 L s 0 0 m 1 0 L 1 .61803 L 0 .61803 L closepath clip newpath % End of Graphics MathPictureEnd \ \>"], "Graphics", Evaluatable->False, AspectRatioFixed->True, ImageSize->{282, 174}, ImageMargins->{{34, 0}, {0, 0}}, ImageRegion->{{0, 1}, {0, 1}}, Background->RGBColor[0, 1, 1], ImageCache->GraphicsData["Bitmap", "\<\ CF5dJ6E]HGAYHf4PAg9QL6QYHg@00ool01000 00000?oo003oo`00000:003oo`0300000000ool00?oo00<00?oo0`00000O003oo`00Z@00ool00`00 00000?oo003oo`0h003oo`0400000000ool00?oo000000/00?oo00<00000003oo`00ool00P00ool0 0`0000000?oo003oo`0O003oo`00ZP00ool00`0000000?oo003oo`0g003oo`0400000000ool00?oo 000000P00?oo00@00000003oo`00ool000001000ool00`0000000?oo003oo`0O003oo`00Z`00ool0 0`0000000?oo003oo`0g003oo`8000002P00ool2000000D00?oo00<00000003oo`00ool07`00ool0 0:`00?oo00<00000003oo`00ool0B@00ool00`0000000?oo003oo`0O003oo`00[@00ool00`000000 0?oo003oo`18003oo`8000008000ool00:h00?oo00<00000003oo`00ool0A`00ool00`0000000?oo 003oo`0O003oo`00[`00ool00`0000000?oo003oo`16003oo`0300000000ool00?oo01l00?oo002_ 003oo`0300000000ool00?oo04H00?oo00<00000003oo`00ool07`00ool00;000?oo00<00000003o o`00ool0A@00ool00`0000000?oo003oo`0O003oo`00/@00ool00`0000000?oo003oo`14003oo`80 00008000ool00;800?oo00<00000003oo`00ool0@`00ool00`0000000?oo003oo`0O003oo`00/`00 ool00`0000000?oo003oo`12003oo`0300000000ool00?oo01l00?oo002d003oo`0300000000ool0 0?oo04400?oo00<00000003oo`00ool07`00ool00;D00?oo00<00000003oo`00ool0@000ool00`00 00000?oo003oo`0O003oo`00]P00ool00`0000000?oo003oo`0o003oo`8000008000ool00;H00?oo 00<00000003oo`00ool0?`00ool00`0000000?oo003oo`0O003oo`00]`00ool00`0000000?oo003o o`0n003oo`0300000000ool00?oo01l00?oo002h003oo`0300000000ool00?oo02X00?oo0P000004 003oo`8000001000ool3000000@00?oo00<00000003oo`00ool07`00ool00;T00?oo00<00000003o o`00ool0:000ool0100000000?oo003oo`000003003oo`8000001@00ool00`0000000?oo003oo`03 003oo`0300000000ool00?oo01l00?oo002j003oo`0300000000ool00?oo02L00?oo00@00000003o o`00ool000002000ool4000000@00?oo0`00000O003oo`00^`00ool00`0000000?oo003oo`0V003o o`0400000000ool00?oo000000P00?oo00<00000003oo`0000001@00ool00`0000000?oo003oo`0O 003oo`00_000ool00`0000000?oo003oo`0U003oo`0400000000ool00?oo000000P00?oo00<00000 003oo`0000001@00ool00`0000000?oo003oo`0O003oo`00_@00ool00`0000000?oo003oo`0U003o o`8000002P00ool2000000D00?oo00<00000003oo`00ool07`00ool00;d00?oo00<00000003oo`00 ool0>000ool2000002000?oo002n003oo`0300000000ool00?oo03L00?oo00<00000003oo`00ool0 7`00ool00;l00?oo00<00000003oo`00ool0=P00ool00`0000000?oo003oo`0O003oo`00`000ool0 0`0000000?oo003oo`0e003oo`0300000000ool00?oo01l00?oo0031003oo`0300000000ool00?oo 03@00?oo00<00000003oo`00ool07`00ool00<800?oo00<00000003oo`00ool0<`00ool200000200 0?oo0033003oo`0300000000ool00?oo03800?oo00<00000003oo`00ool07`00ool00<@00?oo00<0 0000003oo`00ool0<@00ool00`0000000?oo003oo`0O003oo`00a000ool00`0000000?oo003oo`0a 003oo`0300000000ool00?oo01l00?oo0035003oo`0300000000ool00?oo03000?oo00<00000003o o`00ool07`00ool00400?oo00<00000003oo`00ool05000ool2000002000?oo003Q003oo`0300000000ool0 0?oo01@00?oo00<00000003oo`00ool07`00ool00>800?oo00<00000003oo`00ool04`00ool00`00 00000?oo003oo`0O003oo`00h`00ool00`0000000?oo003oo`0B003oo`0300000000ool00?oo01l0 0?oo003S003oo`0300000000ool00?oo01800?oo00<00000003oo`00ool07`00ool00>@00?oo00<0 0000003oo`00ool04@00ool2000002000?oo003U003oo`0300000000ool00?oo01000?oo00<00000 003oo`00ool07`00ool00>H00?oo00<00000003oo`00ool03`00ool00`0000000?oo003oo`0O003o o`00iP00ool00`0000000?oo003oo`07003oo`<000001@00ool00`0000000?oo003oo`0O003oo`00 i`00ool00`0000000?oo003oo`07003oo`0300000000ool00?oo00@00?oo00<00000003oo`00ool0 7`00ool00>P00?oo4`00000O003oo`00l@00ool00`0000000?oo003oo`04003oo`0300000000ool0 0?oo01l00?oo003a003oo`0300000000ool00?oo00@00?oo00<00000003oo`00ool07`00ool00?00 0?oo0P000006003oo`0300000000ool00?oo01l00?oo003o003ooa/00?oo003o003ooa/00?oo003o 003ooa/00?oo003o003ooa/00?oo002:003oo`@000000P00ool4000008H00?oo002;003oo`030000 0000ool00?oo00<00?oo00<00000003oo`00ool0QP00ool004800?oo1@0000001000ool000000000 0000000:003oo`@000000`00ool400000004003oo`0000000000003oo`80000000<00?oo00000000 00000P000003003oo`<000000P00ool4000000P00?oo00@00000003oo`00000000000P00ool50000 0003003oo`0000000000008000000P00ool4000000800?oo1@000008003oo`8000000`00ool40000 00P00?oo1@000002003oo`<00000F`00ool004<00?oo00@00000003oo`00ool000000P00ool00`00 00000?oo003oo`0<003oo`0400000000ool00?oo000000D00?oo00@00000003oo`00ool000001@00 ool0100000000?oo003oo`000003003oo`0400000000ool00?oo000000X00?oo00<00000003oo`00 00000`00ool0100000000?oo003oo`000002003oo`0700000000ool00?oo003oo`0000000?oo0000 0003003oo`0300000000ool0000000800?oo00<00000003oo`00ool01`00ool2000000H00?oo00<0 0000003oo`00ool01P00ool0100000000?oo003oo`000004003oo`0300000000ool00?oo03D00?oo 10000002003oo`@000006`00ool004<00?oo00@00000003oo`00ool000000P00ool3000000T00?oo 10000002003oo`D0000000D00?oo00000000ool00?oo00000002003oo`@000000P00ool01@000000 0?oo003oo`00ool0000000800?oo00<00000003oo`00ool02000ool3000000<00?oo10000002003o o`0700000000ool00?oo003oo`0000000?oo00000003003oo`0300000000ool0000000<000002000 ool0100000000?oo000000000002003oo`@000002000ool0100000000?oo003oo`000004003oo`03 00000000ool00?oo03H00?oo00<00000003oo`00ool00`00ool00`0000000?oo003oo`0K003oo`00 @`00ool00`0000000?oo003oo`03003oo`0400000000ool00?oo000000P00?oo00<00000003oo`00 ool00`00ool01`0000000?oo003oo`00ool00000003oo`0000000P00ool0100000000?oo003oo`00 0005003oo`0500000000ool00?oo003oo`0000000P00ool2000000T00?oo00<00000003oo`000000 1P00ool0100000000?oo003oo`000003003oo`0300000000ool0000000<00?oo00<00000003oo`00 ool00P00ool00`0000000?oo003oo`06003oo`0400000000ool00?oo000000800?oo00<00000003o o`00ool02@00ool0100000000?oo003oo`000003003oo`<00000:`00ool0100000000?oo00000000 0002003oo`D0000000<00?oo0000000000000P000002003oo`@000000P00ool5000001<00?oo0013 003oo`0300000000ool00?oo00<00?oo00@00000003oo`00ool000002000ool4000000<00?oo0`00 00000`00ool0000000000003000000800?oo10000003003oo`<000000P00ool200000003003oo`00 0000000000L00?oo00@00000003oo`00000000000`00ool300000003003oo`000000000000<00000 00<00?oo0000000000000`000003003oo`<000001`00ool2000000800?oo0P0000000`00ool00000 00000002000000L00?oo0P0000000`00ool0000000000003003oo`0300000000ool0000002/00?oo 00<00000003oo`0000000`00ool0100000000?oo003oo`000002003oo`0700000000ool00?oo003o o`0000000?oo00000003003oo`0300000000ool0000000800?oo00<00000003oo`00ool04P00ool0 04800?oo0`000003003oo`@00000<`00ool00`0000000?oo003oo`0i003oo`80000000<00?oo0000 00000000:P00ool3000000<00?oo10000002003oo`0700000000ool00?oo003oo`0000000?oo0000 0003003oo`0300000000ool0000000<000005000ool007h00?oo0P00001Z003oo`0300000000ool0 000000H00?oo00@00000003oo`00ool000000`00ool00`0000000?oo00000003003oo`0300000000 ool00?oo00800?oo00<00000003oo`00ool04P00ool00>X00?oo00@00000003oo`00000000000`00 ool300000003003oo`000000000000<0000000<00?oo0000000000000`000003003oo`<000005000 ool00>X00?oo00<00000003oo`00ool0;@00ool00>T00?oo0P00000_003oo`00o`00oolK003oo`00 o`00oolK003oo`00\ \>"], ImageRangeCache->{{{0, 281}, {173, 0}} -> {-1.34635, -0.536367, 0.00541581, 0.0101708}}] }, Open ]] }, Open ]] }, FrontEndVersion->"4.2 for Macintosh", ScreenRectangle->{{0, 1920}, {0, 1180}}, AutoGeneratedPackage->None, WindowToolbars->{}, CellGrouping->Manual, WindowSize->{1026, 1093}, WindowMargins->{{12, Automatic}, {Automatic, 13}}, PrivateNotebookOptions->{"ColorPalette"->{RGBColor, -1}}, ShowCellLabel->True, ShowCellTags->False, RenderingOptions->{"ObjectDithering"->True, "RasterDithering"->False}, MacintoshSystemPageSetup->"\<\ 00<0001804P000000]P2:?oQon82n@960dL5:0?l0080001804P000000]P2:001 0000I00000400`<300000BL?00400@00000000000000060801T1T00000000000 00000000000000000000000000000000\>" ] (******************************************************************* Cached data follows. If you edit this Notebook file directly, not using Mathematica, you must remove the line containing CacheID at the top of the file. The cache data will then be recreated when you save this file from within Mathematica. *******************************************************************) (*CellTagsOutline CellTagsIndex->{} *) (*CellTagsIndex CellTagsIndex->{} *) (*NotebookFileOutline Notebook[{ Cell[CellGroupData[{ Cell[1776, 53, 344, 9, 72, "Section", Evaluatable->False], Cell[2123, 64, 1062, 18, 336, "Text", Evaluatable->False], Cell[3188, 84, 4543, 116, 280, "Text", Evaluatable->False], Cell[CellGroupData[{ Cell[7756, 204, 1800, 37, 433, "Input", InitializationCell->True], Cell[9559, 243, 194, 4, 22, "Print"], Cell[9756, 249, 164, 3, 43, "Print"] }, Open ]], Cell[9935, 255, 1147, 31, 64, "Text", Evaluatable->False], Cell[CellGroupData[{ Cell[11107, 290, 1265, 29, 193, "Input", InitializationCell->True], Cell[12375, 321, 46, 1, 22, "Print"], Cell[12424, 324, 47, 1, 22, "Print"] }, Open ]], Cell[12486, 328, 2489, 65, 130, "Text", Evaluatable->False], Cell[CellGroupData[{ Cell[15000, 397, 5238, 119, 583, "Input", InitializationCell->True], Cell[20241, 518, 232, 5, 22, "Print"], Cell[20476, 525, 141, 3, 22, "Print"], Cell[20620, 530, 189, 4, 22, "Print"], Cell[20812, 536, 177, 3, 22, "Print"] }, Open ]], Cell[21004, 542, 7035, 141, 808, "Text", Evaluatable->False], Cell[28042, 685, 3036, 59, 688, "Input", InitializationCell->True], Cell[31081, 746, 1161, 32, 64, "Text", Evaluatable->False], Cell[32245, 780, 606, 16, 163, "Input", InitializationCell->True], Cell[32854, 798, 1660, 46, 64, "Text", Evaluatable->False], Cell[34517, 846, 2378, 49, 643, "Input", InitializationCell->True], Cell[36898, 897, 1160, 32, 64, "Text", Evaluatable->False], Cell[38061, 931, 606, 16, 163, "Input", InitializationCell->True], Cell[38670, 949, 1670, 46, 58, "Input"], Cell[40343, 997, 2970, 61, 823, "Input", InitializationCell->True], Cell[43316, 1060, 1174, 32, 64, "Text", Evaluatable->False], Cell[44493, 1094, 655, 17, 178, "Input", InitializationCell->True], Cell[45151, 1113, 1686, 46, 58, "Input"], Cell[46840, 1161, 1641, 37, 478, "Input", InitializationCell->True], Cell[48484, 1200, 316, 8, 46, "Text", Evaluatable->False], Cell[CellGroupData[{ Cell[48825, 1212, 1214, 31, 388, "Input", InitializationCell->True], Cell[50042, 1245, 794, 16, 41, "Print"] }, Open ]], Cell[50851, 1264, 614, 12, 126, "Text", Evaluatable->False], Cell[51468, 1278, 3987, 87, 598, "Input", InitializationCell->True], Cell[55458, 1367, 4195, 106, 340, "Text", Evaluatable->False], Cell[59656, 1475, 12403, 284, 838, "Input", InitializationCell->True], Cell[72062, 1761, 1185, 31, 64, "Text", Evaluatable->False], Cell[CellGroupData[{ Cell[73272, 1796, 2069, 33, 883, "Input"], Cell[75344, 1831, 6144, 85, 1181, "Print", Evaluatable->False], Cell[CellGroupData[{ Cell[81513, 1920, 5128, 493, 198, 4805, 484, "GraphicsData", "PostScript", \ "Graphics", Evaluatable->False, ImageCacheValid->False], Cell[86644, 2415, 5542, 515, 198, 5219, 506, "GraphicsData", "PostScript", \ "Graphics", Evaluatable->False, ImageCacheValid->False] }, Open ]] }, Open ]], Cell[92213, 2934, 1367, 34, 88, "Input"], Cell[CellGroupData[{ Cell[93605, 2972, 1938, 63, 898, "Input"], Cell[95546, 3037, 54, 1, 22, "Print"], Cell[95603, 3040, 251, 5, 22, "Print"], Cell[95857, 3047, 479, 10, 22, "Print"], Cell[96339, 3059, 360, 8, 22, "Print"], Cell[96702, 3069, 357, 7, 22, "Print"], Cell[97062, 3078, 417, 9, 22, "Print"], Cell[97482, 3089, 188, 4, 24, "Print"], Cell[97673, 3095, 186, 4, 22, "Print"], Cell[97862, 3101, 476, 10, 22, "Print"], Cell[98341, 3113, 909, 20, 22, "Print"], Cell[99253, 3135, 54, 1, 22, "Print"], Cell[99310, 3138, 251, 5, 22, "Print"], Cell[99564, 3145, 493, 11, 22, "Print"], Cell[100060, 3158, 360, 8, 22, "Print"], Cell[100423, 3168, 356, 7, 22, "Print"], Cell[100782, 3177, 429, 9, 22, "Print"], Cell[101214, 3188, 175, 3, 24, "Print"], Cell[101392, 3193, 189, 4, 24, "Print"], Cell[101584, 3199, 474, 10, 24, "Print"], Cell[102061, 3211, 904, 20, 24, "Print"], Cell[102968, 3233, 222, 5, 22, "Print"], Cell[103193, 3240, 54, 1, 22, "Print"], Cell[103250, 3243, 251, 5, 22, "Print"], Cell[103504, 3250, 477, 10, 22, "Print"], Cell[103984, 3262, 449, 9, 22, "Print"], Cell[104436, 3273, 357, 7, 22, "Print"], Cell[104796, 3282, 419, 9, 22, "Print"], Cell[105218, 3293, 186, 4, 22, "Print"], Cell[105407, 3299, 186, 4, 22, "Print"], Cell[105596, 3305, 462, 9, 22, "Print"], Cell[106061, 3316, 906, 20, 22, "Print"], Cell[106970, 3338, 54, 1, 22, "Print"], Cell[107027, 3341, 251, 5, 22, "Print"], Cell[107281, 3348, 481, 10, 22, "Print"], Cell[107765, 3360, 449, 9, 22, "Print"], Cell[108217, 3371, 357, 7, 22, "Print"], Cell[108577, 3380, 429, 9, 22, "Print"], Cell[109009, 3391, 188, 4, 24, "Print"], Cell[109200, 3397, 190, 4, 22, "Print"], Cell[109393, 3403, 462, 9, 24, "Print"], Cell[109858, 3414, 906, 20, 24, "Print"], Cell[110767, 3436, 222, 5, 22, "Print"], Cell[110992, 3443, 54, 1, 22, "Print"], Cell[111049, 3446, 251, 5, 22, "Print"], Cell[111303, 3453, 476, 10, 22, "Print"], Cell[111782, 3465, 446, 9, 22, "Print"], Cell[112231, 3476, 357, 7, 22, "Print"], Cell[112591, 3485, 419, 9, 22, "Print"], Cell[113013, 3496, 177, 3, 22, "Print"], Cell[113193, 3501, 186, 4, 22, "Print"], Cell[113382, 3507, 461, 9, 22, "Print"], Cell[113846, 3518, 905, 20, 22, "Print"], Cell[114754, 3540, 54, 1, 22, "Print"], Cell[114811, 3543, 251, 5, 22, "Print"], Cell[115065, 3550, 480, 10, 22, "Print"], Cell[115548, 3562, 446, 9, 22, "Print"], Cell[115997, 3573, 357, 7, 22, "Print"], Cell[116357, 3582, 427, 9, 22, "Print"], Cell[116787, 3593, 188, 4, 24, "Print"], Cell[116978, 3599, 188, 4, 22, "Print"], Cell[117169, 3605, 459, 9, 24, "Print"], Cell[117631, 3616, 901, 20, 24, "Print"], Cell[118535, 3638, 222, 5, 22, "Print"], Cell[118760, 3645, 54, 1, 22, "Print"], Cell[118817, 3648, 251, 5, 22, "Print"], Cell[119071, 3655, 476, 10, 22, "Print"], Cell[119550, 3667, 445, 9, 22, "Print"], Cell[119998, 3678, 357, 7, 22, "Print"], Cell[120358, 3687, 417, 9, 22, "Print"], Cell[120778, 3698, 177, 3, 22, "Print"], Cell[120958, 3703, 186, 4, 22, "Print"], Cell[121147, 3709, 458, 9, 22, "Print"], Cell[121608, 3720, 899, 20, 22, "Print"], Cell[122510, 3742, 54, 1, 22, "Print"], Cell[122567, 3745, 251, 5, 22, "Print"], Cell[122821, 3752, 477, 10, 22, "Print"], Cell[123301, 3764, 445, 9, 22, "Print"], Cell[123749, 3775, 357, 7, 22, "Print"], Cell[124109, 3784, 423, 9, 22, "Print"], Cell[124535, 3795, 175, 3, 24, "Print"], Cell[124713, 3800, 188, 4, 22, "Print"], Cell[124904, 3806, 460, 9, 24, "Print"], Cell[125367, 3817, 903, 20, 24, "Print"], Cell[126273, 3839, 222, 5, 22, "Print"], Cell[126498, 3846, 54, 1, 22, "Print"], Cell[126555, 3849, 251, 5, 22, "Print"], Cell[126809, 3856, 474, 10, 22, "Print"], Cell[127286, 3868, 444, 9, 22, "Print"], Cell[127733, 3879, 351, 7, 22, "Print"], Cell[128087, 3888, 412, 9, 22, "Print"], Cell[128502, 3899, 175, 3, 22, "Print"], Cell[128680, 3904, 186, 4, 22, "Print"], Cell[128869, 3910, 459, 9, 22, "Print"], Cell[129331, 3921, 900, 20, 22, "Print"], Cell[130234, 3943, 54, 1, 22, "Print"], Cell[130291, 3946, 251, 5, 22, "Print"], Cell[130545, 3953, 478, 10, 22, "Print"], Cell[131026, 3965, 444, 9, 22, "Print"], Cell[131473, 3976, 351, 7, 22, "Print"], Cell[131827, 3985, 415, 9, 22, "Print"], Cell[132245, 3996, 133, 3, 22, "Print"], Cell[132381, 4001, 188, 4, 22, "Print"], Cell[132572, 4007, 458, 9, 24, "Print"], Cell[133033, 4018, 901, 20, 24, "Print"], Cell[133937, 4040, 222, 5, 22, "Print"], Cell[134162, 4047, 54, 1, 22, "Print"], Cell[134219, 4050, 251, 5, 22, "Print"], Cell[134473, 4057, 476, 10, 22, "Print"], Cell[134952, 4069, 445, 9, 22, "Print"], Cell[135400, 4080, 351, 7, 22, "Print"], Cell[135754, 4089, 413, 9, 22, "Print"], Cell[136170, 4100, 175, 3, 22, "Print"], Cell[136348, 4105, 177, 3, 22, "Print"], Cell[136528, 4110, 459, 9, 22, "Print"], Cell[136990, 4121, 900, 20, 22, "Print"], Cell[137893, 4143, 54, 1, 22, "Print"], Cell[137950, 4146, 251, 5, 22, "Print"], Cell[138204, 4153, 478, 10, 22, "Print"], Cell[138685, 4165, 445, 9, 22, "Print"], Cell[139133, 4176, 351, 7, 22, "Print"], Cell[139487, 4185, 417, 9, 22, "Print"], Cell[139907, 4196, 175, 3, 24, "Print"], Cell[140085, 4201, 188, 4, 22, "Print"], Cell[140276, 4207, 459, 9, 24, "Print"], Cell[140738, 4218, 902, 20, 24, "Print"], Cell[141643, 4240, 54, 1, 22, "Print"], Cell[141700, 4243, 251, 5, 22, "Print"], Cell[141954, 4250, 478, 10, 24, "Print"], Cell[142435, 4262, 445, 9, 22, "Print"], Cell[142883, 4273, 350, 7, 22, "Print"], Cell[143236, 4282, 420, 9, 24, "Print"], Cell[143659, 4293, 188, 4, 24, "Print"], Cell[143850, 4299, 189, 4, 24, "Print"], Cell[144042, 4305, 461, 9, 24, "Print"], Cell[144506, 4316, 905, 20, 24, "Print"], Cell[145414, 4338, 222, 5, 22, "Print"], Cell[145639, 4345, 54, 1, 22, "Print"], Cell[145696, 4348, 251, 5, 22, "Print"], Cell[145950, 4355, 477, 10, 22, "Print"], Cell[146430, 4367, 445, 9, 22, "Print"], Cell[146878, 4378, 351, 7, 22, "Print"], Cell[147232, 4387, 413, 9, 22, "Print"], Cell[147648, 4398, 175, 3, 22, "Print"], Cell[147826, 4403, 177, 3, 22, "Print"], Cell[148006, 4408, 460, 9, 22, "Print"], Cell[148469, 4419, 901, 20, 22, "Print"], Cell[149373, 4441, 54, 1, 22, "Print"], Cell[149430, 4444, 251, 5, 22, "Print"], Cell[149684, 4451, 468, 9, 22, "Print"], Cell[150155, 4462, 445, 9, 22, "Print"], Cell[150603, 4473, 351, 7, 22, "Print"], Cell[150957, 4482, 417, 9, 22, "Print"], Cell[151377, 4493, 174, 3, 24, "Print"], Cell[151554, 4498, 188, 4, 22, "Print"], Cell[151745, 4504, 462, 9, 24, "Print"], Cell[152210, 4515, 905, 20, 24, "Print"], Cell[153118, 4537, 54, 1, 22, "Print"], Cell[153175, 4540, 251, 5, 22, "Print"], Cell[153429, 4547, 481, 10, 24, "Print"], Cell[153913, 4559, 445, 9, 22, "Print"], Cell[154361, 4570, 350, 7, 22, "Print"], Cell[154714, 4579, 421, 9, 24, "Print"], Cell[155138, 4590, 174, 3, 24, "Print"], Cell[155315, 4595, 189, 4, 24, "Print"], Cell[155507, 4601, 462, 9, 24, "Print"], Cell[155972, 4612, 906, 20, 24, "Print"], Cell[156881, 4634, 222, 5, 22, "Print"], Cell[157106, 4641, 54, 1, 22, "Print"], Cell[157163, 4644, 251, 5, 22, "Print"], Cell[157417, 4651, 466, 9, 22, "Print"], Cell[157886, 4662, 447, 9, 22, "Print"], Cell[158336, 4673, 351, 7, 22, "Print"], Cell[158690, 4682, 415, 9, 22, "Print"], Cell[159108, 4693, 175, 3, 22, "Print"], Cell[159286, 4698, 177, 3, 22, "Print"], Cell[159466, 4703, 461, 9, 22, "Print"], Cell[159930, 4714, 902, 20, 22, "Print"], Cell[160835, 4736, 54, 1, 22, "Print"], Cell[160892, 4739, 251, 5, 22, "Print"], Cell[161146, 4746, 469, 9, 22, "Print"], Cell[161618, 4757, 447, 9, 22, "Print"], Cell[162068, 4768, 351, 7, 22, "Print"], Cell[162422, 4777, 419, 9, 22, "Print"], Cell[162844, 4788, 188, 4, 24, "Print"], Cell[163035, 4794, 187, 4, 22, "Print"], Cell[163225, 4800, 460, 9, 24, "Print"], Cell[163688, 4811, 903, 20, 24, "Print"], Cell[164594, 4833, 222, 5, 22, "Print"], Cell[164819, 4840, 54, 1, 22, "Print"], Cell[164876, 4843, 251, 5, 22, "Print"], Cell[165130, 4850, 474, 9, 22, "Print"], Cell[165607, 4861, 447, 9, 22, "Print"], Cell[166057, 4872, 350, 7, 22, "Print"], Cell[166410, 4881, 419, 9, 22, "Print"], Cell[166832, 4892, 175, 3, 22, "Print"], Cell[167010, 4897, 176, 3, 22, "Print"], Cell[167189, 4902, 473, 9, 22, "Print"], Cell[167665, 4913, 910, 20, 22, "Print"], Cell[168578, 4935, 54, 1, 22, "Print"], Cell[168635, 4938, 251, 5, 22, "Print"], Cell[168889, 4945, 481, 9, 22, "Print"], Cell[169373, 4956, 447, 9, 22, "Print"], Cell[169823, 4967, 351, 7, 22, "Print"], Cell[170177, 4976, 428, 9, 22, "Print"], Cell[170608, 4987, 187, 4, 24, "Print"], Cell[170798, 4993, 175, 3, 22, "Print"], Cell[170976, 4998, 474, 9, 24, "Print"], Cell[171453, 5009, 913, 20, 24, "Print"], Cell[172369, 5031, 54, 1, 22, "Print"], Cell[172426, 5034, 251, 5, 22, "Print"], Cell[172680, 5041, 482, 9, 24, "Print"], Cell[173165, 5052, 447, 9, 22, "Print"], Cell[173615, 5063, 351, 7, 22, "Print"], Cell[173969, 5072, 432, 9, 24, "Print"], Cell[174404, 5083, 188, 4, 24, "Print"], Cell[174595, 5089, 176, 3, 24, "Print"], Cell[174774, 5094, 472, 9, 24, "Print"], Cell[175249, 5105, 912, 20, 24, "Print"], Cell[176164, 5127, 222, 5, 22, "Print"], Cell[176389, 5134, 54, 1, 22, "Print"], Cell[176446, 5137, 251, 5, 22, "Print"], Cell[176700, 5144, 479, 9, 22, "Print"], Cell[177182, 5155, 453, 9, 22, "Print"], Cell[177638, 5166, 351, 7, 22, "Print"], Cell[177992, 5175, 427, 9, 22, "Print"], Cell[178422, 5186, 175, 3, 22, "Print"], Cell[178600, 5191, 177, 3, 22, "Print"], Cell[178780, 5196, 471, 9, 22, "Print"], Cell[179254, 5207, 905, 20, 22, "Print"], Cell[180162, 5229, 54, 1, 22, "Print"], Cell[180219, 5232, 251, 5, 22, "Print"], Cell[180473, 5239, 479, 9, 22, "Print"], Cell[180955, 5250, 453, 9, 22, "Print"], Cell[181411, 5261, 350, 7, 22, "Print"], Cell[181764, 5270, 427, 9, 22, "Print"], Cell[182194, 5281, 187, 4, 24, "Print"], Cell[182384, 5287, 173, 3, 22, "Print"], Cell[182560, 5292, 473, 9, 24, "Print"], Cell[183036, 5303, 912, 20, 24, "Print"], Cell[183951, 5325, 54, 1, 22, "Print"], Cell[184008, 5328, 251, 5, 22, "Print"], Cell[184262, 5335, 481, 9, 24, "Print"], Cell[184746, 5346, 453, 9, 22, "Print"], Cell[185202, 5357, 351, 7, 22, "Print"], Cell[185556, 5366, 433, 9, 24, "Print"], Cell[185992, 5377, 175, 3, 24, "Print"], Cell[186170, 5382, 176, 3, 24, "Print"], Cell[186349, 5387, 475, 9, 24, "Print"], Cell[186827, 5398, 915, 20, 24, "Print"], Cell[187745, 5420, 222, 5, 22, "Print"], Cell[187970, 5427, 2342, 41, 120, "Print"], Cell[190315, 5470, 187, 3, 37, "Message"], Cell[190505, 5475, 26416, 646, 198, 7581, 407, "GraphicsData", "PostScript", \ "Graphics", Evaluatable->False], Cell[216924, 6123, 24795, 642, 198, 7694, 424, "GraphicsData", "PostScript", \ "Graphics", Evaluatable->False] }, Open ]] }, Open ]] } ] *) (******************************************************************* End of Mathematica Notebook file. *******************************************************************)