from sympy import symbols, Function, Eq, dsolve

​

# Define the symbolic variables

x = symbols('x')

y = symbols('y', cls=Function)(x)

​

# Define the differential equation

diff_eq = Eq(y.diff(x, x) + 2*y.diff(x) + y, 0)

​

# Solve the differential equation and obtain the solution with steps

solution = dsolve(diff_eq, y, hint='all', simplify=False)

​

# Print the solution and the steps

for sol in solution:

    print("Solution:")

    print(sol)

    print("Steps:")

    print(sol.simplify())

    print()

import numpy as np

from scipy.integrate import odeint

​

# Define the differential equation

def differential_equation(y, t):

    return 2 * y + 1

​

# Set initial condition

initial_condition = 0

​

# Set time points to solve the equation

time_points = np.linspace(0, 10, 100)

​

# Solve the differential equation

solution = odeint(differential_equation, initial_condition, time_points)

​

# Print the solution

print(solution)