Sysode linear 1storder 3eq

This PR implements linear system of three equations of first order with coefficients being both constant and function of variable. This PR is also about implementing 'n' homogeneous equations of first order with constant coefficients.
Is there a reason this is `_linear_neq_order1_type1` rather than `_linear_3eq_order1_type6`?