I have used comsol for solving PDEs, but I think you can use it for ODEs as well. The trick for COMSOL to adapt your equation to their models, i.e. construct the domain first and give you physical parameters then choose the numerical scheme for your model.

You also use Matlab (ODE45) sutable for solving linear/nonlinear ODEs.

Good luck

You can fit your equation into the predefined models in Comsol by changing their coefficients to be exactly the same as yours and then solve the equation in the domain. 

Actually i'm getting errors like:

Failed to find consistent initial values.

Last time step is not converged.

I don't know what to do here.

HI

I got a same error, when the initial value of dependent variables are so far from the real solution comsol attempt to give you this error

try to change your initial boundary conditions, and try

good luck

you can take a look at my model:

Deleted research item The research item mentioned here has been deleted

Article Impact of Size and Position of Lithium Metal Reference Elect...