(************** 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[ 72967, 2097]*) (*NotebookOutlinePosition[ 74015, 2130]*) (* CellTagsIndexPosition[ 73971, 2126]*) (*WindowFrame->Normal*) Notebook[{ Cell[TextData[ "FormTangentStiff2DTwoNodeBar[XY1_,XY2_,uXY1_,uXY2_,Em_,A0_,s0_]:=\n \ Module[{X1,Y1,X2,Y2,X21,Y21,uX1,uY1,uX2,uY2,\n \ uX21,uY21,L0,L,e,s,ax,ay,Ke}, \n {X1,Y1}=XY1; {X2,Y2}=XY2; X21=X2-X1; \ Y21=Y2-Y1; \n {uX1,uY1}=uXY1; {uX2,uY2}=uXY2; uX21=uX2-uX1; uY21=uY2-uY1;\n\ L0=Sqrt[X21^2+Y21^2]; L=Sqrt[(X21+uX21)^2+(Y21+uY21)^2];\n \ e=(L-L0)*(L+L0)/(2*L0^2); s=s0+Em*e; \n ax=(X21+uX21)/L0; \ ay=(Y21+uY21)/L0;\n Ke=(Em*A0/L0)*{{ ax*ax, ax*ay,-ax*ax,-ax*ay},\n \ { ax*ay, ay*ay,-ay*ax,-ay*ay},\n {-ax*ax,-ay*ax, \ ax*ax, ay*ax},\n {-ax*ay,-ay*ay, ay*ax, ay*ay}} +\n \ (A0*s/L0)*{{1,0,-1,0},{0,1,0,-1},{-1,0,1,0},{0,-1,0,1}};\n \ Return[Simplify[Ke]]\n ];\n \n\ Ke=FormTangentStiff2DTwoNodeBar[{0,-S+H},{0,H},{0,0},{0,uY},Em,A0,0];\n\ Ke=Simplify[Ke/.{Sqrt[S^2]->S,(1/Sqrt[S^2])->1/S}];\n\ Print[\"Ke=\",Ke//InputForm]; \n\ Ke=FormTangentStiff2DTwoNodeBar[{-4,0},{0,3},{0,0},{0,0}, Em,A0,s0];\n\ (*Print[\"Ke=\",Ke];*)\nPrint[Eigenvalues[Ke]];"], "Input", CellFrame->True, CellMargins->{{10, 60}, {Inherited, Inherited}}, CellLabelMargins->{{14, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, Background->RGBColor[1, 0.647715, 0.521981]], Cell[TextData[ "MergeElemIntoMasterStiff[Ke_,eftab_,Km_]:= \n \ Module[{i,j,ii,jj,neldof,K}, K=Km;\n neldof=Dimensions[Ke][[1]]; \n \ For[i=1, i<=neldof, i++, ii=eftab[[i]];\n For[j=i, j<=neldof, j++, \ jj=eftab[[j]];\n If [ii>0 && jj>0, \n \ K[[jj,ii]]=K[[ii,jj]]+=Ke[[i,j]]] \n ]\n ]; Return[K]\n \ ];\n \nAssembleMasterStiffOfShallowArch[S_,H_,Em_,A0_,uX_,uY_]:=\n\ Module[{K1,K2,K},\n K=Table[0,{2},{2}]; \n \ K1=FormTangentStiff2DTwoNodeBar[{-S/2,0},{0,H},\n \ {0,0},{uX,uY},Em,A0,0]; Print[\"K1=\",K1//MatrixForm];\n \ K=MergeElemIntoMasterStiff[K1,{0,0,1,2},K];\n \ K2=FormTangentStiff2DTwoNodeBar[{0,H},{S/2,0},\n \ {uX,uY},{0,0},Em,A0,0]; Print[\"K2=\",K2//MatrixForm];\n \ K=MergeElemIntoMasterStiff[K2,{1,2,0,0},K];\n Return[Simplify[K]]\n \ ];\n \nClearAll[S,H,Em,A0,uX,uY]; S=2; H=1; Em=2; A0=Sqrt[2]; uX=0; \ uY=0; \nK=Simplify[AssembleMasterStiffOfShallowArch[S,H,Em,A0,uX,uY]];\n\ Print[\"Master stiffness matrix = \",K//InputForm];\n(*detK=Simplify[Det[K]]; \ Print[\"det of K=\",detK//InputForm];*) "], "Input", CellFrame->True, CellMargins->{{10, 60}, {Inherited, Inherited}}, CellLabelMargins->{{14, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, Background->RGBColor[1, 0.647715, 0.521981]], Cell[TextData[ "\nK = {{(4*A0*Em*(S^2 + 6*uX^2 + 4*H*uY + 2*uY^2))/(4*H^2 + S^2)^(3/2), \n \ (16*A0*Em*uX*(H + uY))/(4*H^2 + S^2)^(3/2)}, \n {(16*A0*Em*uX*(H + \ uY))/(4*H^2 + S^2)^(3/2), \n (8*A0*Em*(2*H^2 + uX^2 + 6*H*uY + \ 3*uY^2))/(4*H^2 + S^2)^(3/2)}};\nPrint[K//FortranForm];"], "Input", AspectRatioFixed->True], Cell[TextData[ "\nsol=Solve[S^2 + 4*H*uY + 2*uY^2==0,uY];\nPrint[sol];"], "Input", AspectRatioFixed->True], Cell[TextData[ "MergeElemIntoMasterStiff[K_,eftab_,Km_]:= \n Module[{i,j,ii,jj,neldof,K}, \ K=Km;\n neldof=Dimensions[K][[1]];\n For[i=1, i<=neldof, i++, \ ii=eftab[[i]];\n For[j=i, j<=neldof, j++, jj=eftab[[j]];\n \ If [ii>0 && jj>0, \n K[[jj,ii]]=K[[ii,jj]]+=K[[i,j]]] \n \ ]\n ]; Return[K]\n ];\n \n\ AssembleMasterStiffOfShallowArch[S_,H_,Em_,A0_,uX_,uY_]:=\n \ Module[{K1,K2,K3,K},\n K=Table[0,{2},{2}]; \n \ K1=FormTangentStiff2DTwoNodeBar[{-S/2,0},{0,H},\n \ {0,0},{uX,uY},Em,A0,0]; (*Print[K1//TableForm];*)\n \ K=MergeElemIntoMasterStiff[K1,{0,0,1,2},K];\n \ K2=FormTangentStiff2DTwoNodeBar[{0,H},{S/2,0},\n \ {uX,uY},{0,0},Em,A0,0]; (*Print[K2//TableForm];*)\n \ K=MergeElemIntoMasterStiff[K2,{1,2,0,0},K];\n \ K3=FormTangentStiff2DTwoNodeBar[{0,-S+H},{0,H},\n \ {0,0},{uX,uY},Em,A0,0]; (*Print[K2//TableForm];*)\n \ K=MergeElemIntoMasterStiff[K3,{0,0,1,2},K];\n Return[Simplify[K]]\n \ ];\n \nClearAll[S,H,Em,A0,uX,uY]; \n\ K=AssembleMasterStiffOfShallowArch[2,H,1,1,0,uY];\nK=Simplify[Expand[K]];\n\ Print[\"Master stiffness matrix = \",K//InputForm];"], "Input", CellFrame->True, CellMargins->{{11, 56}, {Inherited, Inherited}}, AspectRatioFixed->True, Background->RGBColor[1, 0.647715, 0.521981]], Cell[CellGroupData[{ Cell[TextData[ "Kxx=2/(1 + H^2)^(3/2) + uY/4 + (2*H*uY)/(1 + H^2)^(3/2) + uY^2/16 + \n \ uY^2/(1 + H^2)^(3/2);\na=Coefficient[Kxx,uY,0]; b=Coefficient[Kxx,uY,1]; \ c=Coefficient[Kxx,uY,2];\nd=Simplify[b^2-4*a*c]; Print[\"d=\",d//InputForm]; \ \nPlot[d,{H,0,2},PlotRange->{0,1}];\nPrint[FindRoot[d==0,{H,1.2}]];"], "Input",\ CellMargins->{{10, 62}, {Inherited, Inherited}}, AspectRatioFixed->True, Background->RGBColor[1, 0.647715, 0.521981]], Cell[GraphicsData["PostScript", "\<\ %! %%Creator: Mathematica %%AspectRatio: .61803 MathPictureStart %% Graphics /Courier findfont 10 scalefont setfont % Scaling calculations 0.0238095 0.47619 0 0.618034 [ [(0)] .02381 0 0 2 Msboxa [(0.5)] .2619 0 0 2 Msboxa [(1)] .5 0 0 2 Msboxa [(1.5)] .7381 0 0 2 Msboxa [(2)] .97619 0 0 2 Msboxa [(0.2)] .01131 .12361 1 0 Msboxa [(0.4)] .01131 .24721 1 0 Msboxa [(0.6)] .01131 .37082 1 0 Msboxa [(0.8)] .01131 .49443 1 0 Msboxa [(1)] .01131 .61803 1 0 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 .02381 0 m .02381 .00625 L s P [(0)] .02381 0 0 2 Mshowa p .002 w .2619 0 m .2619 .00625 L s P [(0.5)] .2619 0 0 2 Mshowa p .002 w .5 0 m .5 .00625 L s P [(1)] .5 0 0 2 Mshowa p .002 w .7381 0 m .7381 .00625 L s P [(1.5)] .7381 0 0 2 Mshowa p .002 w .97619 0 m .97619 .00625 L s P [(2)] .97619 0 0 2 Mshowa p .001 w .07143 0 m .07143 .00375 L s P p .001 w .11905 0 m .11905 .00375 L s P p .001 w .16667 0 m .16667 .00375 L s P p .001 w .21429 0 m .21429 .00375 L s P p .001 w .30952 0 m .30952 .00375 L s P p .001 w .35714 0 m .35714 .00375 L s P p .001 w .40476 0 m .40476 .00375 L s P p .001 w .45238 0 m .45238 .00375 L s P p .001 w .54762 0 m .54762 .00375 L s P p .001 w .59524 0 m .59524 .00375 L s P p .001 w .64286 0 m .64286 .00375 L s P p .001 w .69048 0 m .69048 .00375 L s P p .001 w .78571 0 m .78571 .00375 L s P p .001 w .83333 0 m .83333 .00375 L s P p .001 w .88095 0 m .88095 .00375 L s P p .001 w .92857 0 m .92857 .00375 L s P p .002 w 0 0 m 1 0 L s P p .002 w .02381 .12361 m .03006 .12361 L s P [(0.2)] .01131 .12361 1 0 Mshowa p .002 w .02381 .24721 m .03006 .24721 L s P [(0.4)] .01131 .24721 1 0 Mshowa p .002 w .02381 .37082 m .03006 .37082 L s P [(0.6)] .01131 .37082 1 0 Mshowa p .002 w .02381 .49443 m .03006 .49443 L s P [(0.8)] .01131 .49443 1 0 Mshowa p .002 w .02381 .61803 m .03006 .61803 L s P [(1)] .01131 .61803 1 0 Mshowa p .001 w .02381 .02472 m .02756 .02472 L s P p .001 w .02381 .04944 m .02756 .04944 L s P p .001 w .02381 .07416 m .02756 .07416 L s P p .001 w .02381 .09889 m .02756 .09889 L s P p .001 w .02381 .14833 m .02756 .14833 L s P p .001 w .02381 .17305 m .02756 .17305 L s P p .001 w .02381 .19777 m .02756 .19777 L s P p .001 w .02381 .22249 m .02756 .22249 L s P p .001 w .02381 .27193 m .02756 .27193 L s P p .001 w .02381 .29666 m .02756 .29666 L s P p .001 w .02381 .32138 m .02756 .32138 L s P p .001 w .02381 .3461 m .02756 .3461 L s P p .001 w .02381 .39554 m .02756 .39554 L s P p .001 w .02381 .42026 m .02756 .42026 L s P p .001 w .02381 .44498 m .02756 .44498 L s P p .001 w .02381 .46971 m .02756 .46971 L s P p .001 w .02381 .51915 m .02756 .51915 L s P p .001 w .02381 .54387 m .02756 .54387 L s P p .001 w .02381 .56859 m .02756 .56859 L s P p .001 w .02381 .59331 m .02756 .59331 L s P p .002 w .02381 0 m .02381 .61803 L s P P 0 0 m 1 0 L 1 .61803 L 0 .61803 L closepath clip newpath p p .004 w s s s s s s s s s s s s s s s s s s s s s s .56182 0 m .57937 .0322 L s .57937 .0322 m .59921 .06243 L .61905 .08735 L .63889 .10773 L .65873 .12424 L .67857 .13747 L .69841 .14791 L .71825 .156 L .7381 .16212 L .75794 .16657 L .76786 .16826 L .77778 .16963 L .7877 .17072 L .79762 .17154 L .80754 .17212 L .8125 .17233 L .81498 .17242 L .81746 .17249 L .81994 .17255 L .82242 .1726 L .82366 .17262 L .8249 .17263 L .82614 .17265 L .82738 .17266 L .82862 .17266 L .82986 .17267 L .8311 .17267 L .83234 .17267 L .83358 .17267 L .83482 .17266 L .8373 .17264 L .83854 .17263 L .83978 .17262 L .84226 .17258 L .84722 .17247 L .85218 .17233 L .85714 .17215 L .86706 .1717 L .87698 .17113 L .89683 .16967 L .93651 .16581 L .97619 .1611 L s P P % End of Graphics MathPictureEnd \ \>"], "Graphics", CellMargins->{{Inherited, 57}, {Inherited, Inherited}}, Evaluatable->False, AspectRatioFixed->True, ImageSize->{477, 293}, ImageMargins->{{34, Inherited}, {Inherited, Inherited}}, Background->RGBColor[0, 1, 1], ImageCache->GraphicsData["Bitmap", "\<\ CF5dJ6E]HGAYHf4PAg9QL6QYHg/3o`030003o`?o0`3o`030003 o`?o0h3o`030003o`?o0l3o`030003o`?o0<@3o`0090?o 00<000?o0ol0l0?o00<000?o0ol0``?o000T0ol00`000ol3o`3a0ol00`000ol3o`320ol002@3o`03 0003o`?o0?43o`030003o`?o0<83o`0090?o00<000?o0ol0lP?o00<000?o0ol0`@?o000T0ol00`00 0ol3o`3c0ol00`000ol3o`300ol002@3o`030003o`?o0?@3o`030003o`?o0;l3o`0090?o0P00m@?o 00<000?o0ol0_`?o000T0ol00`000ol3o`3e0ol00`000ol3o`2n0ol002@3o`030003o`?o0?H3o`03 0003o`?o0;d3o`0090?o00<000?o0ol0mP?o00<000?o0ol0_@?o000T0ol00`000ol3o`3g0ol00`00 0ol3o`2l0ol002@3o`030003o`?o0?P3o`030003o`?o0;/3o`0090?o00<000?o0ol0n0?o00<000?o 0ol0^`?o000T0ol00`000ol3o`3i0ol00`000ol3o`2j0ol002@3o`030003o`?o0?X3o`030003o`?o 0;T3o`0090?o00<000?o0ol0n`?o00<000?o0ol0^0?o000T0ol2003m0ol00`000ol3o`2g0ol002@3 o`030003o`?o0?d3o`030003o`?o0;H3o`0090?o00<000?o0ol0o@?o00<000?o0ol0]P?o000T0ol0 0`000ol3o`3n0ol00`000ol3o`2e0ol002@3o`030003o`?o0?l3o`030003o`?o0;@3o`0090?o00<0 00?o0ol0o`?o0@?o00<000?o0ol0/`?o000T0ol00`000ol3o`3o0ol20ol00`000ol3o`2b0ol002@3 o`030003o`?o0?l3o`83o`030003o`?o0;83o`0090?o00<000?o0ol0o`?o0`?o00<000?o0ol0/@?o 000T0ol00`000ol3o`3o0ol40ol00`000ol3o`2`0ol002@3o`030003o`?o0?l3o`D3o`030003o`?o 0:l3o`0090?o0P00o`?o1`?o00<000?o0ol0[P?o000T0ol00`000ol3o`3o0ol70ol00`000ol3o`2] 0ol002@3o`030003o`?o0?l3o`P3o`030003o`?o0:`3o`0090?o00<000?o0ol0o`?o2@?o00<000?o 0ol0Z`?o000T0ol00`000ol3o`3o0ol:0ol2002[0ol002@3o`030003o`?o0?l3o``3o`030003o`?o 0:P3o`0090?o00<000?o0ol0o`?o3@?o00<000?o0ol0Y`?o000T0ol00`000ol3o`3o0ol>0ol2002W 0ol002@3o`030003o`?o0?l3oa03o`030003o`?o0:@3o`003P?o0P0010?o00@000?o0ol3o`D000L3 o`030003o`?o0?l3oa43o`030003o`?o0:<3o`003@?o00@000?o0ol000P3o`040003o`?o00070ol0 0`000ol3o`3o0olB0ol00`000ol3o`2R0ol000d3o`040003o`?o00090ol00`000ol3o`070ol4003o 0olB0ol2002R0ol000d3o`040003o`?o000:0ol00`000ol3o`060ol00`000ol3o`3o0olE0ol00`00 0ol3o`2O0ol000d3o`040003o`?o00070ol01@000ol3o`?o00001`?o00<000?o0ol0o`?o5P?o0P00 W`?o000=0ol010000ol3o`001`?o00D000?o0ol3o`0000L3o`030003o`?o0?l3oaP3o`030003o`?o 09`3o`003P?o0P002@?o0`0020?o00<000?o0ol0o`?o6@?o0P00W0?o000T0ol00`000ol3o`3o0olK 0ol00`000ol3o`2I0ol002@3o`030003o`?o0?l3oa`3o`8009T3o`0090?o00<000?o0ol0o`?o7P?o 00<000?o0ol0UP?o000T0ol00`000ol3o`3o0olO0ol2002F0ol002@3o`030003o`?o0?l3ob43o`80 09@3o`0090?o00<000?o0ol0o`?o8`?o00<000?o0ol0T@?o000T0ol2003o0olU0ol3002@0ol002@3 o`030003o`?o0?l3obL3o`<008d3o`0090?o00<000?o0ol0o`?o:P?o0`00RP?o000T0ol00`000ol3 o`3o0ol]0ol300270ol002@3o`030003o`?o0?l3oc03o`<008@3o`0090?o00<000?o0ol0o`?o<`?o 0`00I0?o1@0060?o000T0ol00`000ol3o`3o0olf0ol3001H0ol9000M0ol002@3o`030003o`?o0?l3 ocT3o`@004`3o`P002H3o`0090?o00<000?o0ol0o`?o?@?o1@00?P?o2@00;P?o000T0ol00`000ol3 o`3o0om20ol8000]0ol9000g0ol002@3o`030003o`?o0?l3odX3obd00403o`0090?o0P00o`?o^0?o 000T0ol00`000ol3o`3o0ong0ol002@3o`030003o`?o0?l3okL3o`0090?o00<000?o0ol0o`?o]`?o 000T0ol00`000ol3o`3o0ong0ol002@3o`030003o`?o0?l3okL3o`0090?o00<000?o0ol0o`?o]`?o 000T0ol00`000ol3o`3o0ong0ol002@3o`030003o`?o0?l3okL3o`0090?o00<000?o0ol0o`?o]`?o 000T0ol2003o0onh0ol002@3o`030003o`?o0?l3okL3o`0090?o00<000?o0ol0o`?o]`?o000T0ol0 0`000ol3o`3o0ong0ol002@3o`030003o`?o0?l3okL3o`0090?o00<000?o0ol0o`?o]`?o000T0ol0 0`000ol3o`3o0ong0ol002@3o`030003o`?o0?l3okL3o`0090?o00<000?o0ol0o`?o]`?o000T0ol0 0`000ol3o`3o0ong0ol002@3o`030003o`?o0?l3okL3o`0090?o0P00o`?o^0?o000T0ol00`000ol3 o`3o0ong0ol002@3o`030003o`?o0?l3okL3o`0090?o00<000?o0ol0o`?o]`?o000T0ol00`000ol3 o`3o0ong0ol002@3o`030003o`?o0?l3okL3o`0090?o00<000?o0ol0o`?o]`?o000T0ol00`000ol3 o`3o0ong0ol000h3o`8000@3o`030003o`?o00<3o`<000L3o`030003o`?o0?l3okL3o`003@?o00@0 00?o0ol000X3o`030003o`?o00H3o`030003o`?o0?l3okL3o`003@?o00@000?o0ol000L3o`D000L3 o`030003o`?o0?l3okL3o`003@?o00@000?o0ol000L3o`040003o`?o00080ol4003o0onf0ol000d3 o`040003o`?o00080ol00`000ol000080ol00`000ol3o`3o0ong0ol000d3o`040003o`?o00090ol2 00080ol00`000ol3o`3o0ong0ol000h3o`8000/3o`030003o`?o00H3o`030003o`?o0?l3okL3o`00 90?o00<000?o0ol0o`?o]`?o000T0ol00`000ol3o`3o0ong0ol002@3o`030003o`?o0?l3okL3o`00 90?o00<000?o0ol0o`?o]`?o000T0ol00`000ol3o`3o0ong0ol002@3o`030003o`?o0?l3okL3o`00 90?o00<000?o0ol0o`?o]`?o000T0ol2003o0onh0ol002@3o`030003o`?o0?l3okL3o`0090?o00<0 00?o0ol0o`?o]`?o000T0ol00`000ol3o`3o0ong0ol002@3o`030003o`?o0?l3okL3o`0090?o00<0 00?o0ol0o`?o]`?o000T0ol00`000ol3o`3o0ong0ol002@3o`030003o`?o0?l3okL3o`0090?o00<0 00?o0ol0o`?o]`?o000T0ol00`000ol3o`3o0ong0ol002@3o`030003o`?o0?l3okL3o`0090?o0P00 o`?o^0?o000T0ol00`000ol3o`3o0ong0ol002@3o`030003o`?o0?l3okL3o`0090?o00<000?o0ol0 o`?o]`?o000T0ol00`000ol3o`3o0ong0ol002@3o`030003o`?o0?l3okL3o`0090?o00<000?o0ol0 o`?o]`?o000T0ol00`000ol3o`3o0ong0ol002@3o`030003o`?o0?l3okL3o`0090?o00<000?o0ol0 o`?o]`?o000T0ol2003o0onh0ol002@3o`030003o`?o0?l3okL3o`0090?o00<000?o0ol0o`?o]`?o 000T0ol00`000ol3o`3o0ong0ol002@3o`030003o`?o0?l3okL3o`0090?o00<000?o0ol0o`?o]`?o 000T0ol00`000ol3o`3o0ong0ol002@3o`030003o`?o0?l3okL3o`0090?o00<000?o0ol0o`?o]`?o 000T0ol00`000ol3o`3o0ong0ol002@3o`030003o`?o0?l3okL3o`0090?o0P00o`?o^0?o000T0ol0 0`000ol3o`3o0ong0ol002@3o`030003o`?o0?l3okL3o`0090?o00<000?o0ol0o`?o]`?o000T0ol0 0`000ol3o`3o0ong0ol002@3o`030003o`?o0?l3okL3o`0090?o00<000?o0ol0o`?o]`?o000T0ol0 0`000ol3o`3o0ong0ol000h3o`8000@3o`030003o`?o0083o`<000P3o`030003o`?o0?l3okL3o`00 3@?o00@000?o0ol000L3o`050003o`?o0ol000070ol00`000ol3o`3o0ong0ol000d3o`040003o`?o 00070ol01@000ol3o`?o00001`?o00<000?o0ol0o`?o]`?o000=0ol010000ol3o`001`?o100020?o 1000o`?o]P?o000=0ol010000ol3o`001`?o00<000?o0ol02@?o00<000?o0ol0o`?o]`?o000=0ol0 10000ol3o`0020?o00<000?o0ol020?o00<000?o0ol0o`?o]`?o000>0ol2000:0ol200080ol00`00 0ol3o`3o0ong0ol002@3o`030003o`?o0?l3okL3o`0090?o00<000?o0ol0o`?o]`?o000T0ol00`00 0ol3o`3o0ong0ol002@3o`030003o`?o0?l3okL3o`0090?o00<000?o0ol0o`?o]`?o000T0ol00`00 0ol3o`3o0ong0ol002@3o`030003o`?o0?l3okL3o`0090?o0P00o`?o^0?o000T0ol00`000ol3o`3o 0ong0ol002@3o`030003o`?o0?l3okL3o`0090?o00<000?o0ol0o`?o]`?o000T0ol00`000ol3o`3o 0ong0ol002@3o`030003o`?o0?l3okL3o`0090?o00<000?o0ol0o`?o]`?o000T0ol00`000ol3o`3o 0ong0ol002@3o`030003o`?o0?l3okL3o`0090?o00<000?o0ol0o`?o]`?o000T0ol00`000ol3o`3o 0ong0ol002@3o`800?l3okP3o`0090?o00<000?o0ol0o`?o]`?o000T0ol00`000ol3o`3o0ong0ol0 02@3o`030003o`?o0?l3okL3o`0090?o00<000?o0ol0o`?o]`?o000T0ol00`000ol3o`3o0ong0ol0 02@3o`030003o`?o0?l3okL3o`0090?o00<000?o0ol0o`?o]`?o000T0ol00`000ol3o`3o0ong0ol0 02@3o`030003o`?o0?l3okL3o`0090?o0P00o`?o^0?o000T0ol00`000ol3o`3o0ong0ol002@3o`03 0003o`?o0?l3okL3o`0090?o00<000?o0ol0o`?o]`?o000T0ol00`000ol3o`3o0ong0ol002@3o`03 0003o`?o0?l3okL3o`0090?o00<000?o0ol0o`?o]`?o000T0ol00`000ol3o`3o0ong0ol002@3o`03 0003o`?o0?l3okL3o`0090?o00<000?o0ol0o`?o]`?o000T0ol00`000ol3o`3o0ong0ol002@3o`80 0?l3okP3o`0090?o00<000?o0ol0o`?o]`?o000T0ol00`000ol3o`3o0ong0ol002@3o`030003o`?o 0?l3okL3o`0090?o00<000?o0ol0o`?o]`?o000T0ol00`000ol3o`3o0ong0ol002@3o`030003o`?o 0?l3okL3o`0090?o00<000?o0ol0o`?o]`?o000>0ol200040ol00`000ol3o`020ol300080ol00`00 0ol3o`3o0ong0ol000d3o`040003o`?o00070ol01@000ol3o`?o00001`?o00<000?o0ol0o`?o]`?o 000=0ol010000ol3o`001`?o00D000?o0ol3o`0000L3o`030003o`?o0?l3okL3o`003@?o00@000?o 0ol000P3o`<000P3o`@00?l3okH3o`003@?o00@000?o0ol000L3o`050003o`?o0ol000070ol00`00 0ol3o`3o0ong0ol000d3o`040003o`?o00070ol01@000ol3o`?o00001`?o00<000?o0ol0o`?o]`?o 000>0ol200090ol300080ol00`000ol3o`3o0ong0ol002@3o`030003o`?o0?l3okL3o`0090?o00<0 00?o0ol0o`?o]`?o000T0ol00`000ol3o`3o0ong0ol002@3o`030003o`?o0?l3okL3o`0090?o00<0 00?o0ol0o`?o]`?o000T0ol00`000ol3o`3o0ong0ol002@3o`030003o`?o0?l3okL3o`0090?o0P00 o`?o^0?o000T0ol00`000ol3o`3o0ong0ol002@3o`030003o`?o0?l3okL3o`0090?o00<000?o0ol0 o`?o]`?o000T0ol00`000ol3o`3o0ong0ol002@3o`030003o`?o0?l3okL3o`0090?o00<000?o0ol0 o`?o]`?o000T0ol00`000ol3o`3o0ong0ol002@3o`030003o`?o0?l3okL3o`0090?o00<000?o0ol0 o`?o]`?o000T0ol00`000ol3o`3o0ong0ol002@3o`800?l3okP3o`0090?o00<000?o0ol0o`?o]`?o 000T0ol00`000ol3o`3o0ong0ol002@3o`030003o`?o0?l3okL3o`0090?o00<000?o0ol0o`?o]`?o 000T0ol00`000ol3o`3o0ong0ol002@3o`030003o`?o0?l3okL3o`0090?o00<000?o0ol0o`?o]`?o 000T0ol00`000ol3o`3o0ong0ol002@3o`030003o`?o0?l3okL3o`0090?o0P00o`?o^0?o000T0ol0 0`000ol3o`3o0ong0ol002@3o`030003o`?o0?l3okL3o`0090?o00<000?o0ol0o`?o]`?o000T0ol0 0`000ol3o`3o0ong0ol002@3o`030003o`?o0?l3okL3o`0090?o00<000?o0ol0o`?o]`?o000T0ol0 0`000ol3o`3o0ong0ol002@3o`030003o`?o0?l3okL3o`0090?o00<000?o0ol0o`?o]`?o000T0ol0 0`000ol3o`3o0ong0ol002@3o`800?l3okP3o`0090?o00<000?o0ol0o`?o]`?o000T0ol00`000ol3 o`3o0ong0ol002@3o`030003o`?o0?l3okL3o`0090?o00<000?o0ol0o`?o]`?o000T0ol00`000ol3 o`3o0ong0ol002@3o`030003o`?o0?l3okL3o`0090?o00<000?o0ol0o`?o]`?o000I0ol300080ol0 0`000ol3o`3o0ong0ol001X3o`030003o`?o00L3o`030003o`?o0?l3okL3o`006P?o00<000?o0ol0 1`?o00<000?o0ol0o`?o]`?o000J0ol00`000ol3o`070ol4003o0onf0ol001X3o`030003o`?o0?l3 ol43o`006@?o0P00o`?o``?o000J0ol00`000ol3o`3o0oo10ol00?l3omh3o`00o`?ogP?o0000\ \>"], ImageRangeCache->{{{0, 476}, {292, 0}} -> {-0.174391, -0.0611163, \ 0.00480695, 0.00370372}}] }, Open ]], Cell[TextData[ "ClearAll[S,H,Em,A0,uX,uY];\n\ eigval1[S_,H_,Em_,A0_,uX_,uY_]:=(2*(256*A0*Em*H^8 + 256*A0*Em*H^6*S^2 + \n \ 96*A0*Em*H^4*S^4 + 16*A0*Em*H^2*S^6 + A0*Em*S^8 + \n \ 512*A0*Em*H^6*uX^2 + 384*A0*Em*H^4*S^2*uX^2 + \n 96*A0*Em*H^2*S^4*uX^2 \ + 8*A0*Em*S^6*uX^2 + \n 1024*A0*Em*H^7*uY + 768*A0*Em*H^5*S^2*uY + \n \ 192*A0*Em*H^3*S^4*uY + 16*A0*Em*H*S^6*uY + \n 512*A0*Em*H^6*uY^2 \ + 384*A0*Em*H^4*S^2*uY^2 + \n 96*A0*Em*H^2*S^4*uY^2 + 8*A0*Em*S^6*uY^2 \ - \n (A0^2*Em^2*(4*H^2 + S^2)^6*\n (16*H^4 - 8*H^2*S^2 + S^4 \ + 32*H^2*uX^2 + 8*S^2*uX^2 + \n 16*uX^4 + 64*H^3*uY - 16*H*S^2*uY \ + 64*H*uX^2*uY + \n 96*H^2*uY^2 - 8*S^2*uY^2 + 32*uX^2*uY^2 + \n \ 64*H*uY^3 + 16*uY^4))^(1/2)))/(4*H^2 + S^2)^(9/2); \n\ eigval2[S_,H_,Em_,A0_,uX_,uY_]:=(2*(256*A0*Em*H^8 + 256*A0*Em*H^6*S^2 + \n \ 96*A0*Em*H^4*S^4 + 16*A0*Em*H^2*S^6 + A0*Em*S^8 + \n \ 512*A0*Em*H^6*uX^2 + 384*A0*Em*H^4*S^2*uX^2 + \n 96*A0*Em*H^2*S^4*uX^2 \ + \n 8*A0*Em*S^6*uX^2 + 1024*A0*Em*H^7*uY + \n \ 768*A0*Em*H^5*S^2*uY + 192*A0*Em*H^3*S^4*uY + \n 16*A0*Em*H*S^6*uY + \ 512*A0*Em*H^6*uY^2 + \n 384*A0*Em*H^4*S^2*uY^2 + 96*A0*Em*H^2*S^4*uY^2 \ + \n 8*A0*Em*S^6*uY^2 + \n (A0^2*Em^2*(4*H^2 + S^2)^6*\n \ (16*H^4 - 8*H^2*S^2 + S^4 + 32*H^2*uX^2 + 8*S^2*uX^2 + \n \ 16*uX^4 + 64*H^3*uY - 16*H*S^2*uY + 64*H*uX^2*uY + \n 96*H^2*uY^2 \ - 8*S^2*uY^2 + 32*uX^2*uY^2 + \n 64*H*uY^3 + \ 16*uY^4))^(1/2)))/(4*H^2 + S^2)^(9/2);\n \n \ eig1=eigval1[S,H,1,1,0,uY]; eig2=eigval2[S,H,1,1,0,uY];\n \ eig1=Simplify[eig1]; \n eig2=Simplify[eig2];\n \ Print[\"eig1P=\",eig1//InputForm]; Print[\"eig2P=\",eig2//InputForm];\n \ (*Plot[{eig1,eig2},{uY,0,-2*H},PlotRange->{{-2*H,0},{-0.5,0.5}}];*)\n \n \ uX=Sqrt[-S^2/2-2*H*uY-uY^2];\n eig1=eigval1[S,H,1,1,uX,uY]; \ eig2=eigval2[S,H,1,1,uX,uY];\n eig1=Simplify[eig1]; \n eig2=Simplify[eig2];\n \ Print[\"eig1B=\",eig1//InputForm]; Print[\"eig2B=\",eig2//InputForm];\n \ uYB1=Re[N[-H+Sqrt[H^2-S^2/2]]]; uYB2=Re[N[-H-Sqrt[H^2-S^2/2]]];\n \ (*Plot[{eig1,eig2},{uY,uYB1,uYB2},PlotRange->{{-2*H,0},{-0.5,0.5}}];*)"], "Input", CellFrame->True, CellMargins->{{10, 60}, {Inherited, Inherited}}, CellLabelMargins->{{14, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, Background->RGBColor[1, 0.647715, 0.521981]], Cell[TextData[ "ClearAll[H,S,Em,A0,uY];\ndetK=(32*A0^2*Em^2*(S^2 + 4*H*uY + 2*uY^2)*\n \ (2*H^2 + 6*H*uY + 3*uY^2))/(4*H^2 + S^2)^3;\nsol=Simplify[Solve[detK==0,uY]]; \ \nPrint[sol//InputForm];\nlambdaCP=(8*A0*Em*uY*(H + uY)*(2*H + uY))/(4*H^2 + \ S^2)^(3/2)/.sol;\nlambdaCP=Simplify[Expand[lambdaCP]]; \n\ Print[\"lambdaCP=\",lambdaCP//InputForm];"], "Input", CellFrame->True, CellMargins->{{10, 60}, {Inherited, Inherited}}, CellLabelMargins->{{14, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, Background->RGBColor[1, 0.647715, 0.521981]], Cell[CellGroupData[{ Cell[TextData[ "ClearAll[H,S,Em,A0,uY];\nlambdaL=(16*A0*Em*H^3)/(3*3^(1/2)*(4*H^2 + \ S^2)^(3/2));\nlambdaB=(2*2^(1/2)*A0*Em*S^2*(2*H^2 - S^2)^(1/2))/(4*H^2 + \ S^2)^(3/2);\nset={Em->1,A0->1,S->2};\n\ pL=Plot[lambdaL/.set,{H,0,5},DisplayFunction->Identity,\n \ AxesLabel->{\"H\",\"lambda at CPs\"},\n PlotLabel->\"lambdaCP vs H \ (S=2,E=A0=1)\"];\n\ pB=Plot[lambdaB/.set,{H,Sqrt[2.00001],5},DisplayFunction->Identity];\n\ Show[{pL,pB},DisplayFunction->$DisplayFunction];"], "Input", CellFrame->True, CellMargins->{{10, 60}, {Inherited, Inherited}}, CellLabelMargins->{{14, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, Background->RGBColor[1, 0.647715, 0.521981]], Cell[GraphicsData["PostScript", "\<\ %! %%Creator: Mathematica %%AspectRatio: .61803 MathPictureStart %% Graphics /Courier findfont 10 scalefont setfont % Scaling calculations 0.0238095 0.190476 0.0147151 1.6219 [ [(1)] .21429 .01472 0 2 Msboxa [(2)] .40476 .01472 0 2 Msboxa [(3)] .59524 .01472 0 2 Msboxa [(4)] .78571 .01472 0 2 Msboxa [(5)] .97619 .01472 0 2 Msboxa [(H)] 1.025 .01472 -1 0 Msboxa [(lambdaCP vs H \\(S=2,E=A0=1\\))] .5 .61803 0 -2 0 0 1 Mouter Mrotsboxa [(0.05)] .01131 .09581 1 0 Msboxa [(0.1)] .01131 .17691 1 0 Msboxa [(0.15)] .01131 .258 1 0 Msboxa [(0.2)] .01131 .3391 1 0 Msboxa [(0.25)] .01131 .42019 1 0 Msboxa [(0.3)] .01131 .50129 1 0 Msboxa [(0.35)] .01131 .58238 1 0 Msboxa [(lambda at CPs)] .02381 .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 .21429 .01472 m .21429 .02097 L s P [(1)] .21429 .01472 0 2 Mshowa p .002 w .40476 .01472 m .40476 .02097 L s P [(2)] .40476 .01472 0 2 Mshowa p .002 w .59524 .01472 m .59524 .02097 L s P [(3)] .59524 .01472 0 2 Mshowa p .002 w .78571 .01472 m .78571 .02097 L s P [(4)] .78571 .01472 0 2 Mshowa p .002 w .97619 .01472 m .97619 .02097 L s P [(5)] .97619 .01472 0 2 Mshowa p .001 w .0619 .01472 m .0619 .01847 L s P p .001 w .1 .01472 m .1 .01847 L s P p .001 w .1381 .01472 m .1381 .01847 L s P p .001 w .17619 .01472 m .17619 .01847 L s P p .001 w .25238 .01472 m .25238 .01847 L s P p .001 w .29048 .01472 m .29048 .01847 L s P p .001 w .32857 .01472 m .32857 .01847 L s P p .001 w .36667 .01472 m .36667 .01847 L s P p .001 w .44286 .01472 m .44286 .01847 L s P p .001 w .48095 .01472 m .48095 .01847 L s P p .001 w .51905 .01472 m .51905 .01847 L s P p .001 w .55714 .01472 m .55714 .01847 L s P p .001 w .63333 .01472 m .63333 .01847 L s P p .001 w .67143 .01472 m .67143 .01847 L s P p .001 w .70952 .01472 m .70952 .01847 L s P p .001 w .74762 .01472 m .74762 .01847 L s P p .001 w .82381 .01472 m .82381 .01847 L s P p .001 w .8619 .01472 m .8619 .01847 L s P p .001 w .9 .01472 m .9 .01847 L s P p .001 w .9381 .01472 m .9381 .01847 L s P [(H)] 1.025 .01472 -1 0 Mshowa p .002 w 0 .01472 m 1 .01472 L s P [(lambdaCP vs H \\(S=2,E=A0=1\\))] .5 .61803 0 -2 0 0 1 Mouter Mrotshowa p .002 w .02381 .09581 m .03006 .09581 L s P [(0.05)] .01131 .09581 1 0 Mshowa p .002 w .02381 .17691 m .03006 .17691 L s P [(0.1)] .01131 .17691 1 0 Mshowa p .002 w .02381 .258 m .03006 .258 L s P [(0.15)] .01131 .258 1 0 Mshowa p .002 w .02381 .3391 m .03006 .3391 L s P [(0.2)] .01131 .3391 1 0 Mshowa p .002 w .02381 .42019 m .03006 .42019 L s P [(0.25)] .01131 .42019 1 0 Mshowa p .002 w .02381 .50129 m .03006 .50129 L s P [(0.3)] .01131 .50129 1 0 Mshowa p .002 w .02381 .58238 m .03006 .58238 L s P [(0.35)] .01131 .58238 1 0 Mshowa p .001 w .02381 .03093 m .02756 .03093 L s P p .001 w .02381 .04715 m .02756 .04715 L s P p .001 w .02381 .06337 m .02756 .06337 L s P p .001 w .02381 .07959 m .02756 .07959 L s P p .001 w .02381 .11203 m .02756 .11203 L s P p .001 w .02381 .12825 m .02756 .12825 L s P p .001 w .02381 .14447 m .02756 .14447 L s P p .001 w .02381 .16069 m .02756 .16069 L s P p .001 w .02381 .19312 m .02756 .19312 L s P p .001 w .02381 .20934 m .02756 .20934 L s P p .001 w .02381 .22556 m .02756 .22556 L s P p .001 w .02381 .24178 m .02756 .24178 L s P p .001 w .02381 .27422 m .02756 .27422 L s P p .001 w .02381 .29044 m .02756 .29044 L s P p .001 w .02381 .30666 m .02756 .30666 L s P p .001 w .02381 .32288 m .02756 .32288 L s P p .001 w .02381 .35531 m .02756 .35531 L s P p .001 w .02381 .37153 m .02756 .37153 L s P p .001 w .02381 .38775 m .02756 .38775 L s P p .001 w .02381 .40397 m .02756 .40397 L s P p .001 w .02381 .43641 m .02756 .43641 L s P p .001 w .02381 .45263 m .02756 .45263 L s P p .001 w .02381 .46885 m .02756 .46885 L s P p .001 w .02381 .48507 m .02756 .48507 L s P p .001 w .02381 .51751 m .02756 .51751 L s P p .001 w .02381 .53372 m .02756 .53372 L s P p .001 w .02381 .54994 m .02756 .54994 L s P p .001 w .02381 .56616 m .02756 .56616 L s P p .001 w .02381 .5986 m .02756 .5986 L s P p .001 w .02381 .61482 m .02756 .61482 L s P [(lambda at CPs)] .02381 .61803 0 -4 Mshowa p .002 w .02381 0 m .02381 .61803 L s P P 0 0 m 1 0 L 1 .61803 L 0 .61803 L closepath clip newpath p p p .004 w .02381 .01472 m .02505 .01472 L .02629 .01472 L .02753 .01472 L .02877 .01473 L .03001 .01474 L .03125 .01475 L .03249 .01477 L .03373 .0148 L .03621 .01489 L .03745 .01494 L .03869 .01501 L .04117 .01518 L .04365 .01541 L .04613 .0157 L .04861 .01606 L .05357 .01701 L .05853 .01832 L .06349 .02001 L .07341 .02471 L .08333 .03128 L .09325 .0398 L .10317 .05023 L .14286 .10765 L .18254 .17851 L .22222 .24907 L .2619 .31196 L .30159 .36491 L .34127 .40832 L .38095 .44356 L .42063 .47212 L .46032 .49536 L .5 .51439 L .53968 .53008 L .57937 .54313 L .61905 .55406 L .65873 .56329 L .69841 .57114 L .7381 .57786 L .77778 .58366 L .81746 .58869 L .85714 .59307 L .89683 .59692 L .93651 .60031 L .97619 .60332 L s P P p p .004 w .29318 .01669 m .29407 .08607 L .29496 .11502 L .29674 .15494 L .29852 .18448 L .3003 .2085 L .30386 .24669 L .30741 .27654 L .31453 .32121 L .32164 .35311 L .32876 .37663 L .33587 .39414 L .34299 .40709 L .34654 .4122 L .3501 .41651 L .35366 .42013 L .35722 .42312 L .36077 .42554 L .36433 .42745 L .36789 .42889 L .36967 .42946 L .37145 .42992 L .37322 .43029 L .37411 .43044 L .375 .43056 L .37589 .43067 L .37678 .43076 L .37767 .43082 L .37856 .43086 L .37945 .43089 L .38034 .4309 L .38123 .43088 L .38212 .43085 L .38301 .4308 L .3839 .43073 L .38567 .43055 L .38745 .4303 L .38923 .42999 L .39279 .42919 L .3999 .42697 L .40702 .42402 L .43548 .40727 L .46394 .3861 L .49239 .36339 L .52085 .34068 L .54931 .31879 L .57777 .29811 L .60623 .27884 L .63469 .261 L .66315 .24456 L Mistroke .6916 .22947 L .72006 .21561 L .74852 .20291 L .77698 .19125 L .80544 .18056 L .8339 .17073 L .86236 .16169 L .89081 .15337 L .91927 .1457 L .94773 .13861 L .97619 .13206 L Mfstroke P P P % End of Graphics MathPictureEnd \ \>"], "Graphics", CellMargins->{{10, 60}, {Inherited, Inherited}}, Evaluatable->False, CellLabelMargins->{{14, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, ImageSize->{395, 241}, ImageMargins->{{0, Inherited}, {Inherited, Inherited}}, Background->RGBColor[0, 1, 1], ImageCache->GraphicsData["Bitmap", "\<\ CF5dJ6E]HGAYHf4PAg9QL6QYHgP?o1@00>P?o0`00?0?o0`00>P?o0`008P?o001_0ol00`000ol3o`0j 0ol010000ol3o`00>@?o00D000?o0ol3o`0003`3o`030003o`?o03P3o`050003o`?o0ol0000Q0ol0 06l3o`030003o`?o03/3o`030003o`?o03d3o`030003o`?o03L3o`D003d3o`030003o`?o01l3o`00 K`?o00<000?o0ol0?0?o00<000?o0ol0>P?o0P00>P?o00@000?o0ol003h3o`030003o`?o01l3o`00 K`?o00<000?o0ol0>@?o00D000?o0ol3o`0003d3o`030003o`?o03P3o`030003o`0003X3o`@00283 o`00KP?o0P00>`?o00D000?o0ol3o`0003T3o`050003o`?o0ol0000k0ol2000j0ol00`000ol3o`0S 0ol006l3o`030003o`?o03X3o`<003/3o`<003d3o`030003o`?o03P3o`D00243o`00o`?oS0?o003o 0on<0ol00383o`030003o`?o0?l3oeL3o`000ol00383o`03 0003o`?o0?l3odD3o`040003o`?o000>0ol002/3ool004L000P3o`@000h3o`00`?o00<000?o0ol0>`?o00<000?o0ol0>`?o00<000?o0ol03P?o00@000?o0ol000h3o`00`?o00<000?o0ol0oP?o000b0ol00`000ol3o`0G0ol00`000ol3o`0k 0ol00`000ol3o`3n0ol00383o`030003o`?o01P3o`030003o`?o03X3o`030003o`?o0?h3o`00@?o00<000?o0ol0oP?o000b0ol3000J0ol00`000ol3o`0h0ol0 0`000ol3o`3n0ol00383o`030003o`?o01X3o`030003o`?o03P3o`030003o`?o0?h3o`000?o00<000?o0ol070?o00<000?o0ol0Q@?o0P00 M@?o000b0ol00`000ol3o`0h0ol00`000ol3o`0L0ol00`000ol3o`230ol2001g0ol00383o`030003 o`?o03T3o`030003o`?o01/3o`030003o`?o0843o`8007T3o`00@?o00<000?o0ol06`?o 00<000?o0ol0O`?o0P00N`?o000b0ol00`000ol3o`0j0ol00`000ol3o`0J0ol00`000ol3o`1m0ol2 001m0ol00383o`030003o`?o03X3o`030003o`?o01/3o`030003o`?o07X3o`8007l3o`00`?o00<000?o0ol06P?o00<000?o0ol0N0?o0P00P@?o000b0ol00`000ol3o`0k0ol00`00 0ol3o`0J0ol00`000ol3o`1f0ol200230ol00383o`030003o`?o03`3o`030003o`?o01T3o`030003 o`?o07D3o`030003o`?o08<3o`000ol2002X0ol00383o`030003o`?o04X3o`030003 o`?o00h3o`030003o`?o04d3o`030003o`?o0:P3o`000ol00`000ol3 o`1:0ol00`000ol3o`2Z0ol00383o`030003o`?o04`3o`030003o`?o00d3o`030003o`?o04P3o`80 0:d3o`000ol00`000ol3o`0;0ol00`000ol3o`160ol00`000ol3o`2^0ol00383o`030003o`?o04h3 o`030003o`?o00/3o`030003o`?o04@3o`800;43o`007P?o0P0010?o00@000?o0ol3o`D000D3o`03 0003o`?o04l3o`030003o`?o00/3o`030003o`?o0483o`030003o`?o0;43o`007@?o00@000?o0ol0 00P3o`040003o`?o00050ol00`000ol3o`1@0ol00`000ol3o`0:0ol00`000ol3o`110ol00`000ol3 o`2b0ol001d3o`040003o`?o00090ol00`000ol3o`050ol3001A0ol00`000ol3o`090ol00`000ol3 o`100ol00`000ol3o`2c0ol001d3o`040003o`?o000:0ol00`000ol3o`040ol00`000ol3o`1B0ol0 0`000ol3o`080ol00`000ol3o`0n0ol2002f0ol001d3o`040003o`?o00070ol01@000ol3o`?o0000 1@?o00<000?o0ol0DP?o00<000?o0ol020?o00<000?o0ol0?@?o00<000?o0ol0]P?o000M0ol01000 0ol3o`001`?o00D000?o0ol3o`0000D3o`030003o`?o05<3o`030003o`?o00P3o`030003o`?o03/3 o`030003o`?o0;L3o`007P?o0P002@?o0`001P?o00<000?o0ol0E0?o00<000?o0ol01`?o00<000?o 0ol0>@?o0P00^P?o000b0ol00`000ol3o`1E0ol00`000ol3o`060ol00`000ol3o`0h0ol00`000ol3 o`2j0ol00383o`<005D3o`030003o`?o00H3o`030003o`?o03L3o`030003o`?o0;/3o`00@3o`0003o`00<3o`H006d3o`005`?o00@0 00?o0ol000/3o`030003o`?o00<3o`030003o`?o00<3o`<00>T3o`L006H3o`005`?o00@000?o0ol0 00T3o`8000H3o`030003o`?o00<3o`030003o`?o0?03o`X005`3o`005`?o00@000?o0ol000/3o`03 0003o`0000<000H3o`030003o`?o0?X3o`T005<3o`005`?o00@000?o0ol000L3o`070003o`?o0ol0 00?o00002@?o00<000?o0ol0o`?o10?o1`00C0?o000H0ol200090ol300020ol500050ol00`000ol3 o`3o0ol;0ol:00120ol00383o`030003o`?o0?l3oaD3o`d003D3o`00"], ImageRangeCache->{{{0, 394}, {240, 0}} -> {-0.822088, -0.0313136, \ 0.0161603, 0.00189786}}] }, Open ]], Cell["\<\ ClearAll[H,S,Em,A0,uY]; lambdaL=(16*A0*Em*H^3)/(3*3^(1/2)*(4*H^2 + S^2)^(3/2)); lambdaB=(2*2^(1/2)*A0*Em*S^2*(2*H^2 - S^2)^(1/2))/(4*H^2 + S^2)^(3/2); set={Em->1,A0->1,S->2}; dlambdaBH=Simplify[D[lambdaB/.set,H]]; sol=Solve[dlambdaBH==0,H]; Print[N[(lambdaL/.set)/.H->Sqrt[3]]]; Print[N[(lambdaB/.set)/.H->Sqrt[3]]]; dlambdaBH=Simplify[D[lambdaB,H]*Sqrt[S^2/4+H^2]/H /.set]; Print[\"dlambdaBH=\",dlambdaBH]; sol=Solve[dlambdaBH==0,H]; Print[\"sol=\",sol] Print[N[(lambdaL/.set)/.sol]]; Print[N[(lambdaB/.set)/.sol]]; (*pL=Plot[lambdaL/.set,{H,0,5},DisplayFunction->Identity, AxesLabel->{\"H\",\"lambda at CPs\"}, PlotLabel->\"lambdaCP vs H (S=2,E=A0=1)\"]; pB=Plot[lambdaB/.set,{H,Sqrt[2.00001],5},DisplayFunction->Identity]; Show[{pL,pB},DisplayFunction->$DisplayFunction];*)\ \>", "Input", CellFrame->True, AspectRatioFixed->True, Background->GrayLevel[0.833326]], Cell["\<\ ClearAll[S,H,A0,Em,V0]; Em=1; S=2; A0=1; L0=V0/(2*A0); H=Sqrt[L0^2-S^2/4]; lambdaL=(16*A0*Em*H^3)/(3*3^(1/2)*(4*H^2 + S^2)^(3/2)); lambdaL=Simplify[lambdaL]; Print[\"L0=\",L0];Print[\"H=\",H]; Print[lambdaL];Plot[lambdaL/V0,{V0,2,10}]; lambdaB=(2*2^(1/2)*A0*Em*S^2*(2*H^2 - S^2)^(1/2))/(4*H^2 + S^2)^(3/2); lambdaB=Simplify[lambdaB]; Print[lambdaB];Plot[lambdaB/V0,{V0,2,10}]; Plot[lambdaB,{V0,2,10}];\ \>", "Input", CellFrame->True, AspectRatioFixed->True, Background->GrayLevel[0.833326]], Cell["\<\ K=(EA/L)*{{0,0},{0,1}}+m*g/L*{{1,0},{0,1}}; Ms={{1/Sqrt[m],0},{0,1/Sqrt[m]}}; eig=Simplify[Eigenvalues[Ms.K.Ms]]; Print[\"squared frequencies=\",eig]; \ \>", "Input", CellFrame->True, AspectRatioFixed->True, Background->GrayLevel[0.833326]], Cell[TextData[{ StyleBox["Em=A0=1; S=2; s=Sin[a]; c=Cos[a]; H=(Sin[a]/Cos[a])*S/2;\n\ lambdaL=(16*A0*Em*H^3)/(3*3^(1/2)*(4*H^2 + S^2)^(3/2));\n\ Print[\"lambdaL=\",Simplify[lambdaL]//InputForm];\n", AspectRatioFixed->True], StyleBox["lambdaL=Simplify[(2*Tan[a]^3)/(3*3^(1/2)*(Sec[a]^3))];\n", AspectRatioFixed->True, FontWeight->"Plain"], StyleBox["Print[\"lambdaL=\",lambdaL];", AspectRatioFixed->True] }], "Input", CellFrame->True, AspectRatioFixed->True, Background->GrayLevel[0.833326]], Cell["\<\ AssembleMasterStiffOfQuiz3Arch[S_,H_,Em_,A0_,uX_,uY_]:= Module[{K1,K2,K3,K4,K5,K6,K,S,H}, K=Table[0,{2},{2}]; S=S; H=H; K1=FormTangentStiff2DTwoNodeBar[{-S/2,0},{0,3*H/4}, {0,0},{uX,uY},Em,4*A0,0]; (*Print[K1//TableForm];*) K=MergeElemIntoMasterStiff[K1,{0,0,1,2},K]; K2=FormTangentStiff2DTwoNodeBar[{-S/2,S},{-S/2,H}, {uX,uY},{0,0},Em,4*A0,0]; (*Print[K2//TableForm];*) K=MergeElemIntoMasterStiff[K2,{1,2,0,0},K]; K3=FormTangentStiff2DTwoNodeBar[{0,H},{S/2,3*H/4}, {0,0},{uX,uY},Em,4*A0,0]; (*Print[K1//TableForm];*) K=MergeElemIntoMasterStiff[K3,{0,0,1,2},K]; K4=FormTangentStiff2DTwoNodeBar[{S/2,3*H/4},{S,0}, {uX,uY},{0,0},Em,4*A0,0]; (*Print[K2//TableForm];*) K=MergeElemIntoMasterStiff[K4,{1,2,0,0},K]; K5=FormTangentStiff2DTwoNodeBar[{-S/2,0},{0,H}, {0,0},{uX,uY},Em,A0,0]; (*Print[K1//TableForm];*) K=MergeElemIntoMasterStiff[K5,{0,0,1,2},K]; K6=FormTangentStiff2DTwoNodeBar[{0,H},{S/2,0}, {uX,uY},{0,0},Em,A0,0]; (*Print[K2//TableForm];*) K=MergeElemIntoMasterStiff[K6,{1,2,0,0},K]; Return[Simplify[K]] ];\ \>", "Input", CellFrame->True, AspectRatioFixed->True, Background->GrayLevel[0.833326]], Cell[CellGroupData[{ Cell["\<\ ClearAll[V0,H,S,Em,A0]; (*S=2; Em=1; A0=1;*) H=Sqrt[V0/(2*A0)-S^2/4]; (*L0=Sqrt[S^2/4+H^2]; A0=V0/(2*L0);*) lambdaL=(16*A0*Em*H^3)/(3*3^(1/2)*(4*H^2 + S^2)^(3/2))/V0; lambdaB=(2*2^(1/2)*A0*Em*S^2*(2*H^2 - S^2)^(1/2))/(4*H^2 + S^2)^(3/2)/V0; (*Plot[{lambdaL,lambdaB},{H,0,3},PlotRange->All];*) Plot[{lambdaL,lambdaB},{V0,0,10},PlotRange->All]; (*Plot[{lambdaL,lambdaB},{H,0,3},PlotRange->All];*) V0best=5;Hbest=H/.V0->V0best; Print[{Hbest,N[ArcTan[Hbest]*(180/Pi)]}]; Print[Simplify[D[lambdaL,V0]]]; Print[Simplify[lambdaL]]; \ \>", "Input", CellMargins->{{16, Inherited}, {Inherited, Inherited}}, AspectRatioFixed->True, Background->RGBColor[0, 1, 0]], 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.0238095 0.952381 0.0147151 0.588604 [ [.21429 .00222 -9 -9 ] [.21429 .00222 9 0 ] [.40476 .00222 -9 -9 ] [.40476 .00222 9 0 ] [.59524 .00222 -9 -9 ] [.59524 .00222 9 0 ] [.78571 .00222 -9 -9 ] [.78571 .00222 9 0 ] [.97619 .00222 -3 -9 ] [.97619 .00222 3 0 ] [.01131 .13244 -18 -4.5 ] [.01131 .13244 0 4.5 ] [.01131 .25016 -18 -4.5 ] [.01131 .25016 0 4.5 ] [.01131 .36788 -18 -4.5 ] [.01131 .36788 0 4.5 ] [.01131 .4856 -18 -4.5 ] [.01131 .4856 0 4.5 ] [.01131 .60332 -6 -4.5 ] [.01131 .60332 0 4.5 ] [ 0 0 0 0 ] [ 1 .61803 0 0 ] ] MathScale % Start of Graphics 1 setlinecap 1 setlinejoin newpath 0 g .25 Mabswid [ ] 0 setdash .21429 .01472 m .21429 .02097 L s [(0.2)] .21429 .00222 0 1 Mshowa .40476 .01472 m .40476 .02097 L s [(0.4)] .40476 .00222 0 1 Mshowa .59524 .01472 m .59524 .02097 L s [(0.6)] .59524 .00222 0 1 Mshowa .78571 .01472 m .78571 .02097 L s [(0.8)] .78571 .00222 0 1 Mshowa .97619 .01472 m .97619 .02097 L s [(1)] .97619 .00222 0 1 Mshowa .125 Mabswid .07143 .01472 m .07143 .01847 L s .11905 .01472 m .11905 .01847 L s .16667 .01472 m .16667 .01847 L s .2619 .01472 m .2619 .01847 L s .30952 .01472 m .30952 .01847 L s .35714 .01472 m .35714 .01847 L s .45238 .01472 m .45238 .01847 L s .5 .01472 m .5 .01847 L s .54762 .01472 m .54762 .01847 L s .64286 .01472 m .64286 .01847 L s .69048 .01472 m .69048 .01847 L s .7381 .01472 m .7381 .01847 L s .83333 .01472 m .83333 .01847 L s .88095 .01472 m .88095 .01847 L s .92857 .01472 m .92857 .01847 L s .25 Mabswid 0 .01472 m 1 .01472 L s .02381 .13244 m .03006 .13244 L s [(0.2)] .01131 .13244 1 0 Mshowa .02381 .25016 m .03006 .25016 L s [(0.4)] .01131 .25016 1 0 Mshowa .02381 .36788 m .03006 .36788 L s [(0.6)] .01131 .36788 1 0 Mshowa .02381 .4856 m .03006 .4856 L s [(0.8)] .01131 .4856 1 0 Mshowa .02381 .60332 m .03006 .60332 L s [(1)] .01131 .60332 1 0 Mshowa .125 Mabswid .02381 .04415 m .02756 .04415 L s .02381 .07358 m .02756 .07358 L s .02381 .10301 m .02756 .10301 L s .02381 .16187 m .02756 .16187 L s .02381 .1913 m .02756 .1913 L s .02381 .22073 m .02756 .22073 L s .02381 .27959 m .02756 .27959 L s .02381 .30902 m .02756 .30902 L s .02381 .33845 m .02756 .33845 L s .02381 .39731 m .02756 .39731 L s .02381 .42674 m .02756 .42674 L s .02381 .45617 m .02756 .45617 L s .02381 .51503 m .02756 .51503 L s .02381 .54446 m .02756 .54446 L s .02381 .57389 m .02756 .57389 L s .25 Mabswid .02381 0 m .02381 .61803 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->{525, 323.875}, ImageMargins->{{36, 0}, {0, 0}}, ImageRegion->{{0, 1}, {0, 1}}, Background->RGBColor[0, 1, 1], ImageCache->GraphicsData["Bitmap", "\<\ CF5dJ6E]HGAYHf4PAg9QL6QYHg003oo`00LP00ool0100000000?oo 003oo`000003003oo`8000001000ool00`0000000?oo003oo`1A003oo`0400000000ool00?oo0000 00<00?oo0P000005003oo`0300000000ool00?oo05000?oo00@00000003oo`00ool000000`00ool2 000000<00?oo00@00000003oo`00ool00000D@00ool0100000000?oo003oo`000003003oo`800000 0`00ool0100000000?oo003oo`00001H003oo`0300000000ool00?oo00d00?oo001b003oo`040000 0000ool00?oo000000X00?oo00<00000003oo`00ool0D000ool0100000000?oo003oo`000008003o o`@00000D@00ool0100000000?oo003oo`000008003oo`0400000000ool00?oo000005400?oo00@0 0000003oo`00ool000002000ool0100000000?oo003oo`00001H003oo`0300000000ool00?oo00d0 0?oo001b003oo`0400000000ool00?oo000000/00?oo00<00000003oo`00ool0C`00ool010000000 0?oo003oo`000008003oo`0300000000ool0000005800?oo00@00000003oo`00ool000002000ool3 000005800?oo00@00000003oo`00ool000002@00ool2000005T00?oo00<00000003oo`00ool03@00 ool007800?oo00@00000003oo`00ool000002000ool0100000000?oo003oo`00001A003oo`040000 0000ool00?oo000000P00?oo00<00000003oo`000000DP00ool0100000000?oo003oo`000008003o o`0300000000ool00?oo05800?oo00@00000003oo`00ool000002000ool0100000000?oo003oo`00 001H003oo`0300000000ool00?oo00d00?oo000I003oo`0300000000ool00?oo05L00?oo0P00000: 003oo`800000D`00ool2000000X00?oo0P00001C003oo`8000002P00ool3000005800?oo0P00000: 003oo`800000F000ool2000000l00?oo000I003oo`0300000000ool00?oo0?l00?oolP00ool001T0 0?oo00<00000003oo`00ool0o`00ooob003oo`006@00ool00`0000000?oo003oo`3o003ooo800?oo 000I003oo`0300000000ool00?oo0?l00?oolP00ool001T00?oo00<00000003oo`00ool0o`00ooob 003oo`006@00ool00`0000000?oo003oo`3o003ooo800?oo000I003oo`0300000000ool00?oo0?l0 0?oolP00ool000d00?ooo`00003o000000800?oo000I003oo`0300000000ool00?oo01D00?oo00<0 0000003oo`00ool05@00ool00`0000000?oo003oo`0F003oo`0300000000ool00?oo01D00?oo00<0 0000003oo`00ool05@00ool00`0000000?oo003oo`0F003oo`0300000000ool00?oo01D00?oo00<0 0000003oo`00ool05@00ool00`0000000?oo003oo`0E003oo`0300000000ool00?oo01H00?oo00<0 0000003oo`00ool05@00ool00`0000000?oo003oo`0E003oo`0300000000ool00?oo01D00?oo00<0 0000003oo`00ool05P00ool00`0000000?oo003oo`0E003oo`0300000000ool00?oo01D00?oo00<0 0000003oo`00ool05P00ool00`0000000?oo003oo`0E003oo`0300000000ool00?oo01D00?oo00<0 0000003oo`00ool05@00ool00`0000000?oo003oo`0<003oo`006@00ool00`0000000?oo003oo`0E 003oo`0300000000ool00?oo01D00?oo00<00000003oo`00ool05P00ool00`0000000?oo003oo`0E 003oo`0300000000ool00?oo01D00?oo00<00000003oo`00ool05P00ool00`0000000?oo003oo`0E 003oo`0300000000ool00?oo01D00?oo00<00000003oo`00ool05@00ool00`0000000?oo003oo`0F 003oo`0300000000ool00?oo01D00?oo00<00000003oo`00ool05@00ool00`0000000?oo003oo`0E 003oo`0300000000ool00?oo01H00?oo00<00000003oo`00ool05@00ool00`0000000?oo003oo`0E 003oo`0300000000ool00?oo01H00?oo00<00000003oo`00ool05@00ool00`0000000?oo003oo`0E 003oo`0300000000ool00?oo01D00?oo00<00000003oo`00ool03000ool001T00?oo00<00000003o o`00ool0GP00ool00`0000000?oo003oo`1N003oo`0300000000ool00?oo05h00?oo00<00000003o o`00ool0GP00ool00`0000000?oo003oo`1N003oo`0300000000ool00?oo00`00?oo000I003oo`03 00000000ool00?oo0?l00?oolP00ool001T00?oo00<00000003oo`00ool0o`00ooob003oo`006@00 ool00`0000000?oo003oo`3o003ooo800?oo000I003oo`0300000000ool00?oo0?l00?oolP00ool0 01T00?oo00<00000003oo`00ool0o`00ooob003oo`006@00ool00`0000000?oo003oo`3o003ooo80 0?oo000I003oo`0300000000ool00?oo0?l00?oolP00ool001T00?oo00<00000003oo`00ool0o`00 ooob003oo`006@00ool00`0000000?oo003oo`3o003ooo800?oo000I003oo`0300000000ool00?oo 0?l00?oolP00ool001T00?oo00<00000003oo`00ool0o`00ooob003oo`006@00ool300000?l00?oo lP00ool001T00?oo00<00000003oo`00ool0o`00ooob003oo`006@00ool00`0000000?oo003oo`3o 003ooo800?oo000I003oo`0300000000ool00?oo0?l00?oolP00ool001T00?oo00<00000003oo`00 ool0o`00ooob003oo`006@00ool00`0000000?oo003oo`3o003ooo800?oo000I003oo`0300000000 ool00?oo0?l00?oolP00ool001T00?oo00<00000003oo`00ool0o`00ooob003oo`006@00ool00`00 00000?oo003oo`3o003ooo800?oo000I003oo`0300000000ool00?oo0?l00?oolP00ool001T00?oo 00<00000003oo`00ool0o`00ooob003oo`006@00ool00`0000000?oo003oo`3o003ooo800?oo000I 003oo`0300000000ool00?oo0?l00?oolP00ool001T00?oo00<00000003oo`00ool0o`00ooob003o o`006@00ool00`0000000?oo003oo`3o003ooo800?oo000I003oo`<00000o`00ooob003oo`006@00 ool00`0000000?oo003oo`3o003ooo800?oo000I003oo`0300000000ool00?oo0?l00?oolP00ool0 01T00?oo00<00000003oo`00ool0o`00ooob003oo`006@00ool00`0000000?oo003oo`3o003ooo80 0?oo000I003oo`0300000000ool00?oo0?l00?oolP00ool001T00?oo00<00000003oo`00ool0o`00 ooob003oo`006@00ool00`0000000?oo003oo`3o003ooo800?oo000I003oo`0300000000ool00?oo 0?l00?oolP00ool001T00?oo00<00000003oo`00ool0o`00ooob003oo`006@00ool00`0000000?oo 003oo`3o003ooo800?oo000I003oo`0300000000ool00?oo0?l00?oolP00ool001T00?oo00<00000 003oo`00ool0o`00ooob003oo`006@00ool00`0000000?oo003oo`3o003ooo800?oo000I003oo`03 00000000ool00?oo0?l00?oolP00ool001T00?oo0`00003o003ooo800?oo000I003oo`0300000000 ool00?oo0?l00?oolP00ool001T00?oo00<00000003oo`00ool0o`00ooob003oo`006@00ool00`00 00000?oo003oo`3o003ooo800?oo000I003oo`0300000000ool00?oo0?l00?oolP00ool001T00?oo 00<00000003oo`00ool0o`00ooob003oo`006@00ool00`0000000?oo003oo`3o003ooo800?oo000I 003oo`0300000000ool00?oo0?l00?oolP00ool001T00?oo00<00000003oo`00ool0o`00ooob003o o`006@00ool00`0000000?oo003oo`3o003ooo800?oo000I003oo`0300000000ool00?oo0?l00?oo lP00ool001T00?oo00<00000003oo`00ool0o`00ooob003oo`006@00ool00`0000000?oo003oo`3o 003ooo800?oo0003003oo`8000001000ool2000000<00?oo10000007003oo`0300000000ool00?oo 0?l00?oolP00ool000800?oo00@00000003oo`00ool000000`00ool2000000@00?oo00<00000003o o`00ool01`00ool00`0000000?oo003oo`3o003ooo800?oo0002003oo`0400000000ool00?oo0000 00X00?oo00<00000003oo`00ool01P00ool400000?l00?ool@00ool000800?oo00@00000003oo`00 ool000002`00ool00`0000000?oo003oo`05003oo`0300000000ool00?oo0?l00?oolP00ool00080 0?oo00@00000003oo`00ool000002000ool0100000000?oo003oo`000007003oo`0300000000ool0 0?oo0?l00?oolP00ool000<00?oo0P00000:003oo`8000002000ool00`0000000?oo003oo`3o003o oo800?oo000I003oo`0300000000ool00?oo0?l00?oolP00ool001T00?oo00<00000003oo`00ool0 o`00ooob003oo`006@00ool00`0000000?oo003oo`3o003ooo800?oo000I003oo`0300000000ool0 0?oo0?l00?oolP00ool001T00?oo00<00000003oo`00ool0o`00ooob003oo`006@00ool00`000000 0?oo003oo`3o003ooo800?oo000I003oo`0300000000ool00?oo0?l00?oolP00ool001T00?oo00<0 0000003oo`00ool0o`00ooob003oo`006@00ool00`0000000?oo003oo`3o003ooo800?oo000I003o o`0300000000ool00?oo0?l00?oolP00ool001T00?oo00<00000003oo`00ool0o`00ooob003oo`00 6@00ool300000?l00?oolP00ool001T00?oo00<00000003oo`00ool0o`00ooob003oo`006@00ool0 0`0000000?oo003oo`3o003ooo800?oo000I003oo`0300000000ool00?oo0?l00?oolP00ool001T0 0?oo00<00000003oo`00ool0o`00ooob003oo`006@00ool00`0000000?oo003oo`3o003ooo800?oo 000I003oo`0300000000ool00?oo0?l00?oolP00ool001T00?oo00<00000003oo`00ool0o`00ooob 003oo`006@00ool00`0000000?oo003oo`3o003ooo800?oo000I003oo`0300000000ool00?oo0?l0 0?oolP00ool001T00?oo00<00000003oo`00ool0o`00ooob003oo`006@00ool00`0000000?oo003o o`3o003ooo800?oo000I003oo`0300000000ool00?oo0?l00?oolP00ool001T00?oo00<00000003o o`00ool0o`00ooob003oo`006@00ool00`0000000?oo003oo`3o003ooo800?oo000I003oo`<00000 o`00ooob003oo`006@00ool00`0000000?oo003oo`3o003ooo800?oo000I003oo`0300000000ool0 0?oo0?l00?oolP00ool001T00?oo00<00000003oo`00ool0o`00ooob003oo`006@00ool00`000000 0?oo003oo`3o003ooo800?oo000I003oo`0300000000ool00?oo0?l00?oolP00ool001T00?oo00<0 0000003oo`00ool0o`00ooob003oo`006@00ool00`0000000?oo003oo`3o003ooo800?oo000I003o o`0300000000ool00?oo0?l00?oolP00ool001T00?oo00<00000003oo`00ool0o`00ooob003oo`00 6@00ool00`0000000?oo003oo`3o003ooo800?oo000I003oo`0300000000ool00?oo0?l00?oolP00 ool001T00?oo00<00000003oo`00ool0o`00ooob003oo`006@00ool00`0000000?oo003oo`3o003o oo800?oo000I003oo`0300000000ool00?oo0?l00?oolP00ool001T00?oo0`00003o003ooo800?oo 000I003oo`0300000000ool00?oo0?l00?oolP00ool001T00?oo00<00000003oo`00ool0o`00ooob 003oo`006@00ool00`0000000?oo003oo`3o003ooo800?oo000I003oo`0300000000ool00?oo0?l0 0?oolP00ool001T00?oo00<00000003oo`00ool0o`00ooob003oo`006@00ool00`0000000?oo003o o`3o003ooo800?oo000I003oo`0300000000ool00?oo0?l00?oolP00ool001T00?oo00<00000003o o`00ool0o`00ooob003oo`006@00ool00`0000000?oo003oo`3o003ooo800?oo000I003oo`030000 0000ool00?oo0?l00?oolP00ool001T00?oo00<00000003oo`00ool0o`00ooob003oo`006@00ool0 0`0000000?oo003oo`3o003ooo800?oo0003003oo`8000001000ool2000000@00?oo0`000007003o o`0300000000ool00?oo0?l00?oolP00ool000800?oo00@00000003oo`00ool000000`00ool20000 00D00?oo00<00000003oo`00ool01P00ool00`0000000?oo003oo`3o003ooo800?oo0002003oo`04 00000000ool00?oo000000P00?oo10000007003oo`@00000o`00oooa003oo`000P00ool010000000 0?oo003oo`000008003oo`0300000000ool0000000P00?oo00<00000003oo`00ool0o`00ooob003o o`000P00ool0100000000?oo003oo`000008003oo`0300000000ool0000000P00?oo00<00000003o o`00ool0o`00ooob003oo`000`00ool2000000X00?oo0P000008003oo`0300000000ool00?oo0?l0 0?oolP00ool001T00?oo00<00000003oo`00ool0o`00ooob003oo`006@00ool00`0000000?oo003o o`3o003ooo800?oo000I003oo`0300000000ool00?oo0?l00?oolP00ool001T00?oo00<00000003o o`00ool0o`00ooob003oo`006@00ool00`0000000?oo003oo`3o003ooo800?oo000I003oo`030000 0000ool00?oo0?l00?oolP00ool001T00?oo00<00000003oo`00ool0o`00ooob003oo`006@00ool0 0`0000000?oo003oo`3o003ooo800?oo000I003oo`0300000000ool00?oo0?l00?oolP00ool001T0 0?oo00<00000003oo`00ool0o`00ooob003oo`006@00ool00`0000000?oo003oo`3o003ooo800?oo 000I003oo`<00000o`00ooob003oo`006@00ool00`0000000?oo003oo`3o003ooo800?oo000I003o o`0300000000ool00?oo0?l00?oolP00ool001T00?oo00<00000003oo`00ool0o`00ooob003oo`00 6@00ool00`0000000?oo003oo`3o003ooo800?oo000I003oo`0300000000ool00?oo0?l00?oolP00 ool001T00?oo00<00000003oo`00ool0o`00ooob003oo`006@00ool00`0000000?oo003oo`3o003o oo800?oo000I003oo`0300000000ool00?oo0?l00?oolP00ool001T00?oo00<00000003oo`00ool0 o`00ooob003oo`006@00ool00`0000000?oo003oo`3o003ooo800?oo000I003oo`0300000000ool0 0?oo0?l00?oolP00ool001T00?oo00<00000003oo`00ool0o`00ooob003oo`006@00ool00`000000 0?oo003oo`3o003ooo800?oo000I003oo`0300000000ool00?oo0?l00?oolP00ool001T00?oo0`00 003o003ooo800?oo000I003oo`0300000000ool00?oo0?l00?oolP00ool001T00?oo00<00000003o o`00ool0o`00ooob003oo`006@00ool00`0000000?oo003oo`3o003ooo800?oo000I003oo`030000 0000ool00?oo0?l00?oolP00ool001T00?oo00<00000003oo`00ool0o`00ooob003oo`006@00ool0 0`0000000?oo003oo`3o003ooo800?oo000I003oo`0300000000ool00?oo0?l00?oolP00ool001T0 0?oo00<00000003oo`00ool0o`00ooob003oo`006@00ool00`0000000?oo003oo`3o003ooo800?oo 000I003oo`0300000000ool00?oo0?l00?oolP00ool001T00?oo00<00000003oo`00ool0o`00ooob 003oo`006@00ool00`0000000?oo003oo`3o003ooo800?oo000I003oo`0300000000ool00?oo0?l0 0?oolP00ool001T00?oo00<00000003oo`00ool0o`00ooob003oo`006@00ool300000?l00?oolP00 ool001T00?oo00<00000003oo`00ool0o`00ooob003oo`006@00ool00`0000000?oo003oo`3o003o oo800?oo000I003oo`0300000000ool00?oo0?l00?oolP00ool001T00?oo00<00000003oo`00ool0 o`00ooob003oo`006@00ool00`0000000?oo003oo`3o003ooo800?oo000I003oo`0300000000ool0 0?oo0?l00?oolP00ool001T00?oo00<00000003oo`00ool0o`00ooob003oo`006@00ool00`000000 0?oo003oo`3o003ooo800?oo000I003oo`0300000000ool00?oo0?l00?oolP00ool001T00?oo00<0 0000003oo`00ool0o`00ooob003oo`006@00ool00`0000000?oo003oo`3o003ooo800?oo000I003o o`0300000000ool00?oo0?l00?oolP00ool000<00?oo0P000004003oo`8000001000ool2000000P0 0?oo00<00000003oo`00ool0o`00ooob003oo`000P00ool0100000000?oo003oo`000003003oo`80 00000`00ool0100000000?oo003oo`000007003oo`0300000000ool00?oo0?l00?oolP00ool00080 0?oo00@00000003oo`00ool000002000ool0100000000?oo003oo`000007003oo`@00000o`00oooa 003oo`000P00ool0100000000?oo003oo`000008003oo`<000002000ool00`0000000?oo003oo`3o 003ooo800?oo0002003oo`0400000000ool00?oo000000P00?oo00<00000003oo`00ool02000ool0 0`0000000?oo003oo`3o003ooo800?oo0003003oo`8000002P00ool3000000L00?oo00<00000003o o`00ool0o`00ooob003oo`006@00ool00`0000000?oo003oo`3o003ooo800?oo000I003oo`030000 0000ool00?oo0?l00?oolP00ool001T00?oo00<00000003oo`00ool0o`00ooob003oo`006@00ool0 0`0000000?oo003oo`3o003ooo800?oo000I003oo`0300000000ool00?oo0?l00?oolP00ool001T0 0?oo00<00000003oo`00ool0o`00ooob003oo`006@00ool00`0000000?oo003oo`3o003ooo800?oo 000I003oo`0300000000ool00?oo0?l00?oolP00ool001T00?oo00<00000003oo`00ool0o`00ooob 003oo`006@00ool00`0000000?oo003oo`3o003ooo800?oo000I003oo`0300000000ool00?oo0?l0 0?oolP00ool001T00?oo0`00003o003ooo800?oo000I003oo`0300000000ool00?oo0?l00?oolP00 ool001T00?oo00<00000003oo`00ool0o`00ooob003oo`006@00ool00`0000000?oo003oo`3o003o oo800?oo000I003oo`0300000000ool00?oo0?l00?oolP00ool001T00?oo00<00000003oo`00ool0 o`00ooob003oo`006@00ool00`0000000?oo003oo`3o003ooo800?oo000I003oo`0300000000ool0 0?oo0?l00?oolP00ool001T00?oo00<00000003oo`00ool0o`00ooob003oo`006@00ool00`000000 0?oo003oo`3o003ooo800?oo000I003oo`0300000000ool00?oo0?l00?oolP00ool001T00?oo00<0 0000003oo`00ool0o`00ooob003oo`006@00ool00`0000000?oo003oo`3o003ooo800?oo000I003o o`0300000000ool00?oo0?l00?oolP00ool001T00?oo00<00000003oo`00ool0o`00ooob003oo`00 6@00ool300000?l00?oolP00ool001T00?oo00<00000003oo`00ool0o`00ooob003oo`006@00ool0 0`0000000?oo003oo`3o003ooo800?oo000I003oo`0300000000ool00?oo0?l00?oolP00ool001T0 0?oo00<00000003oo`00ool0o`00ooob003oo`006@00ool00`0000000?oo003oo`3o003ooo800?oo 000I003oo`0300000000ool00?oo0?l00?oolP00ool001T00?oo00<00000003oo`00ool0o`00ooob 003oo`006@00ool00`0000000?oo003oo`3o003ooo800?oo000I003oo`0300000000ool00?oo0?l0 0?oolP00ool001T00?oo00<00000003oo`00ool0o`00ooob003oo`006@00ool00`0000000?oo003o o`3o003ooo800?oo000I003oo`0300000000ool00?oo0?l00?oolP00ool001T00?oo00<00000003o o`00ool0o`00ooob003oo`006@00ool00`0000000?oo003oo`3o003ooo800?oo000I003oo`<00000 o`00ooob003oo`006@00ool00`0000000?oo003oo`3o003ooo800?oo000I003oo`0300000000ool0 0?oo0?l00?oolP00ool001T00?oo00<00000003oo`00ool0o`00ooob003oo`006@00ool00`000000 0?oo003oo`3o003ooo800?oo000I003oo`0300000000ool00?oo0?l00?oolP00ool001T00?oo00<0 0000003oo`00ool0o`00ooob003oo`006@00ool00`0000000?oo003oo`3o003ooo800?oo000I003o o`0300000000ool00?oo0?l00?oolP00ool001T00?oo00<00000003oo`00ool0o`00ooob003oo`00 6@00ool00`0000000?oo003oo`3o003ooo800?oo000I003oo`0300000000ool00?oo0?l00?oolP00 ool001T00?oo00<00000003oo`00ool0o`00ooob003oo`000`00ool2000000@00?oo0P000004003o o`8000002000ool00`0000000?oo003oo`3o003ooo800?oo0002003oo`0400000000ool00?oo0000 00<00?oo0P000003003oo`0400000000ool00?oo000000L00?oo00<00000003oo`00ool0o`00ooob 003oo`000P00ool0100000000?oo003oo`000008003oo`0400000000ool00?oo000000L00?oo1000 003o003ooo400?oo0002003oo`0400000000ool00?oo000000T00?oo0P000008003oo`0300000000 ool00?oo0?l00?oolP00ool000800?oo00@00000003oo`00ool000002000ool0100000000?oo003o o`000007003oo`0300000000ool00?oo0?l00?oolP00ool000<00?oo0P00000:003oo`8000002000 ool00`0000000?oo003oo`3o003ooo800?oo000I003oo`0300000000ool00?oo0?l00?oolP00ool0 01T00?oo00<00000003oo`00ool0o`00ooob003oo`006@00ool00`0000000?oo003oo`3o003ooo80 0?oo000I003oo`0300000000ool00?oo0?l00?oolP00ool001T00?oo00<00000003oo`00ool0o`00 ooob003oo`006@00ool00`0000000?oo003oo`3o003ooo800?oo000I003oo`0300000000ool00?oo 0?l00?oolP00ool001T00?oo00<00000003oo`00ool0o`00ooob003oo`006@00ool00`0000000?oo 003oo`3o003ooo800?oo000I003oo`0300000000ool00?oo0?l00?oolP00ool001T00?oo00<00000 003oo`00ool0o`00ooob003oo`006@00ool300000?l00?oolP00ool001T00?oo00<00000003oo`00 ool0o`00ooob003oo`006@00ool00`0000000?oo003oo`3o003ooo800?oo000I003oo`0300000000 ool00?oo0?l00?oolP00ool001T00?oo00<00000003oo`00ool0o`00ooob003oo`006@00ool00`00 00000?oo003oo`3o003ooo800?oo000I003oo`0300000000ool00?oo0?l00?oolP00ool001T00?oo 00<00000003oo`00ool0o`00ooob003oo`006@00ool00`0000000?oo003oo`3o003ooo800?oo000I 003oo`0300000000ool00?oo0?l00?oolP00ool001T00?oo00<00000003oo`00ool0o`00ooob003o o`006@00ool00`0000000?oo003oo`3o003ooo800?oo000I003oo`0300000000ool00?oo0?l00?oo lP00ool001T00?oo00<00000003oo`00ool0o`00ooob003oo`006@00ool00`0000000?oo003oo`3o 003ooo800?oo000I003oo`<00000o`00ooob003oo`006@00ool00`0000000?oo003oo`3o003ooo80 0?oo000I003oo`0300000000ool00?oo0?l00?oolP00ool001T00?oo00<00000003oo`00ool0o`00 ooob003oo`006@00ool00`0000000?oo003oo`3o003ooo800?oo000I003oo`0300000000ool00?oo 0?l00?oolP00ool001T00?oo00<00000003oo`00ool0o`00ooob003oo`006@00ool00`0000000?oo 003oo`3o003ooo800?oo000I003oo`0300000000ool00?oo0?l00?oolP00ool001T00?oo00<00000 003oo`00ool0o`00ooob003oo`006@00ool00`0000000?oo003oo`3o003ooo800?oo000I003oo`03 00000000ool00?oo0?l00?oolP00ool001T00?oo00<00000003oo`00ool0o`00ooob003oo`006@00 ool00`0000000?oo003oo`3o003ooo800?oo000I003oo`0300000000ool00?oo0?l00?oolP00ool0 01T00?oo0`00003o003ooo800?oo000I003oo`0300000000ool00?oo0?l00?oolP00ool001T00?oo 00<00000003oo`00ool0o`00ooob003oo`006@00ool00`0000000?oo003oo`3o003ooo800?oo000I 003oo`0300000000ool00?oo0?l00?oolP00ool001T00?oo00<00000003oo`00ool0o`00ooob003o o`006@00ool00`0000000?oo003oo`3o003ooo800?oo000I003oo`0300000000ool00?oo0?l00?oo lP00ool001T00?oo00<00000003oo`00ool0o`00ooob003oo`006@00ool00`0000000?oo003oo`3o 003ooo800?oo000I003oo`0300000000ool00?oo0?l00?oolP00ool001T00?oo00<00000003oo`00 ool0o`00ooob003oo`006@00ool00`0000000?oo003oo`3o003ooo800?oo000>003oo`<000002000 ool00`0000000?oo003oo`3o003ooo800?oo000?003oo`0300000000ool00?oo00L00?oo00<00000 003oo`00ool0o`00ooob003oo`003`00ool00`0000000?oo003oo`07003oo`@00000o`00oooa003o o`003`00ool00`0000000?oo003oo`07003oo`0300000000ool00?oo0?l00?oolP00ool000l00?oo 00<00000003oo`00ool01`00ool00`0000000?oo003oo`3o003ooo800?oo000>003oo`8000002@00 ool00`0000000?oo003oo`3o003ooo800?oo000I003oo`0300000000ool00?oo0?l00?oolP00ool0 01T00?oo00<00000003oo`00ool0o`00ooob003oo`006@00ool00`0000000?oo003oo`3o003ooo80 0?oo000I003oo`0300000000ool00?oo0?l00?oolP00ool00001\ \>"], ImageRangeCache->{{{0, 524}, {322.875, 0}} -> {-0.0523433, -0.0512331, 0.00206008, 0.00333328}}] }, Open ]] }, FrontEndVersion->"4.2 for Macintosh", ScreenRectangle->{{0, 1920}, {0, 1180}}, WindowToolbars->{}, CellGrouping->Manual, WindowSize->{827, 717}, WindowMargins->{{20, Automatic}, {Automatic, 3}}, 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[1754, 51, 1286, 21, 358, "Input"], Cell[3043, 74, 1370, 22, 433, "Input"], Cell[4416, 98, 328, 5, 102, "Input"], Cell[4747, 105, 108, 2, 57, "Input"], Cell[4858, 109, 1370, 23, 478, "Input"], Cell[CellGroupData[{ Cell[6253, 136, 452, 8, 118, "Input"], Cell[6708, 146, 14935, 554, 317, 3798, 411, "GraphicsData", "PostScript", \ "Graphics", Evaluatable->False] }, Open ]], Cell[21658, 703, 2448, 37, 658, "Input"], Cell[24109, 742, 570, 10, 148, "Input"], Cell[CellGroupData[{ Cell[24704, 756, 692, 13, 163, "Input"], Cell[25399, 771, 22736, 777, 265, 6085, 565, "GraphicsData", "PostScript", \ "Graphics", Evaluatable->False] }, Open ]], Cell[48150, 1551, 888, 23, 313, "Input"], Cell[49041, 1576, 507, 13, 163, "Input"], Cell[49551, 1591, 257, 10, 118, "Input"], Cell[49811, 1603, 518, 13, 102, "Input"], Cell[50332, 1618, 1264, 28, 388, "Input"], Cell[CellGroupData[{ Cell[51621, 1650, 668, 17, 223, "Input"], Cell[52292, 1669, 20659, 425, 348, 2862, 199, "GraphicsData", "PostScript", \ "Graphics", Evaluatable->False] }, Open ]] } ] *) (******************************************************************* End of Mathematica Notebook file. *******************************************************************)