Groups | Search | Server Info | Keyboard shortcuts | Login | Register


Groups > comp.soft-sys.math.mathematica > #16818

Re: Three masses and four springs

From Bob Hanlon <hanlonr357@gmail.com>
Newsgroups comp.soft-sys.math.mathematica
Subject Re: Three masses and four springs
Date 2014-04-17 09:10 +0000
Message-ID <lio5p8$5jl$1@smc.vnet.net> (permalink)
References <20140416074027.11A846A2B@smc.vnet.net>
Organization Time-Warner Telecom

Show all headers | View raw


Use FullSimplify


DSolve[{
   -2*x1[t] + x2[t] == x1''[t],
   -2*x2[t] + x3[t] + x1[t] == x2''[t],
   -2*x3[t] + x2[t] == x3''[t],
   x1[0] == -1, x2[0] == 2, x3[0] == -1,
   x1'[0] == 0, x2'[0] == 0, x3'[0] == 0},
  {x1[t], x2[t], x3[t]}, t] //
 FullSimplify


{{x1[t] -> (1/2)*
         ((-1 + Sqrt[2])*
              Cos[Sqrt[2 - Sqrt[2]]*
                  t] - (1 + Sqrt[2])*
              Cos[Sqrt[2 + Sqrt[2]]*
                  t]),
  x2[t] -> (1/2)*
    ((-(-2 + Sqrt[2]))*
              Cos[Sqrt[2 - Sqrt[2]]*
                  t] + (2 + Sqrt[2])*
              Cos[Sqrt[2 + Sqrt[2]]*
                  t]),
  x3[t] -> (1/2)*
    ((-1 + Sqrt[2])*
              Cos[Sqrt[2 - Sqrt[2]]*
                  t] - (1 + Sqrt[2])*
              Cos[Sqrt[2 + Sqrt[2]]*
                  t])}}



Bob Hanlon




On Wed, Apr 16, 2014 at 3:40 AM, Robert Jenkins <dale.jenkins8@gmail.com>wrote:

> The instruction
> DSolve[{-2*x1[t] + x2[t] == x1''[t], -2*x2[t] + x1[t] == x2''[t],
>   x1[0] == -1, x2[0] == 2, x1'[0] == 0, x2'[0] == 0}, {x1, x2}, t]
> produces a simple solution. But I am surprised to find the three-mass
> version produces a mass of complication. Have I made a mistake?
> DSolve[{-2*x1[t] + x2[t] == x1''[t], -2*x2[t] + x3[t] + x1[t] ==
>    x2''[t], -2*x3[t] + x2[t] == x3''[t], x1[0] == -1, x2[0] == 2,
>   x3[0] == -1, x1'[0] == 0, x2'[0] == 0, x3'[0] == 0}, {x1, x2, x3},
>   t]
>
>

Back to comp.soft-sys.math.mathematica | Previous | Next | Find similar


Thread

Re: Three masses and four springs Bob Hanlon <hanlonr357@gmail.com> - 2014-04-17 09:10 +0000

csiph-web