MathGroup Archive 2013

[Date Index] [Thread Index] [Author Index]

Search the Archive

Help needed on how plot a stereographic projection

  • To: mathgroup at smc.vnet.net
  • Subject: [mg130966] Help needed on how plot a stereographic projection
  • From: "Eduardo M. A. M. Mendes" <emammendes at gmail.com>
  • Date: Thu, 30 May 2013 06:14:24 -0400 (EDT)
  • Delivered-to: l-mathgroup@mail-archive0.wolfram.com
  • Delivered-to: l-mathgroup@wolfram.com
  • Delivered-to: mathgroup-outx@smc.vnet.net
  • Delivered-to: mathgroup-newsendx@smc.vnet.net

Hello

Although I have been using Mathematica for more than year, I feel that I haven't barely scratched the surface of what Mathematica can do.  

The following example gives the result that I need but the outcome is ugly and slow.

ClearAll[stereographicProjection];

stereographicProjection::usage="stereographicProjection[complexnumber] will return the stereoprojection of a complex point considering the Riemann sphere";

SyntaxInformation[stereographicProjection]={"ArgumentsPattern"->{_}};

stereographicProjection[complexnumber_]:=
Module[{a1,a2,a3},
If[ComplexExpand[Abs[complexnumber]]==Infinity,
a1=0;a2=0;a3=1,
=
a1=ComplexExpand[Re[complexnumber]]/(1+ComplexExpand[Abs[complexnumber]]^2);

a2=ComplexExpand[Im[complexnumber]]/(1+ComplexExpand[Abs[complexnumber]]^2);
=
a3=ComplexExpand[Abs[complexnumber]]^2/(1+ComplexExpand[Abs[complexnumber]]^2)];
{a1,a2,a3}
]

tab3=Table[stereographicProjection[(s+1)/(s^2 (s-1))/.{s-> I \[Omega]}],{\[Omega],-1000,1000,0.1}];

=
Show[ContourPlot3D[x^2+y^2+(z-1/2)^2==(1/2)^2,{x,-1,1},{y,-1,1},{z,0,1},Mesh->Automatic,AxesLabel-> =
{"x","y","z"},BoxRatios->{1,1,1/2},ImageSize-> =
Large],ListPointPlot3D[tab3,PlotStyle->Directive[PointSize[Large],Magenta],ImageSize-> =
Large],ListPointPlot3D[{stereographicProjection[-1]},PlotStyle->Directive[PointSize[0.02],Red]],ImageSize-> Large]

a) Is there another way of getting the same plot?
b) How to get the points of tab3 connected? 
c) How to change the opacity of the sphere?


Improvements, suggestion and critiscims are welcome.  

Many thanks

Ed




  • Prev by Date: Re: Warsaw Univ. course, was Re: Work on Basic Mathematica Stephen!
  • Next by Date: Re: Warsaw Univ. course, was Re: Work on Basic Mathematica Stephen!
  • Previous by thread: Re: Optimize evaluation of symbolic expressions
  • Next by thread: Re: Help needed on how plot a stereographic projection