Allow passing parameters to commensurate lyapunov computing

[ErikQQY]

Right now, the commensurate Lyapunov exponents computing routine doesn't support passing parameters, we need to set that up!

FractionalDiffEq.jl/src/lyapunov.jl

Lines 142 to 150 in f21d3e2

	function extend_fun(t, temp)
	temp=reshape(temp, ne, ne+1)
	result = zeros(ne)
	fun(result, temp[:, 1], t)
	for i=2:ne+1
	result = [result; Jfdefun(t, temp[:, 1])*temp[:, i]]
	end
	return result
	end