Use tiny() instead of == 0
Add test_linregress.f90
Test the results
Add another test
Implement stderr_intercept and test it
Add another test
Tests
Refactor
Update analyze_md
Print diffusion