Re: Re: problem getting the area under a parmetric curve
- To: mathgroup at smc.vnet.net
- Subject: [mg52898] Re: [mg52860] Re: problem getting the area under a parmetric curve
- From: Daniel Lichtblau <danl at wolfram.com>
- Date: Thu, 16 Dec 2004 03:40:34 -0500 (EST)
- References: <cpauof$ipf$1@smc.vnet.net> <200412150926.EAA10615@smc.vnet.net>
- Sender: owner-wri-mathgroup at wolfram.com
Roger Bagula wrote:
> Another of these distribution like pairs is:
> {(1-t^7)/(1+t^7),2^1/6*t*(7+35*t^14+21*t^28+t^42)^1/7/(1+t^7)}
> When fitted in the previous area programs,
> my computer doesn't compute the area of the normalization factor.
> n=6, 8 are again a squared circles.
> It appears that these can be done for general n case
> and the odd ones give the distributioin like parametric.
> [...]
I find it quite difficult to figure out what you want to do with this
pair of functions in 't', or what 'n' might refer to. I'm guessing you'd
like an implicit form? That's easy. As per another recent thread, just do
exprs = {(1-t^7)/(1+t^7),2^(1/6)*t*(7+35*t^14+21*t^28+t^42)^(1/7)/(1+t^7)};
InputForm[imp = First[GroebnerBasis[{x,y} - exprs, {x,y}, t]]]
Out[33]//InputForm= 2^(1/6) - 2^(1/6)*x^7 - y^7
Note that I altered your second expression so that the fractional
exponents would be parsed in the way I believe you desired.
The following plots should convince you that the implicit form is correct.
ParametricPlot[exprs, {t,0,4}]
ContourPlot[imp, {x,-1,1}, {y,0,1.2}, Contours->{0},
ContourShading->False]
It might make it easier for other readers to offer ideas (if that's what
you hope to obtain) if you clarify some remarks. For example:
(i) "When fitted in the previous area programs..."
When what is fitted? Fitted how? What previous programs?
(ii) "[M]y computer doesn't compute the area of the normalization factor."
What exactly are you asking of your computer? (Show a small code snippet
and explain clearly what you expect vs. what you obtain). Generally
speaking your computer will do as it is told, and without knowing that
detail it's difficult to comment or offer suggestions. I should also
point out that if you are using Mathematica and your expressions are
constructed in such a way that they will parse with operator precedence
other than what you intend, that falls under the heading of user error.
(iii) "n=6, 8 are again a squared circles"
Again, what is 'n. Do you mean
Clear[exprs]
exprs[n_] := {(1-t^n),
2^(1/(n-1))*t*Sum[Binomial[n,j]*t^(n*(j-1)),{j,1,n,2}]^(1/n)} /
(1+t^n)
If you do not clearly state details such as this, most readers will be
unable and/or unwilling to guess.
By the way:
ParametricPlot[exprs[6], {t,-10,10}]
InputForm[imp6 = First[GroebnerBasis[{x,y} - exprs[6], {x,y}, t]]]
Out[68]//InputForm= 2^(1/5) - 2^(1/5)*x^6 - y^6
ContourPlot[imp6, {x,-1,1}, {y,-2,2}, Contours->{0},
ContourShading->False]
The plots indicate the closed boxy curve you seem to expect. This closed
curve arises for even values of the exponent parameter 'n' because we
lose the pole at t=-1 (so the curve is continuous) and the limits of the
first coordinate are both -1 as t->+-Infinity, while the limits of the
second are both 0 (so the curve is closed).
My overall comment is this: the key to good posts is clarity (I don't
claim to be a paragon in this virtue, by the way). If you do not provide
sufficient code, or show code that is not going to parse the way you
like, or use terminology that is not standard and do not explain it,
then it is quite difficult to see whatever it is you wish to make known.
If in addition you have questions, then any lack of clarity in the post
decreases the likelihood that people might offer useful assistance.
Daniel Lichtblau
Wolfram Research
- References:
- Re: problem getting the area under a parmetric curve
- From: Roger Bagula <tftn@earthlink.net>
- Re: problem getting the area under a parmetric curve