MathGroup Archive 2011

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

Search the Archive

Integrating a rotated function inbetween two circles

  • To: mathgroup at smc.vnet.net
  • Subject: [mg122989] Integrating a rotated function inbetween two circles
  • From: "Dr. Peter Klamser" <klamser at googlemail.com>
  • Date: Sat, 19 Nov 2011 06:47:59 -0500 (EST)
  • Delivered-to: l-mathgroup@mail-archive0.wolfram.com
  • References: <ja5f65$6a7$1@smc.vnet.net> <201111181249.HAA07910@smc.vnet.net>

Hi,

when you execute the following code you get a rotated parabolic function 
a x^2 rotated around the centre of coordinate system:

becher =.
becher[a_, x_] = a x^2;
rot = RotationTransform[(\[Theta]), {0, 0}] // FullSimplify
bechergedreht =.
bechergedreht[\[Theta]_, a_, r_, x_] = rot[{(x), becher[a, x - r]}]
x$fuer$r[a_, b_,
r_] = (x /. Solve[becher[a, x - r]^2 + x^2 == (r + b)^2, x])[[2]] //
Simplify;
ParametricPlot[
Table[bechergedreht[2 \[Pi] i/10, 1, 3, x], {i, 0, 9}], {x, 3,
x$fuer$r[1, 2, 3]}, PlotRange -> 1.01 {{-5, 5}, {-5, 5}},
Epilog -> {Circle[{0, 0}, 3], Circle[{0, 0}, 5]}]
{xstrich, ystrich} = bechergedreht[\[Theta], a, r, x];
loesungx = x /. Solve[xstrich == xs, x] // TrigFactor;
bechergedreht$xstrich[\[Theta]_, a_, r_, xs_] =
ystrich /. x -> loesungx // FullSimplify;
Print["Mit n wird für Table die Anzahl der Schritte festgelegt."]
Print["Mit n = 4 bekomme ich das erwartete Ergebnis."]
n = 4;
Plot[Table[
bechergedreht$xstrich[2 \[Pi] i/10, 1, 3, x], {i, 1, n}], {x, -4,
4}, AspectRatio -> 1, Epilog -> {Circle[{0, 0}, 3]},
PlotRange -> 1.01 {{-5, 5}, {-5, 5}}]
Print["Mit n = 5 bekomme ich keinen Plot mehr."]
n = 5;
Plot[Table[
bechergedreht$xstrich[2 \[Pi] i/10, 1, 3, x], {i, 1, n}], {x, -4,
4}, AspectRatio -> 1, Epilog -> {Circle[{0, 0}, 3]},
PlotRange -> 1.01 {{-5, 5}, {-5, 5}}]

If you fill water into the cups on the right side you get a water wheel.

How do I manage the rotated function most efficiently, that is the well 
known pair of functions regrading to the rotation vector?

As you see I developed the function bechergedreht$xstrich, solving the 
x´(x) part of the rotating function to x(x´) and inserting this into the 
y´(x) part of the rotation vector. But this does not handle to well.

How can I calculate the surface oft the cups on the right side?

Kind regards

Peter



  • Prev by Date: Re: How to force integers
  • Next by Date: What is the point of having Initializations in DynamicModule and Manipulate?
  • Previous by thread: Re: x and y plot finding maximum
  • Next by thread: Re: x and y plot finding maximum