Hello, for a PDE of this type du/dt + (A-B*x)*du/dx = 0 , where variable u(x) describes some probability distribution (i.e. must be > 0) what is the proper numerical scheme to integrate it? When I set u to normal distribution and let the system evolve, soon I get u(x) < 0 for some x. I need Lax-like method (http://www.cmth.ph.ic.ac.uk/people/a.mackinnon/Lectures/cp3/node29.html). Is this available in FiPy?
Thank you. _______________________________________________ fipy mailing list [email protected] http://www.ctcms.nist.gov/fipy [ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ]
