begin coupling coefficient for 1973 marcuse paper
add more parameters to kbar function
try writing the RHS function
expand complex equations to fully real expressions, i don't know how to cope with antik's foreign arrays
implement vandermonde problem, apply-evolution actually only makes a single step
code to check that the jacobian for the vandermonde equation is correct
simulation starts to resemble something that i would expect
something must be wrong, energy is not conserved
bug in e-mu-nu-m. there is no coupling between incident cos and spurious sin when fiber is bent (a cos distortion)
use macro to define an optimized function to calculate derivatives. it doesnt compute values with k=0. for 59 modes the calculation for 7mm fiber is speeded up from 153.7s to 19.9s
there are still issues, when the c have non-zero imaginary value
i didnt implement the complex multiplication between c_nu and exp(i (beta_mu-beta_nu) z) correctly
add complex multiplication to macro
save an animation of the field
make a second animation with bigger bend radius