# How can I obtain a unitary matrix of quantum circuit with angles as parametric?

I need to draw a quantum circuit in Clifford+T library and obtain automatically its transformation matrix. (I received a response to this part of my question before in the link: How can I obtain transformation matrix of a quantum circuit in the Qiskit?) and now I want to do this work as a parametric. For example, I want to define gate T as $$\begin{pmatrix} 1 & 0 \\ 0 & w \end{pmatrix}$$ where $$w=\frac{(1+i)}{\sqrt{2}}$$, and I want to obtain the transformation matrix of the circuit with parameter $$w$$ and not as numeric. How can I do this?

That being said, you can of course write your own script for a symbolic simulator: all you have to do is replace matrix multiplications with SymPy multiplications. A good starting point could be the quantum_info.Operator that uses straightforward NumPy multiplications and could probably be generalized easily.