Dear guys,

I am resolving a system of ODE with Open Modelica through Python. All variables are zero at t=0, and V200=1 if t=365. Here one of the equations

der(V200) = k*(pi2*I200-(c2+epsilon0*G200)*V200+D2*(V100-V200));

simulate(model, startTime=0, stopTime=315360000, numberOfIntervals = 500, outputFormat="mat")

However, V200 = 1e-22 if t<365, when it should be zero.

I have changed the interpolation method and the number of intervals and the best approximation is dassl with 500.

Can you help me please?

Thank you!