differential equations - Analytical solution of system of ODE's python

Question

I am aware that is quite easy to find the solution of ODE's when you have some initial conditions in python. My question is, is it possible to find the general solution of a system of ODE's in python and plot it? Kind of what Mathematica would do? A minimum reproducible example, (but it requires intial conditions and we don't want that) :

import numpy as np
import matplotlib.pyplot as plt
from scipy.integrate import odeint

def tank(s,t):
    dhdt1 = s[0]
    dhdt2 = s[1]
    dhdt1 = y - x
    dhdt2 = 100*y
    dhdt = [dhdt1,dhdt2]
    return dhdt


t = np.linspace(0,1)
s0=[0,2]

s = odeint(tank,s0,t)


# plot results
plt.figure(1)
plt.plot(t,y[:,0],'b-')
plt.plot(t,y[:,1],'r--')
plt.xlabel('Time (hrs)')
plt.ylabel('Height (m)')
plt.legend(['h1','h2'])
plt.show()

question from:https://stackoverflow.com/questions/65870079/analytical-solution-of-system-of-odes-python

Answer 1

Waitting for answers