Skip to main content
3 events
when toggle format what by license comment
Jul 27, 2014 at 16:59 comment added unlikely @Szabolcs I think this approach is the most general we can use to allow a discontinuous function inside an NDSolve, if we don't know in advance what discontinuos functions exactly can be involved (UnitStep, HeavisideTheta, Boole, ...). Thanks. There are some issue? I'm curious too to understand why Mathematica behave this way...
Jul 27, 2014 at 15:50 comment added Szabolcs it's also possible to just make a numerical black box out of UnitBox, as in f[x_?NumericQ] := UnitBox[(x - 0.5) 0.5/0.2]. Then it works. But I'd like to understand why Mathematica thinks that it shouldn't attempt solving this equation when it knows that there are variables taking discrete values.
Jul 27, 2014 at 15:47 history answered Michael E2 CC BY-SA 3.0