bec_t.py¶

from __future__ import print_function

import numpy as np

from bec import solve_bec


def main():
    tmin = 0
    tmax = 3
    n = 301
    for t in np.linspace(tmin, tmax, n):
        alpha, c = solve_bec(t)
        print(f"{t:.4f} {alpha:.5e} {c:.5e}")


if __name__ == '__main__':
    main()