MathGroup Archive 1997

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

Search the Archive

Re: Plot3D precision limits?

  • To: mathgroup at
  • Subject: [mg6588] Re: [mg6504] Plot3D precision limits?
  • From: "Preferred Customer" <sherman.reed at>
  • Date: Fri, 4 Apr 1997 02:11:37 -0500 (EST)
  • Sender: owner-wri-mathgroup at

Although not specific answers to your question, I found that using
ContourPlot was useful in looking at the shape of the solution.

DensityPlot was not as useful.

I did not find a way to have the entire x or y scale printed out, but did
verify that the input xmax, xmin, etc. were observed.

sherman reed

> From: Jim Hicks <jim at>
> To: mathgroup at
> Subject: [mg6588] [mg6504] Plot3D precision limits?
> Date: Thursday, March 27, 1997 1:42 AM
> I have two questions regarding the statements listed below that find the
> numerical solution to a system of two equations and two unknowns (x and
> y).  This solution happens to be a maximum likelihood solution to the
> problem of estimating two unknown parameters, x and y.  The solution
> provided by FindRoot using Mathematica 3.0 under Solaris 2.5 for SPARC
> is the following:
> 	{x-> -0.237575444848816, y-> -0.0531098274657849}
> I am interested in visualizing the surface of the log-likelihood
> function in the immediate vicinity of the optimal solution.  As you can
> see, I have used Plot3D to plot this function near the solution.
> Question 1:
> As I increase the number of digits specified for the ranges of x and y
> to zoom in closer to the actual solution, Plot3D produces a non-smooth
> surface.  Plot3D produces a smooth surface if fewer digits are specified
> following the decimal for the plot ranges.  I assume that I am crossing
> some machine precision threshold.  Is this true and is there anyway to
> overcome it so that I can see an accurate representation of the surface
> near the solution given by x and y above?
> Question 2:
> If you reproduce the plot that I am considering, you will notice that
> all the axis labels are printed with up to 6 digits after the decimal. 
> Is it possible to change this by some option such that for example, as
> many as 15 digits would be displayed?
> Thank you very much for your consideration of my questions.
> Jim
> je-hicks at
> c={0,0,1,0,0,1,1,0,0,0,0,1,1,0,1,1,0,1,1,0,1}
> p=1/(1+Exp[-x+y(b-a)])
> FindRoot[{Sum[c[[i]]-p[[i]],{i,1,21}],
>           Sum[(c[[i]]-p[[i]])(a[[i]]-b[[i]]),{i,1,21}]},{x,0},{y,0}]
> Plot3D[Sum[c[[i]]Log[p[[i]]]+(1-c[[i]])Log[1-p[[i]]],{i,1,21}],
>           {x,-.23757544485,-.23757544484},{y,-.05310983,-.05310982},
>           PlotPoints->15, BoxRatios->{1.0,1.0,0.4}]

  • Prev by Date: Mathlink and VisualBASIC
  • Next by Date: Re: Default notebook window size
  • Previous by thread: Mathlink and VisualBASIC
  • Next by thread: Re: What is the funny-looking "i" in the palette?