Services & Resources / Wolfram Forums
MathGroup Archive
*Archive Index
*Ask about this page
*Print this page
*Give us feedback
*Sign up for the Wolfram Insider

MathGroup Archive 1992

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

Search the Archive

solving large (53000 equations) systems of linear ODEs with Mma

  • To: mathgroup at
  • Subject: solving large (53000 equations) systems of linear ODEs with Mma
  • From: Roque Donizete de Oliveira <oliveria at>
  • Date: Wed, 25 Nov 92 14:29:50 -0500

I have a systems of 53000 linear ODEs:
                                       P'(t) = A P(t)
where P(t) is a vector of length 53000, and A is a very large matrix
(constant) with only 5-diagonal of non-zero elements.
I've solved this problem using standard Runge-Kutta (in C and Fortran)
and using the Krylov approximation method for the matrix exponential operator.
I'd like to do the same thing (i.e., solve the system above) using Mma, if possible.
So I ask you fellow Mma users: Has anyone used Mma for this purpose ? Can Mma
solve (numerically, of course) such a large number of linear ODEs ?

Thanks for any information.

   oliveria at
   Nuclear Eng Dept, The U. of Michigan

  • Prev by Date: time consumption FindRoot
  • Previous by thread: time consumption FindRoot