Services & ResourcesWolfram ForumsMathGroup Archive

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

Re: why the inputs are not ODEs in NDSolve function

• To: mathgroup at smc.vnet.net
• Subject: [mg102679] Re: why the inputs are not ODEs in NDSolve function
• From: dr DanW <dmaxwarren at gmail.com>
• Date: Tue, 18 Aug 2009 06:10:05 -0400 (EDT)
• References: <h6b4i1$rmk$1@smc.vnet.net>

Typos, dude, typos.

First off, e is not the same as E, so:

Simgmoi[x_?(Length@# == 5 &)] := {1/(1 + e^-x[[1]]), 1/(
  1 + e^-x[[2]]), 1/(1 + e^-x[[3]]), 1/(1 + e^-x[[4]]), 1/(
  1 + e^-x[[5]])};

Should be:

Simgmoi[x_?(Length@# == 5 &)] := {1/(1 + E^-x[[1]]), 1/(
  1 + E^-x[[2]]), 1/(1 + E^-x[[3]]), 1/(1 + E^-x[[4]]), 1/(
  1 + E^-x[[5]])};

Also, in your initialization statement:

init = {W == ( {
      {0.2, 0.3, 0.1},
      {0.5, 0.7, 0.3},
      {0.5, 0.2, 0.1},
      {0.3, 0.2, 0.6},
      {0.7, 0.5, 0.1}
     } ), e[0] == {0.2, 0.3, 0.4}, r[0] == {0.1, 0.4, 0.5}};

Should be

init = {W[0] == ( { ...

This should get you going again.  I think AccuracyGoal->Infinity is
inadvisable, you might want to scale back on your expectations.

Daniel

• Follow-Ups:
  • Re: Re: why the inputs are not ODEs in NDSolve function
    • From: Haibo Min <yshdfeng@gmail.com>