Services & ResourcesWolfram ForumsMathGroup Archive

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

Re: How can I "perturbate" a NDSolve expresion?

• To: mathgroup at smc.vnet.net
• Subject: [mg102035] Re: [mg102000] How can I "perturbate" a NDSolve expresion?
• From: Iván Lazaro <gaminster at gmail.com>
• Date: Mon, 27 Jul 2009 05:55:20 -0400 (EDT)
• References: <200907260754.DAA18931@smc.vnet.net>

Thanks for the answers!

I need a more specific perturbation. I need something like this (with this I
expect to manipulate the frequency of the perturbation)


If[RandomReal[]<0.001,a=RandomReal[],a=0]


 NDSolve[{x1'[t] == (Sin[t] + a)*x2[t],
  x2'[t] == (Cos[t] - a*x1[t]^2)*x1[t], x1[0] == 1,
  x2[0] == 0}, {x1[t], x2[t]}, {t, 0, 10}]

But I need that the If condition evaluates in each time step. But I haven't
been able to introduce the conditional in the ndsolve statment.

• References:
  • How can I "perturbate" a NDSolve expresion?
    • From: Iván Lazaro <gaminster@gmail.com>