03-093147.py

from astropy import coordinates as coords, time, units as u
from datetime import datetime, timezone
import numpy as np


now = time.Time.now()
wageningen = coords.EarthLocation(lat=51.9692*u.deg, lon=5.6654*u.deg, height=9*u.m)

for mins in 15 * np.arange(40):
    timept = now + mins * u.minute
    sun = coords.get_sun(timept).transform_to(coords.AltAz(obstime=timept, location=wageningen))
    print(timept, '{:.2f}, {:.2f}'.format(sun.az.degree, sun.alt.degree))