MaxNeg8

import scipy as sc
import numpy as np
from scipy import integrate
from matplotlib import pyplot as plt
from matplotlib.ticker import (MultipleLocator, FormatStrFormatter, AutoMinorLocator, AutoLocator)


def osz(y, t, gamma1, gamma2, omega0):
    x, p = y
    dydt = [p, -(-gamma1 + gamma2 * x**2)*p - omega0**2 * x]