added large chunk of explanation of sympy functionality.  some trouble with solving for constants with ICs