get solve() working again, add VoltageDerivativeConstraints