Wrote flux_in_line, which seems to work (to first approximation).
Refactored some code relating to constants for 1a.
I wrote problem_1c()! It's got a break statement in it right now

but mostly works otherwise.

DAMN slow, though. Maybe the vectorize is slow?