Solve a differential equation with contraints

[Geoffrey Mégardon]

Hi!

I'm trying Sympy this morning and I tried to do something like this:

equations = [
Eq(a(t).diff(t), (I(t) - gamma*a(t))/theta ),
Eq(a(0), 0),
Eq(I(t), 30) ]

dsolve(equations, a(t))

I thought it would solve the first equation with the constraints defined in the two other equations.
But apparently, that doesn't work as I expected, it told me:

AttributeError: 'list' object has no attribute 'atoms'

I tried with the function solve(), but it told me it can't solve the third equation.

So, I'm wondering how that works when we want to solve an differential equation with constraints?
Can we then plot the solution?


--
MEGARDON Geoffrey