[x,y] = ellipse(N10,delN10,N26,delN26,0,plotString) [h1] = ellipse(N10,delN10,N26,delN26,1,plotString) [h1,h2] = ellipse(N10,delN10,N26,delN26,2,plotString)
Plots an uncertainty ellipse on the [Be] - [Al] / [Be] diagram (the `Lal-Klein-Nishiizumi Al-Be diagram').
delN26 are the Be and Al concentrations and their uncertainties (atoms g).
The third argument can be 0, in which case the function returns the x and y coordinates of the uncertainty ellipse; 1, in which case the function plots the uncertainty ellipse and returns its handle; or 2, in which case the function plots and uncertainty ellipses and returns their handles. The default is 1.
plotString is an optional string telling the function what line type and color to use to draw the ellipse. See the documentation for the MATLAB function
plot for more information.
What this function actually does is calculate the joint probability on a mesh in :
It then creates a normalized cumulative distribution for to determine the values of that account for 68% and 95% of the total probability, and draws contours of the 2-d PDF at those values.