MathGroup Archive 2007

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

Search the Archive

Re: Foucault pendulum

  • To: mathgroup at smc.vnet.net
  • Subject: [mg80147] Re: Foucault pendulum
  • From: Jens-Peer Kuska <kuska at informatik.uni-leipzig.de>
  • Date: Tue, 14 Aug 2007 06:49:19 -0400 (EDT)
  • Organization: Uni Leipzig
  • References: <f9mq8f$r63$1@smc.vnet.net> <f9p4db$qdh$1@smc.vnet.net>
  • Reply-to: kuska at informatik.uni-leipzig.de

Hi,

sorry the last version has a bug
here is a nicer version

pendelPos[{x_, y_}] := {x, y, -Sqrt[10 - x^2 - y^2]}

pendel[tau_?NumericQ,
   sol_] := ({AbsoluteThickness[3],
      Line[{{0, 0, 0}, pendelPos[{x[t], y[t]}]}],
      Sphere[pendelPos[{x[t], y[t]}], 0.15]} /. sol[[1]]) /. t -> tau

Manipulate[
  DynamicModule[{traj, sol, fed, fulltraj},
   fde =
    {x''[t] == -\[Omega]^2*x[t] +
       2 \[CapitalOmega]*Sin[\[Phi]]*y'[t] ,
     y''[t] == -\[Omega]^2*y[t] -
       2 \[CapitalOmega]*Sin[\[Phi]]*x'[t]};
   sol = NDSolve[
     Join[fde, {x[0] == 2, y[0] == 2, x'[0] == 0, y'[0] == 0}],
     {x[t], y[t]}, {t, 0, 64 Pi}, MaxSteps -> Infinity];
   fulltraj = ParametricPlot3D[
     {x[t], y[t], -4} /. sol[[1]], {t, 0, 64 Pi},
     PlotPoints -> 1024
     ];
   traj = ParametricPlot3D[
     pendelPos[{x[t], y[t]} /. sol[[1]]], {t, t1, t1 + 4 Pi},
     PlotStyle -> RGBColor[1, 0, 0]
     ];
   Graphics3D[
    {traj[[1]], fulltraj[[1]], pendel[t1, sol]},
    PlotRange -> {{-2.5, 2.5}, {-2.5, 2.5}, {-4, 0}}
    ]], {{\[Omega], 1, "Pendel Frequence"}, 0.1,
   10}, {{\[CapitalOmega], 1/16, "Earth Circum Frequence"}, 0,
   1/2}, {{\[Phi], Pi/3, "Latitude"}, 0, Pi/2}, {{t1, 0, "time"}, 0,
   60 Pi}
  ]

Regards
   Jens

dimitris wrote:
> On 12    , 14:17, dimitris <dimmec... at yahoo.com> wrote:
>> Hello.
>> Does anyone have notebooks
>> demonstrating Foucault's pendulum?
>>
>> Thanks
>> Dimitris
> 
> I mean the motion of a Foucalt's pendulum of course.
> 
> Dimitris
> 
> 


  • Prev by Date: Re: Re: Reading Coordinates from Plot
  • Next by Date: Re: Indexed element treated as number?
  • Previous by thread: Re: Foucault pendulum
  • Next by thread: Re: Foucault pendulum