Quarto - A Hello World

print("Hello World!")
print("Here is Quarto Jupyter!")

import numpy as np
import matplotlib.pyplot as plt
from scipy.stats import norm

# Parameters for the Gaussian distribution
mu = 0       # mean
sigma = 1    # standard deviation

# Generate a range of x values
x = np.linspace(mu - 4*sigma, mu + 4*sigma, 1000)

# Compute the PDF using scipy.stats.norm
pdf = norm.pdf(x, mu, sigma)

# Plot
plt.figure(figsize=(8, 5), dpi=300)
plt.plot(x, pdf, label=f'N({mu}, {sigma}²)', color='blue')
plt.title('Probability Density Function of a Gaussian Distribution')
plt.xlabel('x')
plt.ylabel('Density')
plt.grid(True)
plt.legend()
plt.tight_layout()
plt.show()