import numpy as np
from functools import wraps, lru_cache
from abc import ABCMeta, abstractmethod
from scipy.integrate import solve_ivp
from scipy.stats import levy_stable
from dataclasses import dataclass
from typing import Protocol, runtime_checkable

class QuantizedHamiltonian(type):
    def __new__(cls, name, bases, dct):
        dct['_energy_spectrum'] = property(lambda self: np.diag(self.H))...