New feature needed: get the unitary in the computational basis only

[Belaloui]

As the title says, we need a method that returns the unitary representation of the circuit in the computational basis.

This will only be needed in future updates.