substitute lambda calls for arrays of values. Recursion broken, but values diverge rapidly
avoid singularity
various typo corrections; add final computation step for plotting purposes
add cauchy value for F
fix sqrt(n) mistake, adds high frequency divergences in qt
add power conservation check; replace sum with np.sum
begin redraft
pythonic_main now runs, need to do sanity checks on all functions
add pudb debugger
tested integration
replace sum integration with trapz
forgot to specify polarization in pythonic_main mode solver
flip Bt in definition of F - should refer to second p argument
with sum method for integrals, both pythonic and normal method produce identical results
replaced F with cruedest solution to integral - system converges kind of but to wrong value
thorough testing of matching between fast and slow method
method working, adapting for multiple thicknesses
add dimension to support multiple thicknesses
remove test_beta
remove old methods