MathGroup Archive 2007

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

Search the Archive

Re: pursuit curve (differential equations)

  • To: mathgroup at smc.vnet.net
  • Subject: [mg72858] Re: [mg72843] pursuit curve (differential equations)
  • From: "Josef Otta" <josef.otta at gmail.com>
  • Date: Tue, 23 Jan 2007 04:23:04 -0500 (EST)

Hi,
i think that you have some missing arguments in your equation. I tried to
repair it (q->q[t], p->p[t] etc.)and here is the output:

p[t_] := Sin[t];
q[t_] := Cos[t];
k = 10.;
t0 = 0;
t1 = 6Pi;
poc0 = {1, 2};

soln = {x[t],
    y[t]} /. NDSolve[{x'[t] == k*
      Sqrt[p'[t]^2 + q'[
          t]^2]*(p[t] - x[t])/
              Sqrt[(p[t] - x[t])^2 + (q[t] - y[t])^2], y'[t] == k*
                Sqrt[p'[t]^2 + q'[t]^2]*(y[t] - q[t])/
                    Sqrt[(p[t] - x[t])^2 + (q[t] - y[
                      t])^2], x[t0] == poc0[[1]], y[t0] ==
            poc0[[2]]}, {x[t], y[t]}, {t, t0, t1}][[1]]

Regards,
Josef Otta
http://home.zcu.cz/~jotta

2007/1/22, Trijezni Pijanac <trijezni.pijanac at gmail.com>:
>
> hi i am doing pursuit curve in mathematica.. for instance, fox is chasing
> a
> rabbit - rabbit has a certain defined path (for example - a circle
> {cos(t),sin(t)}. fox always heads directly toward the rabbit.
>
> k - relative speed fox/rabbit
> p,q - rabbit's path (for instance a circle {cos(t),sin(t)}
> x,y - fox's path
> t - time :)
>
> soln = NDSolve[
>       {
>         x'[t] == k.Sqrt[p'[t]^2 + q'[t]^2].(p - x[t])
>                 /Sqrt[(p - x[t])^2 + (q - y[t])^2],
>         y'[t] == k.Sqrt[p'[t]^2 + q'[t]^2].(y - q[t])
>                 /Sqrt[(p - x[t])^2 + (q - y[t])^2],
>         x[0] == poc0[[1]],
>         y[0] == poc0[[2]]
>         },
>       {x[t], y[t]}, {t, t0, t1}];
>
> but this wont work, any suggestions?
>
>


  • Prev by Date: Re: pursuit curve (differential equations)
  • Next by Date: Re: How install tools for mmade?
  • Previous by thread: Re: pursuit curve (differential equations)
  • Next by thread: Re: pursuit curve (differential equations)