implement trapezoidal method
update test/transient
add sparse support for transient analysis

BUG: different results