Re: 2nd attempt at post with corrected code
- To: mathgroup at smc.vnet.net
- Subject: [mg69304] Re: 2nd attempt at post with corrected code
- From: Jean-Marc Gulliet <jeanmarc.gulliet at gmail.com>
- Date: Wed, 6 Sep 2006 04:28:16 -0400 (EDT)
- Organization: The Open University, Milton Keynes, UK
- References: <edjg8p$lq8$1@smc.vnet.net>
- Sender: owner-wri-mathgroup at wolfram.com
T Harris wrote: > Hello, > > I have tried everything I know to get this code to look as it does in the > notebook. I can't do it. However, once it is copied and pasted in a > notebook, it works on my machine and looks like it originally did. If it > doesn't work for anyone else, thanks for the try. I can't remove my last > message. This is my second attempt at this. I am sorry others have tried > and been unable to help because of missing code. If it doesn't work now, > thanks, and I won't try again to repost this. There should be three parts > to the code. If there is any type of error, please forget about it. > Thanks. > > The two questions I had were (1) How do I find the planar area of the > ellipsoid? and (2) How do I identify what plane it plots out on? > > In[1222]:= > \!\(\(Clear[perpframe, r];\)\n > \(r = Random[Real, {\[Pi]\/4, \[Pi]\/2}];\)\n > \(s = Random[Real, {0, \[Pi]\/2}];\)\n > \(\(t = Random[Real, {\(-\(\[Pi]\/2\)\), \[Pi]\/2}];\)\(\n\) > \)\n > \(Clear[perpframe];\)\n > \(\({perpframe[1], perpframe[2], > perpframe[3]} = {{Cos[r]\ Cos[t] - Cos[s]\ Sin[r]\ Sin[t], > Cos[t]\ Sin[r] + Cos[r]\ Cos[s]\ Sin[t], > Sin[s]\ Sin[t]}, {\(-Cos[s]\)\ Cos[t]\ Sin[r] - Cos[r]\ Sin[t], > Cos[r]\ Cos[s]\ Cos[t] - Sin[r]\ Sin[t], > Cos[t]\ Sin[s]}, {Sin[r]\ Sin[s], \(-Cos[r]\)\ Sin[s], > Cos[s]}};\)\(\n\) > \)\n > \(ranger = 1.0;\)\n > \(frameplot = > Show[Table[ > Arrow[perpframe[k], Tail \[Rule] {0, 0, 0}, > VectorColor \[Rule] Indigo], {k, 1, 3}], > Graphics3D[Text["\<perpframe[1]\>", 0.5\ perpframe[1]]], > Graphics3D[Text["\<perpframe[2]\>", 0.5\ perpframe[2]]], > Graphics3D[ > Text["\<perpframe[3]\>", > 0.5\ perpframe[3]]], \[IndentingNewLine]Axes3D[ranger], > PlotRange \[Rule] {{\(-ranger\), ranger}, {\(-ranger\), > ranger}, {\(-ranger\), ranger}}, Boxed \[Rule] False, > Axes \[Rule] True, ViewPoint \[Rule] CMView, > AxesLabel \[Rule] {"\<x\>", "\<y\>", "\<z\>"}];\)\) --> Show::"gcomb" : "An error was encountered in combining the graphics objects in << 1 >>. More... > In[1230]:= > > > {xstretch,ystretch,zstretch} = {1.5,0.8,0}; > ranger=1.5; > {slow,shigh} = {0,Pi}; > {tlow,thigh} = {0,2 Pi}; > > hungellipsoidplot = > ParametricPlot3D[ xstretch Sin[s] Cos[t] perpframe[1] + > ystretch Sin[s] Sin[t] perpframe[2] + > zstretch Cos[s]perpframe[3], {s,slow,shigh},{t,tlow, > thigh}, > > PlotRange->{{-ranger,ranger},{-ranger,ranger},{-ranger,ranger}}, > Axes->True,AxesLabel->{"x","y","z"}, > Boxed->False,ViewPoint->CMView,DisplayFunction->Identity]; > > setup =Show[hungellipsoidplot,frameplot,Axes3D[ranger], > DisplayFunction->$DisplayFunction]; --> ViewPoint::"nlist3" : "(CMView) is not a list of three numbers. More... --> Show::"gcomb" : "An error was encountered in combining the graphics objects in [...snipped...] > In[1236]:= > > > Show[setup,ViewPoint->12 perpframe[1]]; --> Show::"gtype" : "(Show) is not a type of graphics. More... Jean-Marc