[WIP] Add Common Interface Bindings

[ChrisRackauckas]

This adds LSODA to the common interface. For example, as seen in the tests, problems on things like 2-dimensional matrix are supported through solve:

using LSODA, DiffEqProblemLibrary
prob = prob_ode_2Dlinear
sol = solve(prob,LSODAAlg(),save_timeseries=false,saveat=[1/2])

Additionally, this allows LSODA to be the backend solver for sensitivity analysis, parameter estimation, etc.

Currently there are some details to fix up:

• If the package is named LSODA, then the algorithm type cannot be named LSODA. This is a bummer. Currently I named it LSODAAlg. Is there a better name?
• Setup with save_timeseries. Currently you can only use the saveat style of saving instead of letting it save at each timepoint it steps to. Can someone help me find out how to get the value of each timestep @sdwfrost ?

[coveralls]

Coverage decreased (-38.6%) to 52.817% when pulling bcfdcdf on ChrisRackauckas:master into 913574f on rveltz:master.

[coveralls]

Coverage decreased (-38.6%) to 52.817% when pulling bcfdcdf on ChrisRackauckas:master into 913574f on rveltz:master.

[ChrisRackauckas]

Switching the PR branch.