MathGroup Archive 2005

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

Search the Archive

Re: Problem with singularities in Plot3D

  • To: mathgroup at smc.vnet.net
  • Subject: [mg55354] Re: Problem with singularities in Plot3D
  • From: Peter Pein <petsie at arcor.de>
  • Date: Sun, 20 Mar 2005 04:11:47 -0500 (EST)
  • References: <d1gtdv$jh$1@smc.vnet.net>
  • Sender: owner-wri-mathgroup at wolfram.com

Steeve Brechmann (schumi) wrote:
> Hi,
>  
> I have a little problem with the Plot3D command of Mathematica.
> I want to plot the 3D graphics of :
>  
> eqt[q_,N_,M_]:= - (3*q/(8*Pi))*(Sum[(Cos[4*Pi*n*d]/(2*n*d))-(Sin[4*Pi*n*d]/(8*Pi*(n*d)^2))-(Cos[4*Pi*n*d]/(32*Pi^2*(n*d)^3)),{n,1,N}] - Sum[(Cos[4*Pi*Abs[z0-m*d]]/Abs[z0-m*d])-(Sin[4*Pi*Abs[z0-m*d]]/(4*Pi*Abs[z0-m*d]^2))-(Cos[4*Pi*Abs[z0-m*d]]/(16*Pi^2*Abs[z0-m*d]^3)),{m,-M,M}])
>  
> Plot3D[Evaluate[eqt[1,10,10]],{z0,0,3},{d,0.1,6},PlotPoints->50];
>  
> Like you see, my expression eqt[q,N,M] have some singularities at d = 0 and Abs[z0-m*d] = 0.
>  
> How can i say to Mathematica that it must ignore the singularities Abs[z0-m*d] = 0 ?
>  
> Thanks for your help.
>  
> Steeve Brechmann von Brechenmacher
> Physics Department, Laval University
> steevebrechmann at yahoo.ca
> 
> 

If it is sufficient, to plot for z0 from 10^-6 to 3, try
Plot3D[Evaluate[eqt[1,10,10]],{z0,10^-6,3},{d,1/10,6},PlotPoints->50];

-- 
Peter Pein
Berlin


  • Prev by Date: Re: filling in an array
  • Next by Date: Re: J/Link problem on Mac OS X
  • Previous by thread: Problem with singularities in Plot3D
  • Next by thread: String to numbers