Hello everyone,

I am having trouble with this a system of 12 differential equations with the respective 12 initial conditions. I try to usen dsolve with the "numeric" parameter but it stays "evaluating" forever. I tried a system of 10 equation and I got the correct answer, although the program stayed calculating for several minutes (more or less 5 minutes). Is the size of this system too much to handle for the program? 

This is the system of differential equations (sys) and the initial conditions (ics) that I am trying to solve:

maple_primes_question.mw

Download maple_primes_question.mw

What am I doing wrong? or is it that the program cannot handle such a large system? In that case, that would be dissappointing since I have to solve a similar system but with up to 22 differential equations and initial conditions so I am stuck...

Anyway, thank you very much for any guidance, and happy new year.

Hi everyone,

I have been trying to solve a system of 8 differential equations (1st order) with 8 different initial conditions. The system in question is the following.

sys:={

diff(a[0](t), t) - diff(a[1](t), t) + diff(a[2](t), t) - diff(a[3](t), t) = 0,

diff(b[0](t), t) + diff(b[1](t), t) + diff(b[2](t), t) + diff(b[3](t), t) = 0,

diff(a[0](t), t) + diff(a[1](t), t) + diff(a[2](t), t) + diff(a[3](t), t) - diff(b[0](t), t) + diff(b[1](t), t) - diff(b[2](t), t) + diff(b[3](t), t) = 0,

diff(a[1](t), t) + 4*diff(a[2](t), t) + 9*diff(a[3](t), t) - 3*diff(b[1](t), t) + 12*diff(b[2](t), t) - 27*diff(b[3](t), t) = 0,

diff(a[0](t), t) - diff(a[1](t), t)/2 - diff(a[2](t), t)/2 - 16*a[2](t) + diff(a[3](t), t) + 48*a[3](t) = 0,

diff(a[0](t), t) + diff(a[1](t), t)/2 - diff(a[2](t), t)/2 - 16*a[2](t) - diff(a[3](t), t) - 48*a[3](t) = 0,

diff(b[0](t), t) - diff(b[1](t), t)/2 - diff(b[2](t), t)/2 - 48*b[2](t) + diff(b[3](t), t) + 144*b[3](t) = 0,

diff(b[0](t), t) - diff(b[1](t), t)/2 - diff(b[2](t), t)/2 - 48*b[2](t) + diff(b[3](t), t) + 144*b[3](t) = 0}

and the initial conditions are

ics:={a[0](0) = 1/8, a[1](0) = 1/4, a[2](0) = 1/4, a[3](0) = 1/8, b[0](0) = 23/24, b[1](0) = 1/12, b[2](0) = -1/12, b[3](0) = 1/24}

When I run 

sols := dsolve(sys union ics, numeric)

I get the error message

Error, (in DEtools/convertsys) ODE system is insufficient to determine values for all dependent variables in the system

What is happening?

Thanks for the help.

Hello everyone,

I am trying to solve a system of six equations through a matrix and a vector. The matrix is 6x6, so the function "LinearSolve" should find a solution 

The matrix is

Matrix(6, 6, [[1., -1., 1., -1., 1., -1.], [1., 1., 1., 1., 1., 1.], [-2., 1.618033989, 30.94427190, -153.8246851, 371.1559479, -572.9674774], [-2., 0.6180339876, 22.94427191, -41.15905356, -37.04759741, 149.3606798], [-2., -0.6180339876, 13.05572810, 22.82468509, -20.15594802, -81.03252254], [-2., -1.618033989, 5.055728096, 28.15905368, 68.04759752, 104.6393203]])

and the vector associated with the system is

Vector[column](6, [-1/3, -1, 0, 0, 0, 0])

I am trying to solve this system with

a := LinearSolve(M3, v);

but it stays evaluating.

I have solved a similar system (5x5) with this function (LinearSolve) and it took less than a second, so I dont understand why it takes so long in this case.

Thank you in advance.

Hello

I am trying to calculate a definite integral of an absolute value function. I should get a positive result, but I end up getting a negative result. Why is this?

The line I am trying to run is this

And the result I get is:
 ET2 := -0.056854377998556975271421429744140962019176108843917

What am I missing?