# Basic usage:
#
import numpy as np
from matplotlib import pyplot as plt
from nustattools import plotting as nuplt
rng = np.random.default_rng()
x = np.linspace(0, 10, 5)
u = x[:,np.newaxis] / 4
u[-2] *= -1
cov = np.eye(5) + u@u.T
err = np.sqrt(np.diag(cov))
y = rng.multivariate_normal(np.zeros(5), cov)
up = nuplt.wedgeplot(x, y, err, color="C2")
down = nuplt.wedgeplot(x, y, -err, color="C3")
down.set_facecolor("C1")