FAULT TOLERANT PREPARATION OF QUANTUM POLAR CODES DESCRIPTION FIELD OF THE INVENTION The present invention concerns the field of quantum computation and more particularly of fault tolerant preparation of quantum error correcting code states. It relates more specifically to a system and a method of preparation of quantum polar code states for fault tolerant quantum computation. BACKGROUND OF THE INVENTION Quantum computers make use of quantum phenomena such as superposition and entanglement to perform computation. Through precise control of these phenomena, it is in principle possible for quantum computers to outperform their classical counterparts. Quantum computation is based on the manipulation of quantum bits or “qubits” which can be regarded as a superposition of the 1 and 0 states of a quantum physical variable. A qubit is the basic unit of quantum information, also called quantum state |^^ which corresponds to a superposition of basis states |0^ and |1^, as follows: | _^ ^{^} _{= ^} ^| ₀ ^{^} _{+ ^} ^| ₁ ^{^} ₍₁₎ where ^ and ^ are complex numbers satisfying the normalization constraint, ^| ^ ^{|^} + ^| ^ ^{|^} = 1. The quantum state |^^ is a vector in a complex linear vector space, known as the Hilbert space. The set { ^| 0 ^{^} , ^| 1 ^{^} } is an orthogonal basis of the Hilbert state, known as the computational basis. This basis is not unique. For example, another important basis is the phase basis, which corresponds to the set { ^| + ^{^} , ^| − ^{^} }, where ^| + ^{^} and ^| − ^{^} are orthogonal vectors defined as follows: A qubit ^| ^ ^{^} can also be written as a superposition of basis states ^| + ^{^} and ^| − ^{^} . The quantum state or qubit ^| ^ ^{^} in Eq. (1) can be expressed in the phase basis, as follows: The basis states can be extended to ^ qubits, where ^ > 1, by tensor product. For a binary vector ^ = ⁽ ^ _^ , … , ^ _^ ⁾ ∈ ^{ 0, 1 ^{}^} , let ^| ^ ^{^} = | ⁽ ^ _^ , … , ^ _^ ⁾ ^ ∶= ^| ^ _^ ^{^} . Then, the set { ^| ^ ^{^ |} ^ ∈ {0, 1} ^{^} } is the computational basis on ^ qubits. Any ^ qubit quantum state can be written as a superposition of the ^ qubit computational basis states: where ^ _^ are complex numbers satisfying the normalization condition ^∑ _^ ^| ^ _^ ^{|^} = 1. Similarly, the quantum state ^| ^ ^{^} can also be written as superposition of ^ qubit phase basis states, which are tensor products of ^| + ^{^} and ^| − ^{^} states. A quantum state ^| ^ ^{^} on ^ qubits is said to be entangled if it is not possible to write ^| ^ ^{^} as a tensor product of ^ single qubit quantum states, that is: |^^ ≠ |^ _^ ^ ⊗ |^ _^ ^ ⊗ ⋯ ⊗ |^ _^ ^ (5) where ^| ^ _^ ^{^} , ^| ^ _^ ^{^} , … , ^| ^ _^ ^{^} are single qubit states. Hence, entanglement refers to correlation between parts of a quantum system. It is a peculiar property of quantum systems, as it does not have a classical counterpart. The processing of quantum information is performed by applying quantum gates on qubits. Some examples of these quantum gates are Pauli, Hadamard, and CNOT gates. Pauli gates are a set of four quantum gates, denoted by ^, ^, ^, ^, which act on a single qubit. Their action in the computational basis is as follows: ^|0^ = |0^, ^|1^ = |1^ (6) ^ ^| 0 ^{^} = ^| 1 ^{^} , ^ ^| 1 ^{^} = ^| 0 ^{^} (7) ^ ^| 0 ^{^} = ^| 0 ^{^} , ^ ^| 1 ^{^} = − ^| 1 ^{^} (8) ^ ^| 0 ^{^} = ^ ^| 1 ^{^} , ^ ^| 1 ^{^} = −^ ^| 0 ^{^} (9) Note that ^ = ^ ^^, and ^ ^{^} = ^ ^{^} = ^ ^{^} = ^, where ^ is the identity operator in Eq. (6). From Eq. (7), Pauli ^ acts like a NOT gate in the computational basis. From Eq. (8), it also follows that ^ ^| + ^{^} = ^| − ^{^} and ^ ^| − ^{^} = ^| + ^{^} , hence Pauli ^ acts like a NOT gate in the phase basis. Pauli gates are extended to act on ^ qubits by tensor product. For example, ^ ⊗ ^ is a Pauli operator on two qubits. Let ^ _^ ^̅ be the set of all Pauli operators on ^ qubits, and ^ _^ ≔ {±1, ±^} × ^ _^ ̅ , the set of Pauli operators possibly with a ±1 or ±^ sign. Then ^ _^ is a mathematical group, which is referred to as the ^-qubit Pauli group. It is worth noting here that any two elements ^ _^ , ^ _^ ∈ ^ _^ either commute, that is, [^ _^ , ^ _^ ] ≔ ^ _^ ^ _^ − ^ _^ ^ _^ = 0, or anti-commute, that is, {^ _^ , ^ _^ } ≔ ^ _^ ^ _^ + ^ _^ ^ _^ = 0. The Hadamard gate, denoted by ^, acts on a single qubit. It maps a computational basis state to a phase basis state and vice-versa. Its action in the computational basis is as follow: ^ ^| ₀ ^{^} ₌ ^| ₊ ^{^} _{, ^} ^| ₁ ^{^} ₌ ^| ₋ ^{^} ₍₁₀₎ Fig. 1 represents a Controlled-NOT (CNOT) gate the CNOT _^→^ . The gate takes as input two qubits with control on the second qubit and target on the first qubit. Each horizontal wire carries a single qubit from left to right. In the computational basis, the action of the gate takes place at the target qubit while the control qubit is unaffected by this action. Precisely, the gate CNOT _^→^ acts similarly to the classical reversible XOR gate, denoted XOR _^→^ , in the computational basis, that is: CNOT _^→^ (|^^ ⊗ |^^) = |^ ⊕ ^^ ⊗ |^^ , ^, ^ ∈ {0, 1} (11) where ^ ⊕ ^ denotes the XOR (sum modulo 2) of binary values ^ and ^. The CNOT _^→^ gate and the classical XOR gate are represented by the same circuit except to the fact that the CNOT gate acts on two qubits, while the XOR gate acts on two bits. It is worth noting that CNOT _^→^ acts as the XOR _^→^ gate (with reversed control and target) in the phase basis. Precisely, let ^| 0 ^{^^} ≔ |+^ and ^| 1 ^{^^} ≔ ^| − ^{^} , then we have the following: Figs. 2A and 2B represent quantum measurement circuits on a qubit in computational and phases bases. The single wire on the input carries a single qubit while the double wire on the output carries a single classical bit. In general, a quantum measurement on a qubit is performed with respect to an orthogonal basis, and the measurement outcome gives classical information. After the measurement, the qubit collapses randomly into one of the basis states, depending on the measurement outcome. In particular, Fig.2A represents the measurement of the qubit ^| ^ ^{^} = ^ ^| 0 ^{^} + ^ ^| 1 ^{^} in the computational basis. The measurement output is 0 with probability ^| ^ ^{|^} and 1 with probability ^| ^ ^{|^} . If the measurement outcome is 0, then the output state is ^| 0 ^{^} and if the measurement outcome is 1, then the output state is |1^. The computational basis measurement is also known as the Pauli ^ measurement as |0^ and |1^ are eigenstates of the Pauli gate ^. Fig. 2B represents the measurement of the qubit ^| ^ ^{^} = ^ ^| 0 ^{^} + ^ ^| 1 ^{^} in the phase basis { ^| + ^{^} , ^| − ^{^} }. The measurement circuit is equivalent to first applying the Hadamard gate on ^| ^ ^{^} , and then measuing it in the computational basis. The phase basis measurement is also known as the Pauli ^ measurement as ^| + ^{^} and ^| − ^{^} are eigenstates of the Pauli gate ^. Although various technologies exist for implementing quantum computers, they all share the same shortcomings, namely that the qubits are affected by external noise and decoherence. Whereas bits in classical computers are materialized at the physical level by on/off states of transistor switches with high error margins, there is no such security for qubits. Indeed, the fragile superposition of states of a qubit may easily be disturbed by its environment and collapse, resulting in a loss of information. Quantum computers therefore fundamentally require error correction codes and fault tolerance at the physical level. Quantum error correcting codes entangle several physical qubits, which act as a logical qubit. Entanglement between physical qubits is used to protect the logical information from error. More precisely, entanglement defines a correlation between physical qubits in terms of their Pauli operators, and an error happening on physical qubits changes the correlation. It is possible to detect this change in correlation by doing joint quantum measurements, in a way that the logical information is not collapsed. The classical information learned by doing this measurement is called a ‘syndrome’. The extracted syndrome is given as an input to a classical decoder, which generates an estimate of the error that has happened. There are different types of quantum error correcting codes such as stabilizer codes and Calderbank-Steane-Shor (CSS) codes. A stabilizer code on ^ physical qubits is defined using a subgroup ^ of the ^-qubit Pauli group ^ _^ . The codespace corresponding to the stabilizer code is the subset of the ^-qubit Hilbert space, stabilized by the subgroup ^. A Pauli operator ^ ∈ ^ _^ stabilizes a quantum state ^| ^ ^{^} , if it is an eigenstate of ^ with the eigenvalue 1, that is: ^ ^| ^ ^{^} = ^| ^ ^{^} (13) A subset ^ of quantum states is said to be stabilized by a subgroup ^ ⊂ ^ _^ if every element ^ ∈ ^ stabilizes every quantum state ^| ^ ^{^} ∈ ^. Note that for ^ to be non-empty, it is sufficient to have −^ ∉ ^ and all the elements in ^ commute with each other, meaning that, for any two elements ^ _^ , ^ _^ ∈ ^, we have: [^ _^ , ^ _^ ^] = 0 (14) The subgroup ^ can be completely specified by a generating set ^ = {^ _^ , ^ _^ , … , ^ _^ }. A generating set is independent if any ^ _^ ∈ ^ cannot be written as a product of elements from The size of an independent generating set determines the number of logical qubits encoded by the stabillizer code. Precisely, if the number of elements in an independent generating set is equal to ^ − ^, then the stabilizer code encodes ^ qubits. When ^ = 0, the code does not encode any quantum information, as it has only one fixed quantum state in its codespace, called a stabilizer state. Stabilizer codes are suitable for detecting Pauli errors. Consider a code state |^^, on which a random ^-qubit Pauli error ^ ∈ ^ _^ happens. Since any two elements of ^ _^ either commute or anti-commute, then, for any ^ _^ ∈ ^, we have the following: ^ _^ ^|^^ = (−1) ^{^} ^|^^ (15) where ^ = 0 if ^ _^ commutes with ^, and ^ = 1 if ^ _^ anti-commutes with ^. Therefore, if ^ _^ anti- commutes with ^, Eq. (15) implies that the error corrupted state ^ ^| ^ ^{^} is an eigenstate of ^ _^ , with eigenvalue −1. This means that we can detect the error by doing the Pauli measurement corresponding to the generator ^ _^ . The syndrome measurement of stabilizer codes corresponds to measuring all the generators ^ _^ , ^ _^ , … , ^ _^ . Based on the extracted syndrome, an estimate ^ ^{^} of ^ is then generated using a classical decoder. The CSS codes are an important subclass of stabilizer codes. A stabilizer code is a CSS code if there exists a generating set ^ = ^ _^ ∪ ^ _^ of the stabilizer group, such that any ^ _^ ∈ ^ _^ can be written as a tensor product of ^ and ^, that is, ^ _^ = ^ ^{^} ∶= ^ ^{^^} ⊗ ^ ^{^^} ⊗ ⋯ ⊗ ^ ^{^^} , for some ^ = ⁽ ^ _^ , ^ _^ , ⋯ , ^ _^ ⁾ ∈ ^{ 0, 1 ^{}^} , and similarly any ^ _^ ∈ ^ _^ can be written as a tensor product of ^ and ^, that is, ^ _^ = ^ ^{^} ∶= ^ ^{^^} ⊗ ^ ^{^^} ⊗ ⋯ ⊗ ^ ^{^^} for some ^ = ⁽ ^ _^ , ^ _^ , ⋯ , ^ _^ ⁾ ∈ {0, 1} ^{^} . Since, ^ _^ and ^ _^ must commute with each other, this imposes the following constraint: The CSS code may be associated with two classical codes on ^ bits, with parity check matrices ^ _^ and ^ _^ , where ^ _^ is a binary matrix whose rows are vectors ^ ∈ {0, 1} ^{^} such that ^ ^{^} ∈ ^ _^ , and ^ _^ is a binary matrix whose rows are vectors ^ ∈ {0, 1} ^{^} , such that ^ ^{^} ∈ ^ _^ . Then, Eq. (13) is equivalent to: ^ _^ ^^ ^ = 0 (mod 2) ⁽¹⁷⁾ Similarly, CSS quantum polar codes are constructed on the basis of classical polar codes. Figs.3A and 3B represent the encoding of classical polar codes using reversible XOR gates. The encoding of classical polar codes is done by applying the reversible XOR gate recursively on an ^-bit input ^ = ⁽ ^ _^ , ^ _^ , ⋯ , ^ _^ ⁾ ∈ ^{ 0, 1 ^{}^} , where ^ = 2 ^{^} , ^ > 0. For a set of positions ℱ ⊂ {1, … , ^ − 1}, the corresponding component ^ _^ ∈ ^{ 0, 1 ^}|ℱ| of the input vector ^ is frozen (i.e. fixed). We may take ^ _^ to be any vector in ^{ 0, 1 ^}|ℱ| , but it should be known to both the encoder and decoder. The set ℱ is called the ‘frozen set’. The remaining positions ℐ ∶= {1, … , ^ − 1} ∖ ℱ are used to encode bits. The set ℐ is called the ‘information set’. In the following, we denote by ^ ⁽ ^, ^, ^ _^ ⁾ , the classical polar code of codelength ^, frozen positions ^, and frozen vector ^ _^ ∈ ^{ 0, 1 ^}|ℱ| . The action of the reversible XOR gate XOR _^→^ on ^ = (^ _^ , ^ _^ ) ∈ {0, 1} ^{^} gives ^ ^{^} = ⁽ ^ _^ ⊕ ^ _^ , ^ _^ ⁾ . The vector ^ ^{^} can be expressed as ^ ^{^} = ^ _^ ^, where ^ _^ is the following matrix: Classical polar transform, that is, the recursive application of XOR _^→^ on ^ = 2 ^{^} qubits, is given by the matrix ^ ⊗^ ^ = ^ _^ . We note that the action of the opposite XOR, i.e., XOR _^→^ is described by ^^ ^ , i.e. the transpose of ^ _^ . Hence, the recursive application of XOR _^→^ is described by ^^ ^ . For any vector of information bits ^ _ℐ ∈ ^{ 0, 1 ^{}|^^ |} , ^ = ^ _^ ⁽ ^ _^ , ^ _ℐ ⁾ ∈ ^{ 0, 1 ^{}^} is a codeword of the polar code ^ ⁽ ^, ^, ^ _^ ⁾ . Let ^(^) ^ be the ^th column of ^ _^ , for ^ ∈ {1, … , ^ − 1} and let ^ _^ be the all zero vector. Then, the classical polar code ^ ⁽ ^, ^, ^ _^ ⁾ is generated by the set ^ = {^(^) ^ ∣ ^ ∈ ℐ}. Classical polar codes have an efficient decoder known as ‘successive cancellation’ (SC) decoder. SC decoder takes as input the frozen vector ^ _^ , and a noisy version of the codeword ^ ⊕ ^, where ^ = ^ _^ ⁽ ^ _^ , ^ _ℐ ⁾ is a codeword, and ^ ∈ ^{ 0, 1 ^{}^} is a random error, and it outputs an estimate ^ _ℐ of ^ _ℐ . The polar transform ^ _^ , for any ^ = 2 ^{^} , ^ ≥ 0, relates to ^^, as depicted in Fig. 3A. ^ Similarly, ^^ relates to ^^, and ^^ relates to ^^ and so on. Hence, doing this recursively, we can ^ ^ ^ ^ express ^ _^ only in terms of the XOR gate. Fig. 3B represents the polar transform of a classical polar code of codelength ^ = 2 ^{^} encoding 5 bits and where 3 bits are frozen. In this example, the frozen set is ℱ = {1,2,3} and the frozen vector is ^ _^ = ⁽ 0,0,0 ⁾ . In what follows, we use the notation |0 ^{^} ^ ≔ |+^ and |1 ^{^} ^ ≔ |−^ to denote phase basis states. For ^ = (^ _^ , … , ^ _^ ) ∈ {0, 1} ^{^} , we define |^^ = |^^ _^ ^^ ⊗ |^^ _^ ^^ ⊗ ⋯ ⊗ |^^^ _^ ^^. Moreover, we define ^ ^{^} ≔ ^ ^{^^} ⊗ ^ ^{^^} ⊗ ⋯ ⊗ ^ ^{^^} , and similarly ^ ^{^} ≔ ^ ^{^^} ⊗ ^ ^{^^} ⊗ ⋯ ⊗ ^ ^{^^} . Fig. 4 represents the encoding of CSS quantum polar codes using the quantum CNOT gates. The encoding of CSS quantum polar codes is done by applying the quantum CNOT gate recursively on an ^ qubit quantum state ^| ^ ^{^} _^ , For a subset of positions ^ ⊂ ^, the input quantum state is frozen to a computational basis quantum state ^| ^ ^{^} _^ , where ^ ∈ ^{ 0, 1 ^{}|^|} , and for another subset ^ ⊂ ^, it is frozen to a phase basis state ^| ^ ^{^} _^ , where ^ ∈ ^{ 0, 1 ^{}|^|} . For ^ and ^, we may take any vectors in ^{ 0, 1 ^{}|^|} and ^{ 0, 1 ^{}|^|} , respectively, but they should be known to both the encoder and decoder. The remaining subset ^ ≔ ^ ∖ (^ ∪ ^) is used to encode a quantum state ^| ^ ^{^} _^ . Hence, the uncoded quantum state |^^ _^ can be written as: where ^| ^ ^{^} _^ is an arbitrary quantum state that we want to encode. Then, the encoded quantum state is given by ^ _^ ^| ^ ^{^} _^ , where ^ _^ denotes the quantum polar transform on ^ qubits, that is, the quantum operator on ^ qubits defined by the recursive application of the CNOT gate. A quantum polar code on ^ qubits, with frozen indices ^ and ^ corresponding to the computational and phase bases, respectively, and with frozen quantum states |^^ _^ and |^^ _^ is denoted by ^ ⁽ ^, ^, ^, ^| ^ ^{^} _^. , ^| ^ ^{^} _^ ⁾ . In particular, the circuit represented in Fig. 4 concerns a quantum polar code ^ ⁽ ^, ^, ^, ^| ^ ^{^} _^. , ^| ^ ^{^} _^ ⁾ , with ^ = 2 ^{^} , ^ = ^{ 1,2,3 ^} , ^ = ^{ 6,7,8 ^} , ^| ^ ^{^} _^ = ^| 0,0,0 ^{^} , and ^| ^ ^{^} _^ = ^{| ^} 0 ^{^} , ^{^} 0 ^{^^} ,0 ^{^} _^. It is known that CNOT _^→^ acts as the reversible XOR gate XOR _^→^ in the computational basis, while it acts as XOR _^→^ in the phase basis. Hence, the quantum polar transform acts as classical polar transform in the computational basis, while it acts as the opposite polar transform in the phase basis. The polar transform and the opposite polar transform are described by the matrices ^ _^ and ^^ ^ , respectively. Hence, for computational and phase basis states corresponding to ^ ∈ {0,1} ^{^} , the encoded quantum state ^ _^ ^| ^ ^{^} _^ can be expressed as: ^ _^ ^| ^ ^{^} _^ = ^ ^{^} ^ ^{^} _^ ^{^ ^^} _^ ^{^} _^ ^ (21) To define the ^ and ^ type stabilizer generators of the polar code ^ ⁽ ^, ^, ^, ^| ^ ^{^} _^. , ^| ^ ^{^} _^ ⁾ , let ^^ ^ ∈ ^{ 0, 1 ^{}^} be the binary vector, with 0 everywhere except 1 at the ^th position. Hence, the uncoded quantum state ^| ^ ^{^} _^ in Eq. (19) is stabilized by the following Pauli operators: Then, it follows that the encoded quantum state ^ _^ ^| ^ ^{^} _^ is stabilized by the following operators: (−1 ^{)^^} ^ _^ ^ ^{^} ^ ^ ^ _^ , ∀^ ∈ ^ (24) (−1 ^{)^^} ^ _^ ^ ^{^} ^ ^ ^ _^ , ∀^ ∈ ^ (25) The above operators are generators of the stabilizer group of the code ^(^, ^, ^, |^^ _^. , |^^ _^ ). The sandwiching actions of the CNOT gate CNOT _^→^ on two qubit Pauli ^ and ^ type operators are expressed as follows: CNOT _^→^ ^ ^{^} CNOT _^→^ = ^ ^{^^^} , ^ ∈ {0, 1} ^{^} (27) Since ^ _^ is simply the recursion of CNOT _^→^ , the sandwiching action of ^ _^ is described by classical polar transforms ^^ ^ and ^ _^ , respectively, on ^ and ^ type operators. This means that the ^ and ^ type generators of ^(^, ^, ^, |^^ _^. , |^^ _^ ) in Eqs. (24) and (25) can be written as follows: ^ ^{^} ^ ^ ( ₋₁ ^{) ^} _^ ^{^^} _{^ , ^ ∈ ^} ⁽²⁸⁾ (−1) ^{^} ^ ^{^} ^ ^ ^{^^^} _{, ^ ∈ ^} ⁽²⁹⁾ Therefore, the ^ type stabilizers of the quantum polar code ^ ⁽ ^, ^, ^, ^| ^ ^{^} _^ , ^| ^ ^{^} _^ ⁾ are given by the columns corresponding to the set ^, and the ^ type stabilizers are given by the columns of ^ _^ , corresponding to the set ^. Note from Eqs. (28) and (29) that the vectors ^, ^ cause only a sign factor in the stabilizer generators. Hence, two quantum polar codes that are defined by the same frozen set ^ (respectively, ^), but by different frozen states ^| ^ ^{^} _^ (respectively, ^| ^ ^{^} _^ ) have the same ^ (respectively, ^) type generating set, up to a sign factor. For this reason, two stabilizer codes having the same sets ^, ^, are said to be ‘equivalent stabilizer codes’. A quantum polar code ^(^, ^, ^, |^^ _^. , |^^ _^ ) as defined above can be used to encode information. However, the preparation of a polar code state according to the implementation of the quantum polar transform as described in relation to Fig. 4 is not fault tolerant, under the effect of noise. In fact, an error caused by a noise on a qubit may propagate to many qubits through the CNOT gates, applied during the implementation. Therefore, the implementation of the quantum polar transform according to the prior art cannot be used for fault tolerant quantum computations. An object of the present invention is to remedy the aforementioned drawbacks by proposing a circuit for preparing quantum polar code states in a simple manner that can be efficiently used for fault tolerant quantum computations. BRIEF DESCRIPTION OF THE INVENTION The present invention concerns a method of preparation of quantum polar code states in a quantum computing system materializing qubits, quantum gates and quantum circuits, for fault tolerant quantum computation, comprising the following steps: -prepare an initial quantum system ^ = {1, … , ^}, of ^ = 2 ^{^} single-qubit states associated to an initial quantum base, -prepare a desired quantum polar code state ^| ^ _^ ^{^^} _^ of codelength ^ = 2 ^{^} , referenced by predetermined first and second sets of frozen indices ^ = with respect to first and second quantum bases respectively, said quantum polar code state ^| ^ _^ ^{^^} _^ being prepared by recursively using two qubit Pauli measurement ^ ⊗ ^ circuits wherein, at each recursive level ^, where ^ = 1 ^^ ^, prepare a set of 2 ^{^⁄} 2 ^{^} quantum polar code states = 1 ^^ 2 ^{^⁄} 2 ^{^ ^} of codelengths 2 ^{^} , referenced by corresponding sets ^ _^(^) of indices comprising first and second sets of frozen indices = ^(^ + 1) _^(^) , … , 2 ^ ^ _(^) ^, each quantum polar code state ^^ _^ ^{^} ^ ^{^} ^(^)being prepared by applying two qubit Pauli measurement ^ ⊗ ^ circuits on two equivalent polar code states ^^ ^ ^ ^{^^^} belonging to the output of the antecedent recursive level ^ − 1 and referenced by two corresponding sets of indices ^ ^ ^ _(^^^) and ^ ^ ^ _(^^^) each of which comprising first and second sets of frozen indices ^ _^(^^^) and corresponding to first and second quantum bases respectively, the antecedent of the first recursive level ^ = 1 being said initial quantum system of single qubit states. The method according to the present invention enables to prepare in a simple and robust manner, polar code states that can be used in a noiseless scenario as well as in a noisy scenario. According to the present method of preparation of quantum polar code states, any eventual error on a quantum state at any recursive level do not propagate to an increased number of qubits. Various preferred embodiments of the present invention are defined in the appended claims. In particular, the fault tolerance under the effect of noise is achieved using two different estimation techniques. The first one is based on an error detection technique and the second is based on an error correction technique. The fault tolerant quantum polar code states are important to quantum computation. BRIEF DESCRIPTION OF THE DRAWINGS The present invention will be better understood from the description of the following embodiments, by way of illustration and in no way limitative thereto: Fig. 1, already described, represents a Controlled-NOT gate the CNOT _^→^ known from the prior art; Figs.2A and 2B, already described, represent quantum measurement circuits on a qubit in computational and phases bases, known from the prior art; Figs.3A and 3B, already described, represent the encoding of classical polar codes using reversible XOR gates, known from the prior art; Fig.4, already described, represents the encoding of CSS quantum polar codes using the quantum CNOT gates, known from the prior art; Fig. 5 schematically represents a quantum computing system for the preparation of quantum polar code states, according to a preferred embodiment of the present invention; Figs.6A and 6B schematically represent Pauli ^ ⊗ ^ and ^ ⊗ ^ measurement circuits on two qubits ^ _^ and ^ _^ , used in Fig.5; Fig. 7 schematically represents a method of preparation of quantum polar code states according to a preferred embodiment of the present invention; Figs. 8A and 8B, schematically represent the preparation of polar code states of length ^, given two equivalent polar code states of length ^/2, according to an embodiment of the present invention; Figs. 9A, 9B, 9C and 9D schematically represent different examples of fault tolerant preparation of quantum polar code states, according to an embodiment of the present invention; and Fig.10 schematically represents a computing system according to an embodiment of the present invention. DETAILED DISCLOSURE OF PARTICULAR EMBODIMENTS The present invention proposes a measurement-based preparation of fault tolerant quantum polar code states, using only two qubit quantum measurements. Fig. 5 schematically represents a quantum computing system for the preparation of quantum polar code states to be used for fault tolerant computation, according to a preferred embodiment of the present invention. The quantum computing system 1 comprises different types of quantum gates 3 that are interconnected by wires 5 to form quantum circuits 7 and in particular, Pauli measurement circuits. The wires 5 carry qubits around the circuits, while the quantum gates 3 execute some operations on the qubits to make quantum computations. This quantum computing system uses single Pauli measurement circuits for example of the ^ and ^ types (see Figs.2A and 2B) as well as two qubit Pauli ^ ⊗ ^ measurement circuits for example of the ^ ⊗ ^ type or ^ ⊗ ^ type (see Figs.6A and 6B). As explained below in relation to Fig. 7, single-qubit measurements are used at an initialization step ^ = 0, as a way to prepare an initial quantum system of ^ = 2 ^{^} single-qubit states associated to an initial quantum base. Two-qubit measurements are used at each subsequent step ^ = 1 ^^ ^, as a way to prepare a set of quantum polar code states, of codelength 2 ^{^} . Figs.6A and 6B schematically represent Pauli ^ ⊗ ^ and ^ ⊗ ^ measurement circuits on two qubits ^ _^ and ^ _^ . In Fig. 6A the right-hand side is a quantum circuit used to implement the Pauli ^ ⊗ ^ measurement on two qubits and ^ _^ and the left-hand side is simply a shorthand notation for the circuit on the right-hand side. In practice, Pauli ^ ⊗ ^ measurement may be implemented by the circuit on the right hand side, using an ancilla qubit ^. Here, qubits ^ _^ and ^ _^ are first entangled with A, by applying the CNOT gates CNOT _^^→^ and CNOT _^^→^ . Then, a Pauli-^ measurement is done on ^. The resulting procedure is a Pauli ^ ⊗ ^ measurement on qubits ^ _^ and ^ _^ . In Fig.6B the left-hand side is a shorthand notation for the circuit on the right-hand side which is used to implement the Pauli ^ ⊗ ^ measurement on two qubits ^ _^ and ^ _^ . In practice, Pauli ^ ⊗ ^ measurement may be implemented by the circuit on the right hand side, using an ancilla qubit ^. At first, a Hadamard gate is applied on ^. Then it is entangled with qubits ^ _^ and ^ _^ , by applying the CNOT gates CNOT _^→^^ and CNOT _^→^^ . Then another Hadamard gate is applied back on ^, which is then followed by a Pauli-^ measurement. The resulting procedure is a Pauli ^ ⊗ ^ measurement on qubits ^ _^ and ^ _^ . Various technologies exist for materializing or implementing qubits, quantum gates, quantum circuits and quantum computers. One technology is based on the energy levels of ions trapped in an electric or magnetic field at a temperature near absolute zero using also laser pulses, optical pumping, etc. Another technology may use nuclear magnetic resonance where transformations may be constructed from magnetic field pulses applied to spins in a strong magnetic field, etc. Other technologies use physical systems based on small semiconductors called quantum dots bounding the spin of electrons. Other systems may take advantage of electrons or ions trapped in synthetic diamonds. Different examples of physical systems materializing qubits, quantum gates and quantum circuits can be found in the reference book entitled “Quantum computation and quantum information” authored by M.A. Nielsen and I.L. Chuang, Cambridge University Press, 2016. The quantum circuits according to the present invention, comprise different sets of Pauli measurement circuits configured to recursively prepare a fault tolerant quantum polar code |^ _^ ^{^} ^ _^ of codelength ^ = 2 ^{^} , where ^ = {1, … , ^} denotes a system of ^ qubits. These sets are organized in a sequential order of different levels adapted to recursively prepare a quantum polar code ^| ^ _^ ^{^^} _^ of codelength ^ = 2 ^{^} . At the first level ^ = 0, a set of single qubit Pauli measurement circuits 3 are used to prepare an initial quantum system of ^ = 2 ^{^} single-qubit states, where each qubit is prepared in an initial quantum base. For example, each qubit is prepared as a ^| 0 ^{^} in a Z-quantum-basis. At each subsequent level ^, a set of two qubit Pauli measurement ^ ⊗ ^ circuits are used to prepare a set of quantum polar code states, of codelength 2 ^{^} , ^ = 1 ^^ ^. Lastly, at the final level, the quantum polar code ^| ^ _^ ^{^^} _^ of codelength ^ = 2 ^{^} is prepared. The details of this preparation is explained in relation to the method of Fig.7. Fig. 7 schematically represents a method of preparation of quantum polar code states according to a preferred embodiment of the present invention. Step E1 concerns the preparation of an initial quantum system ^ of ^ = 2 ^{^} single-qubit states ^{| ^ _^ ^{^} _^ , … , ^| ^ _^ ^{^} _^ ^} associated to an initial quantum base. The initial quantum base may for example, be a computational base ^| ^ ^{^} where ^ ∈ ^{ 0,1 ^{}^} or a phase base ^| ^ ^{^} where ^ ∈ ^{ 0,1 ^{}^} . The initial quantum system of ^ single-qubit states may be prepared by applying a single qubit Pauli ^ measurement circuit on each qubit of an ^ qubit quantum state. Steps E2-E5 describe a recursive preparation of a desired quantum polar code |^ _^ ^{^} ^ _^ by using two qubit Pauli measurement ^ ⊗ ^ circuits. The quantum polar code |^ _^ ^{^} ^ _^ is referenced by predetermined first and second sets of frozen indices ^ and ^with respect to computational and phase bases, with corresponding frozen values ^| ^ ^{^} _^ and ^| ^ ^{^} _^ . As an embodiment of the present invention, CSS polar codes ^ ⁽ ^, ^, ^, ^| ^ ^{^} _^. , ^| ^ ^{^} _^ ⁾ that encodes only one information qubit, that is, ^| ^ ^| + ^| ^ ^| = ^ − 1, are considered. Such a code is constructed as follows. An index ^ ∈ ^ ≔ {1, … , ^} is selected to place the information qubit. All the indices before it belong to the frozen set ^, that is, ^ = {1, … , ^ − 1}, and all the indices after it belong to the frozen set ^, that is, ^ ≔ {^ + 1, … , ^}. The present invention is mainly interested in preparing a logical code state, and therefore, the information qubit corresponding to the selected index ^ ∈ ^ may be fixed in either one of the bases. The present invention does not depend on which basis the information qubit is fixed. For example, it may be chosen to fix the information index ^ ∈ ^ in the computational basis, i.e., the information qubit may be either |0^ or |1^. In this case, it may be assumed that ^ ∈ ^, meaning that, ^| ^ ^| + ^| ^ ^| = ^. Thus, according to this example, the indices in the set ^ = {1, … , ^ − 1, ^} are frozen in some computational basis state ^| ^ ^{^} _^ , ^ ∈ {0, 1} ^{^} , and the indices ^ ≔ {^ + 1, … , ^} are frozen in some phase basis state ^| ^ ^{^} _^ , ^ ∈ ^{ 0, 1 ^{}^^^} . In other words, the desired quantum polar code to be prepared may be defined as follows: where ^ ∈ ^{ 0,1 ^{}^} and ^ ∈ ^{ 0,1 ^{}^^^} are known. It is to be noted that according to Eqs. (28) and (29) described in the background section, the quantum state given in Eq. (30), is a stabilizer state, whose ^ and ^ type stabilizer groups are given by: The above Eqs. (31) and (32) imply that any two polar code states, |^^ ^ ^ ^ _^ and |^ _^ ^ _^ defined by the same frozen sets ^, ^ ⊂ ^ = {1, … , ^}, have the same stabilizer groups ^ _^ and ^ _^ , modulo a sign factor even though their corresponding states ( ^| ^ _^ ^{^} _^ , ^| ^ _^ ^{^} _^ ) and ( ^| ^ _^ ^{^} _^ , ^| ^ _^ ^{^} _^ ) may be different. Hence, |^^ ^ ^ ^ _^ and |^ _^ ^ _^ are equivalent as stabilizer states. Step E2 is a loop counter indexed by an integer ^, that controls the recursive levels of a loop from ^ = 1 to ^ = ^. The output of step E1, that is, the initial quantum system of single- qubit states may be considered as a preliminary level ^ = 0. At each recursive level ^, a set of 2 ^{^⁄} 2 ^{^} equivalent quantum polar code states ^^^ _^ ^{^} ^ ^{^ ^ ^} _^ ( _^ ), = 1 ^^ 2 ^⁄ 2 ^ of codelengths ^ = 2 ^{^} are prepared. In particular, step E3 concerns the preparation at the recursive level ^, of a set of ^ ^⁄ ^ = 2 ^{^⁄} 2 ^{^} quantum polar code states = 1 ^^ ^ ^⁄ ^ ^ of codelengths ^ = 2 ^{^} . This set of ^ ^⁄ ^ quantum polar code states is referenced by corresponding sets ^ _^(^) of indices comprising first and second sets of frozen indices = are specified for ^ = 1 ^^ ^, as follows: For the last recursion level, ^ = ^, the corresponding index = ^, where ^ ∈ ^ = {1, … , ^} is the selected information index used to encode quantum information, hence ^ _^(^) = ^ and ^ _^(^) = ^ from Eq. (30). For ^ < ^, the value of ∈ ^{ 1, … , ^ = 2 ^{^}} is determined from that of ∈ ^{ 1, … ,2^ = 2 ^{^^^}} , according to the following rule: ^ _^(^) = ^ _^(^^^) − ^, if ^ _{^(^ )} > ^ (33) ^ _^^ ^ _^(^) = ^ _^(^^^) , otherwise Hence, ^ _(^) = ^1 _^(^) , … , ^ _^(^) ^ is the set of the first indices of the set ^ _^(^) , ^ _^(^) is the set of the last 2 ^{^} − ^ _^(^) indices of the set ^ _^(^) . At the recursive level ^, each quantum polar code state ^| ^ _^ ^{^} _^^(^) (i.e. ^^ _^ ^^ _^^(^) ) is prepared by applying two qubit Pauli measurement ^ ⊗ ^ circuits on two equivalent polar code states belonging to the of the antecedent recursive level ^ − 1 and referenced by two corresponding sets of indices . Each set of indices comprises first and second sets of frozen indices ^ _^(^^^) and corresponding to the first and second quantum bases respectively, which example, are the computational and phase bases states respectively. The type of the two-qubit Pauli measurement ^ ⊗ ^ circuits applied on the two equivalent polar code states depends on the last index of the first set of frozen indices ^ _^(^) , as described below in relation to Figs.8A and 8B. Step E4 is a test, such that as long as ^ < ^, the loop counter is incremented ^ = ^ + 1 before looping back to step E3. Otherwise, i.e., when ^ = ^, we go to step E5. At step E5, we have the desired quantum polar code |^ _^ ^{^} ^ _^ of codelength ^ = 2 ^{^} . The above method shows how to prepare polar code states of length ^, given two equivalent polar code states of length ^⁄ 2 , and so on. For ^ = 1, the method is initialized by for example, preparing all the qubits in the computational basis, which can be done for instance by performing single qubit Pauli-^ measurements. Hence, applied recursively, the method is adapted to prepare polar code states of any length ^. Figs. 8A and 8B, schematically represent the preparation of polar code states of length ^, given two equivalent polar code states of length ^/2, according to an embodiment of the present invention. In fact, the represented construction applies to any polar code of length ^ = 2 ^{^} , ^ = 1 ^^ ^, given two equivalent polar code states of length ^/2. Furthermore, for any given ^ = 2 , ^ > 0, two cases are to be considered, according to whether ≤ ^ represented in Figs.8A and 8B respectively. For simplicity of notation, we shall consider hereafter, the preparation of the polar code state of length ^, given the output at the recursive level ^ − 1, of two equivalent polar code states of length ^/2. Hence, since = ^, the two cases to be distinguished are ^ > ^/2 and ^ ≤ ^/2. as represented in Figs.8A and 8B respectively. In particular, the example in Fig.8A schematically represents the preparation of polar code state C1 of length ^ using Pauli ^ ⊗ ^ measurement circuits 71, given two equivalent polar code states C2, C3 of length ^/2, when ^ > ^/2. The shorthand notation from Fig.6A is used to represent Pauli ^ ⊗ ^ measurement circuits. et ^ ^{^} ^ L ≔ ^ − ^{^ ^} ^, and ^′ ≔ {1, … , ^ }, and ^′ ≔ ^{ ^ + 1, … , ^/2 ^} . Assume that the following two equivalent polar code states of length ^/2 are already determined: ^^ ^ ^ _/^ ^ ^{^} _^ ≔ ^ _^/^ ^| ^ _^ , ^ _^ ^{^} _^^ (34) ^^ ^ ^ _/^ ^ ^{^} _^ ≔ ^ _^/^ |^ _^ , ^ _^ ^ _^^ (35) ^ where binary vectors ^ _^ , ^ _^ ∈ {0, 1} ^{^^} and binary vectors ^ _^ , ^ _^ ∈ {0, 1} ^ _^^ ^{^} are known. The method according to the invention corresponds to performing qubit wise Pauli ^ ⊗ ^ measurements on the corresponding qubits in systems and ^ _^ . The measurement outcome ^ is as follows: ^ = ^ _^/^ ⁽ ^′, ^ ⁾ ∈ ^{ 0,1 ^{}^/^} (36) ^ where binary vector ^′ = ^ ⊕ ^ ∈ 0, 1 ^{^^} , and for some ra ^{^} ^ _{^^ ^ ^} ^{{ }} ndom binary vector ^ ∈ ^{ 0,1 ^} . After the measurements, we are left with the following polar code state of length ^ on the joint system where binary vector ^′ = ^ _^ ⊕ ^ _^ . To prove that the polar code state ^| ^ _^ ^{^} as expressed in Eq. (37) is obtained out of the Pauli ^ ⊗ ^ measurements on polar code states we express these polar code states in the computational basis, as follows: where the normalization factors are ignored and where the third equality follows from Eq. (20). On the other hand, the polar code state also be expressed in the computational basis similarly to Eq. (38). As described in relation to Fig.6A (where and ^ _^ are denoted as ^ _^ and ^ _^ ), the qubit wise Pauli ^ ⊗ ^ measurement is done according to the following two steps: At the first step, an ^/2 qubit ancilla state |^^ _^ ≔ |0^ ⊗ ⋯ ⊗ |0^ is defined and then qubit wise CNOT gates ^^ _^^→^ and ^^ _^^→^ are performed. This gives the following joint quantum state Then, at the second step, single qubit Pauli-^ measurement is done on each qubit in the ancilla system ^. The measurement outcome gives ^ = ^ _^/^ (^′, ^), where ^′ = ^ _^ ⊕ ^ ∈ {0, 1} ^{^^} ^ , and for som { } ^{^} ^ _{^^ ^} e random ^ ∈ 0,1 . The joint state of systems and ^ _^ after the measurements is as follows: The above quantum state |^′^ _^^^^ corresponds to the polar code state |^ _^ ^ in Eq.(37), as shown by the following equations: where the third equality follows from the recursion of polar transform given in Fig.3A, and where ^′ = ^ _^ ⊕ ^ _^ in the fourth equality. Finally, ^ can be determined from the measurement outcome. In fact, Eq. (36), implies (^′, ^) = ^ _^/^ (^) (using ^^ ^ = ^), and therefore, ^ is determined as follows: ^ = ^ _^/^ (^)| _^ ^ (42) Hence, for the state ^| ^ _^ ^{^} in Eq. (37), we have ^ ^{^} ^ + ^, and the corresponding frozen states |^^ _^ = |(^ ^{^} , ^, ^ _^ )^ _^ , and |^^ _^ ≔ |^ ^{^} ^ _^ . Fig.8B schematically represents the preparation of polar code state C11 of length ^, given two equivalent polar code states C12, C13 of length ^/2, when ^ ≤ ^/2, according to an example of the present invention. Let assume that we have been given the following equivalent polar code states of length ^/2: where binary vectors ^ _^ , ^ _^ ∈ {0, 1} ^{^^} and ^ _^ , ^ _^ ∈ {0, 1} ^{^/^^^^} are known. In this case, the preparation method corresponds to performing qubit wise Pauli ^ ⊗ ^ measurement circuits 72 on the corresponding qubits in and _^ Similar to the previous case, the outcome of Pauli ^ ⊗ ^ measurements gives: ^ = ^ _^ ^{^} ( _{^, ^′} ) _∈ { _{0, 1} } ^{^/^ (45)} ^ where vector ^′ = ^ _^ ⊕ ^ _^ and ^ ∈ ^{ 0,1 ^{}^^} is a random vector. After the measurements, we get the following polar code state: where ^′ = ^ _^ ⊕ ^ _^ , and ^ is determined from the measurement outcome in Eq. (45), as follows: ^ = ^ _^ ^{^} / _{^ (^)|^^} ⁽⁴⁷⁾ For the state |^ _^ ^ in Eq. (46), we have the corresponding frozen states In the general case of preparing any polar code of length ^ = 2 ^{^} , ^ = 1 ^^ ^, given two equivalent polar code states of length ^/2, index ^ = ^ _^(^) is to be replaced by and index ^ ^{^} = is to be replaced by ^ _^(^^^) . Furthermore, frozen sets ^ = ^ _^(^) , ^ = ^ _^(^) , ^ ^{^} = ^ _^(^^^) and ^ ^{^} = ^ _^(^^^) are to be replaced by ^ _^(^) , ^ _^(^) , ^ _^(^^^) , and ^ _^(^^^) , respectively. The two equivalent polar code states on which Pauli ^ ⊗ ^ or Pauli ^ ⊗ ^ measurement circuits are applied, are to be replaced by two equivalent polar code states Hence, for any level of recursion ^, the quantum polar code is associated the one hand, to the frozen states when the Pauli ^ ⊗ ^ measurement circuits were applied on the corresponding qubits at the recursive level ^ − 1, and on the other hand, to the frozen states when the Pauli ^ ⊗ ^ measurement circuits were applied on the corresponding qubits of the precedent recursive level ^ − 1. The above frozen states are defined such that, ^′ = 2 ^{^^^} , for Pauli ^ ⊗ ^ measurements and with = ^ _^(^) , for Pauli ^ ⊗ ^ The two equivalent polar code states from the precedent level ^ − 1, on which the Pauli ^ ⊗ ^ or Pauli ^ ⊗ ^ measurement circuits were applied, were associated with frozen states respectively. Vectors ^ and ^ are estimated on the basis of the measurement outcome ^ of the corresponding Pauli measurement circuits and corresponding polar transform ^ ^ ^ _/^ or ^ _^/^ , using ^ = when the Pauli ^ ⊗ ^ measurement circuits were applied on the qubits at the recursive level ^ − 1, and ^ = when the Pauli ^ ⊗ ^ circuits were applied on the corresponding qubits of the precedent recursive level ^ − 1. In the particular case ^ = ^ (i.e. ^ = ^), the quantum polar code |^ _^ ^ is associated on the one hand, to the frozen states ^| ^ ^{^} _^ = ^|( ^ ^{^} , ^, ^ _^ ^{)^} _^ , and ^| ^ ^{^} _^ ≔ ^| ^ ^{^^} _^ when the Pauli ^ ⊗ ^ measurement circuits were applied on the corresponding qubits of the precedent recursive level ^ = ^ − 1, and on the other hand, to the frozen states ^| ^ ^{^} _^ = ^| ^ ^{^^} _^ , and ^| ^ ^{^} _^ = ^ ^{^(^} ^ ^{^^} _^ ^{^} , ^{^} ^ ^{^^} , ^{^} ^ ^{^^ ^} _)^ ^ when the Pauli ^ ⊗ ^ measurement circuits were applied on the corresponding qubits of the precedent recursive level ^ = ^ − 1. In this case, ^′ = ^ _^ ⊕ ^ _^ ∈ {0, 1} ^{^^} , ^′ = ^ _^ ⊕ ^ _^ , and ^ ^{^} = ^ − ^/2 for Pauli ^ ⊗ ^ measurements and with ^ ^{^} = ^, for Pauli ^ ⊗ ^ measurements. Vectors ^ and ^ are estimated on the basis of the measurement outcome ^ of the corresponding Pauli measurement circuits according to Eqs. (43) and (48), respectively. In the following, we restrict again the discussion to the case ^ = ^, thus ^ = ^, but again, this is only for simplicity of notation and consistency with the previous detailed description of Figs.8A and 8B. As described above, the binary vectors ^ and ^ are estimated on the basis of the measurement outcome ^ of the corresponding Pauli measurement circuits, and corresponding polar transform, either ^ ^ ^ _/^ or ^ _^/^ , according to Eqs. (42) and (47), respectively. In fact, the fault tolerance according to the method of the present invention comes down to estimating the binary vectors ^ and ^ in Eqs. (42) and (47), respectively, from the measurement outcomes of Pauli ^ ⊗ ^, or Pauli ^ ⊗ ^ measurements. In what follows, for a binary vector ^ = (^ _^ , … , ^ _^ ) ∈ ^{ 0, 1 ^{}^} , we define sup (^) ≔ {^ ∣ ^ _^ = 1, ^ ∈ ^{ 1, … , ^ ^} }, and wt(^) ≔ ^| sup ⁽ ^ ^)| . Moreover, for a set ℐ ⊂ {1, … , ^}, we define ^| _ℐ ∈ to be the part of ^ corresponding to the indices in ℐ. To consider the preparation procedure of polar code states, under the effect of noise, we suppose that noise acts on each qubit independently, and it causes random Pauli errors on each qubit, with some probability. The preparation method is said to be fault tolerant if for any ^, the errors on quantum states not propagate to an increased number qubits in |^ _^ ^. In the case ^ > ^/2, instead of Eqs. (43) and (44), we have the following error corrupted polar code states of length ^/2: where binary ^ , ^ _^ ∈ ^{ 0, 1 ^{}^} and ^ _^ , ^ _^ ∈ ^{ 0, 1 ^{}^/^^^^} are known, ^ binary vector errors ^ _^^ , ^ _^^ , ^ _^^ , ^ _^^ ∈ ^{ _{0, 1} ^{} ^} are unknown. For simplicity, we start by assuming that Pauli ^ ⊗ ^ measurements can be performed perfectly, without errors. The qubit wise Pauli ^ ⊗ ^ measurements on the corresponding qubits of and ^ _^ of Eqs. (48) and (49), give the following measurement ^: ^ ^ = ^ _^/^ ⁽ ^′, ^ ⁾ ⊕ ^ _^ ∈ ^{ 0, 1 ^} ^ (50) ^ where ^′ = ^ _^^ ^{^} ^ ⊕ ^ _^ is known, and ^ ∈ ^{ 0,1 ^} ^ , ^ _^ ∶= ^ _^^ ⊕ are unknown. After the measurements, the error corrupted polar code state on ^ qubits, depends on vector ^ according to the following expression: The ^ (respectively, ^) error on ^| ^ _^ ^{^} _^^^^ is simply equal to the ^ (respectively, ^) errors on and thus, the method of preparation according to the present invention is fault tolerant. However, in order to know the prepared polar code state, the vector ^ in Eq. (51) has to be estimated. Note that, in case no error occurs, that is ^ _^ = 0 in Eq. (50), then the vector ^ can be determined as in Eq. (42). However, if an unknown error ^ _^ ≠ 0 occurs in Eq. (50), the equality in Eq. (42) for ^ does not hold anymore. To estimate the vector ^, two techniques can be used, one based on error detection, and the other one based on error correction. The error-detection based technique exploits the redundancy of measurements in the described preparation procedure. As described above in relation to Fig. 8A, let ^′ ≔ {1, … , ^ and ^′ ≔ {^′ + 1, … , ^}. After getting the measurement outcome ^ in Eq. (50), it is verified whether the equality ^ _^/^ ⁽ ^ ⁾ | _^ ^{^} = ^ ^{^} holds or not. If ^ _^/^ ⁽ ^ ⁾ | _^ ^{^} = ^ ^{^} , then the syndrome of the error ^ _^ is zero, that is ^ _^/^ ⁽ ^ _^ ⁾ | _^ ^{^} = ^. Hence, it is assumed that no error occurred, and the value of vector ^ is estimated according to Eq. (42), that is ^ = ^ _^/^ (^)| _^ ^. If ^ _^/^ (^)| _^ ^ ≠ ^ ^{^} , then it is deduced that some error has occurred. Hence, the procedure is discarded and restarted, by taking fresh polar code states of length ^/2. In the general case, i.e., at any level ^, when the Pauli measurement circuits were applied on the corresponding qubits at the recursive level ^ − 1, it is verified whether the equality ^ _^/^ ⁽ ^ ⁾ | _^^(^^^) = holds or not. If ^{^} = ^ , then the vector ^ is estimated using ^ = ^ _^/^ (^)| _^^(^^^) , or otherwise it is discarded and the procedure is restarted by taking fresh polar code states of length ^/2. The error-correction based technique also exploits the redundancy of measurements in the described preparation procedure, but it determines the vector ^ using a maximum likelihood decision rule as follows. As described above in relation to Fig.8A, and going back to the specific case ^ = ^ (for ^ the simplicity and consistency of notation), let ^′ ≔ {1, … , ^′}, and ^′ ≔ {^′ + 1, … , ^}. Further, for an ^ ^{^^} such that ^′ < ^ ^{^^} ^ ≤ , we define ^′ ^{^} ≔ {^′ + 1, ⋯ , ^ ^{^^} } ^{^^ ^ ^} ^ ⊆ ^′, and ^ _^ ≔ ^ ∖ ^′ . Hence, the binary vector ^ can be defined as follows: The ^| _^^ ^{^^} part of ^ in Eq. (52) is estimated from the measurement outcome ^ in Eq. (50), as follows: ^| _^^ ^^ = ^ _{^ ^ |} ^^ (53) ^ ^{( )} _^^ A maximum likelihood decision for the ^| _^ ^{^^} part of ^ in Eq. (52), is then determined as follows: When ^| _^ ^{^^} is not unique, meaning that there are several values ^| _^ ^{^^} giving the minimum weight ^ ⊕ ^^ ^^, the procedure is discarded and restarted again, by ^ taking fresh polar code states of length ^/2. Finally, the above assumption that Pauli ^ ⊗ ^ measurements are perfect is not needed and may be relaxed. In fact, when Pauli ^ ⊗ ^ measurements are noisy, an extra error ^^ ^ on the measurement outcome ^ is introduced in Eq. (50), which becomes ^ = ^ _^/^ (^′, ^) ⊕ ^ _^ ⊕ ^^ ^ . However, ^| ^ _^ ^{^} _^^^^ in Eq. (51) remains the same, except possibly an extra error introduced upon it by the noisy measurement. In the case ^ ≤ ^/2, instead of Eqs. (43) and (44), we have the following error corrupted polar code states of length ^/2: where ^ _^ , ^ _^ ∈ {0, 1} ^{^^} and ^ _^ , ^ _^ ∈ {0, 1} ^{^/^^^^} are known, with ^ ^{^} = ^, and the errors unknown. The qubit wise Pauli ^ ⊗ ^ measurements on the corresponding qubits of Eqs. (55) and (56) give the following: ^ where ^′ = ^ _^ ⊕ ^ _^ ∈ {0,1} ^ _^^ ^{^} ^ ^ is known, and ^ ∈ {0,1} ^{^} , ^ _^ = ^ _^^ ⊕ ^ _^^ ∈ {0, 1} ^{^} are unknown. After the measurements, the error corrupted polar code state on ^ qubits depends on vector ^ according to the following expression: Similarly, in order to know the prepared polar code state, the vector ^ has to be estimated. To estimate the vector ^, two methods can be used, one based on error detection, and the other one based on error correction. The error-detection based method exploits the redundancy of measurements in the described preparation procedure. As described above in relation to Fig. 8B, let ^′ ≔ {1, … , ^′}, ^ and ^′ ≔ {^′ + 1, … , ^}. After getting the measurement outcome ^ in Eq. (57), it is verified if the equality ^ _^ ^{^} / _^ (^)| _^ ^ = ^ ^{^} holds. If ^ _^ ^{^} / _^ ⁽ ^ ⁾ | _^ ^{^} = ^ ^{^} , then the syndrome of the error ^ _^ is zero, that is ^ _^ ^{^} / _^ ⁽ ^ _^ ⁾ | _^ ^{^} = ^. Hence, it is assumed that no error occurred, and the value of vector ^ is estimated according to Eq. (47), that is ^ = ^ _^ ^{^} / _^ ⁽ ^ ⁾ | _^ ^{^} . If ^ _^ ^{^} / _^ (^)| _^ ^ ≠ ^ ^{^} , then it is deduced that some error has occurred. Hence, the procedure is discarded and restarted, by taking fresh polar code states of length ^/2. In the general case, i.e., at any level ^, when the Pauli ^ ⊗ ^ measurement circuits were applied on the corresponding qubits at the recursive level ^ − 1, it is verified whether the equality ^ _^ ^{^} / _^ ⁽ ^ ⁾ | _^^(^^^) = ^ ^{^} holds or not. If ^ _^ ^{^} / _^ ⁽ ^ ⁾ | _^^(^^^) = ^ ^{^} , then the vector ^ is estimated using ^ = ^ _^ ^{^} / _^ ⁽ ^ ⁾ | _^^(^^^) , or otherwise it is discarded and the procedure is restarted by taking fresh polar code states of length ^/2. The error-correction based method also exploits the redundancy of measurements in the described preparation procedure, but it determines the vector ^ using a maximum likelihood decision rule. As described above with relation to Fig. 8B, let ^′ ≔ {1, … , ^ ^{^} }, and + ^ 1, … , ^{^^ ^^ ^ ^^} ^^. Further for an ^ , such that 1 ≤ ^ ≤ ^ + 1, let ^ ≔ {^ ^{^^} , … , ^ ^{^} } ⊆ ^′, and ^ _^ ^{^^} ≔ ^ ^{^} ∖ ^′ ^{^} . Hence, the binary vector ^ can be defined as follows: The ^| _^^ ^{^^} part of ^ in Eq. (59), is estimated from the measurement outcome ^ in Eq. (57), as follows: ^| _^^ ^{^^} = ^ _^ ^{^} / _^ ⁽ _^ ⁾ _|^^ ^^ (60) A maximum likelihood decision for the ^| _^ ^{^^} part of ^ in Eq. (59), is then determined as follows: When ^| _^ ^{^^} is not unique, meaning that there are several ^| _^ ^{^^} giving the minimum weight wt(^ ⊕ ^ _^/^ ^^| _^^ ^^, ^| _^ ^^ , ^^), the procedure is discarded and restarted again, by taking fresh polar code states of length ^/2. In the general case, i.e., at any level ^, the fault tolerant preparation of quantum polar code states is represented by the estimation of vectors ^ and ^. The binary vector ^ is defined as defined as ^ = ^ ^| _^^ ^{^^} , ^| _^ ^{^^} ^, with respect to the sets ^ _^(^^^) = {1, … , ^ _^(^^^) }, ^′′ ≔ Figs. 9A, 9B, 9C and 9D schematically represent different examples of fault tolerant preparation of quantum polar code states, according to an embodiment of the present invention. The four examples concern fault tolerant preparation of polar code states using Pauli ^ ⊗ ^ and ^ ⊗ ^ measurements, recursively. For these examples, the sets ^′ and ^′, are specified at each level of recursion. To ensure fault-tolerance, one of the two techniques described above can be used, based on either error detection or error correction. To implement the error detection based technique, only the sets ^′ and ^′, and the measurement outcome ^, have to be specified. Further, for the error-correction based technique, when qubit wise Pauli ^ ⊗ ^ measurements are performed, the set ^′ ^{^} (hence, also ^ _^ ^{^^} ) is specified, and thus, Eqs. (53) and (54) are used to estimate vector ^. When qubit wise Pauli ^ ⊗ ^ measurements are performed, the set ^′ ^{^} (hence, also ^ _^ ^{^^} ) is specified, and thus, Eqs. (60) and (61) are used to estimate vector ^. In the description of Figs. 9A, 9B, 9C and 9D, we consider the example of error- correction based technique. In particular, Fig.9A concerns the recursive preparation of a quantum polar code ^| ^ _^ ^{^} of length ^ = 8 with ^ = {1,2,3} and ^ = {4,5,6,7,8}. A single qubit Pauli ^ measurement is performed on each qubit such to initialize them into a computational basis state. Alternatively, any other method to initialize the qubits in the computational basis can be used. Then, at each level of recursion, polar code states of lengths ^ = 2, 4, 8 are prepared, as follows: At the first level of recursion, four polar code states ^ ^{^} ^ _^ ^{^^} _, ^{^} _^ ^ ^ ^^, ^{^} ^ _^ ^, ^| ^ _^ ^{^^} ^ of length ^ = 2, with ^ = {1} and ^ = {2} are prepared represented by the four blocks B2. Here, ^ = {1} and ^ = ^{ 2 ^} refer to the first and second qubit of each of the four blocks B2. Here ^ ^{^} = {1}, ^′ = ∅, and ^ ^{^^} = ∅, hence ^ _^ ^{^^} = ^ ^{^} . The estimate of ^ ∈ {0,1} is determined from Eq. (60), as follows: At the second level of recursion, two polar code states of length ^ = 4, with ^ = ^{ 1,2,3 ^} and ^ = ^{ 4 ^} are prepared, represented by the two blocks B4. Here, ^ = ^{ 1,2,3 ^} and ^ = ^{ 4 ^} refer to the first, second, third, and fourth qubit of each of the two blocks B4. Then, ^′ = ^{ 1 ^} , ^′ = {2}, and ^′ ^{^} = ∅, hence ^ _^ ^{^^} = ^ ^{^} . The estimate of ^ ∈ {0,1} is determined from Eq. (53), as follows: ^ = ^ _^ ⁽ ^ ⁾ | _^ ^{^} = ^ _^ ⊕ ^ _^ , where ^ = (^ _^ , ^ _^ ) ∈ ^{ 0,1 ^{}^} (63) Finally, at the third level of recursion, the desired polar code state |^ _^ ^ of length ^ = 8, with ^ = {1,2,3} and ^ = {4,5,6,7,8} is prepared, represented by the single blocks B8. Here ^ ^{^} = {1,2,3}, ^ ^{^} = {4}, and ^ ^{^^} = ∅, hence ^ _^ ^{^^} = ^ ^{^} . The estimate of ^ ∈ {0,1} ^{^} is determined from Eq. (60), as follows: where ^ = (^ _^ , ^ _^ , ^ _^ , ^ _^ ) ∈ {0,1} ^{^} . Fig.9B concerns the recursive preparation of a quantum polar code ^| ^ _^^ ^{^} of length ^ = 16, with ^ = ^{ 1, … , 7 ^} and ^ = ^{ 8, … ,16 ^} . The first two levels of recursion are the same as the first two level of recursions in Fig. 9A, except that there are twice as many polar code states of length ^ = 2 represented by the eight blocks B12, and twice as many polar code states of length ^ = 4 represented by the four blocks B14. Thus, only the third and fourth levels of recursion will be described hereafter. At the third level of recursion, two polar code states of length ^ = 8, with ^ = ^{ 1, … , 7 ^} and ^ = ^{ 8 ^} are prepared, represented by the two blocks B18. Here, ^′ = ^{ 1,2,3 ^} and ^′ = ^{ 4 ^} , then, ^ ^{^^} = ^ ^{^} = {4}, hence ^ _^ ^{^^} = ∅. The estimate of ^ ∈ {0,1} is then determined from Eq. (54), as follows: ^ = argmin _^∈{^,^} wt (^ ⊕ ^ _^ ⁽ ^′, ^ ⁾ ) (65) Recall that if ^ is not unique, the prepared state of length ^ = 8 is discarded without proceeding to the next level of recursion and the whole procedure is restarted back from the initialization. Finally, at the fourth level of recursion, the desired polar code state ^| ^ _^^ ^{^} of length ^ = 16, with ^ = ^{ 1, … , 7 ^} and ^ = ^{ 8, … ,16 ^} is prepared, represented by the single block B116. We have ^ ^{^} = {1, … , 7} , ^ ^{^} = {8}, and ^ ^{^^} = ∅, hence ^ _^ ^{^^} = ^ ^{^} . The estimate of ^ ∈ {0,1} ^{^} is then determined from Eq. (60), as follows: where ^ = ⁽ ^ _^ , … , ^ _^ ⁾ ∈ ^{ 0,1 ^{}^} . Fig.9C concerns the recursive preparation of a quantum polar code ^| ^ _^^ ^{^} of length ^ = 16, with ^ = {1, … , 11} and ^ = {12, … ,16}. At the first two levels of recursion, eight polar code states of length ^ = 2 and four polar code states of length ^ = 4 are prepared, represented by eight blocks B22 and four blocks B24, respectively. At the third level of recursion, two instances of polar code states of length ^ = 8, represented by two blocks B28, with ^ = ^{ 1,2,3 ^} and ^ = ^{ 4,5,6,7,8 ^} , are prepared. This third level is similar to the third level of the first example depicted in Fig. 9A which was already described. Hereafter, we describe the fourth level of recursion. At the fourth level of recursion, the desired polar code state ^| ^ _^^ ^{^} of length ^ = 16, with ^ = {1, … , 11} and ^ = {12, … ,16} is prepared, represented by a single block B216. Here, ^′ = {1,2,3} and ^′ = {4,5,6,7,8}. We take ^ ^{^^} = {4}, and hence ^ _^ ^{^^} = {5,6,7,8}. From Eq. (53), the ^| _^^ ^{^^} ∈ ^{ 0,1 ^{}^} part of ^ is estimated, as follows: From Eq. (54), the ^| _^ ^^ ∈ {0,1} part of ^ is estimated, as follows: If ^| _^ ^{^^} is not unique, the prepared state of length ^ = 16 is discarded and the whole procedure is restarted back from the initialization. Fig.9D concerns the recursive preparation of a quantum polar code ^| ^ _^^ ^{^} of length ^ = 64, with ^ = {1, … , 23} and ^ = {24, … ,64}. Here, after the fourth level of recursion, we have four blocks B316 of polar code states of length ^ = 16, with ^ = {1, … , 7} and ^ = {8, … ,16}, as in the second example described in relation to Fig.9B. Below, only the remaining two levels of recursion are described. At the fifth level of recursion, two polar code states of length ^ = 32, with ^ = ^{ 1, … , 23 ^} and ^ = ^{ 24, … ,32 ^} , represented by the two blocks B332 are prepared. Here, we have ^′ = ^{ 1, … , 7 ^} and ^′ = ^{ 8, … ,16 ^} . We take ^ ^{^^} = ^{ 8 ^} , and hence ^ _^ ^{^^} = ^{ 9, … , 16 ^} . From Eq. (53), the ^| _^^ ^ ∈ {0,1} ^{^} part of ^ is estimated, that is, ^| _^^ ^^ = ^ _^^ (^)| _^^ ^^, where ^ ∈ {0,1} ^{^^} . Further, from Eq. (54), we get the estimate ^| _^ ^^ = argmin _^∈{^,^} wt (^ ⊕ At the sixth level of recursion, the desired polar code state |^ _^^ ^ of length ^ = 64, with ^ = {1, … , 23} and ^ = {24, … ,64} is prepared, represented by the ingle block B364. Here, we have ^′ = ^{ 1, … , 23 ^} and ^′ = ^{ 24, … ,32 ^} . We take ^ ^{^^} = ∅, and hence ^ _^ ^{^^} = ^′. From Eq. (60), we get the estimate of ^ ∈ {0,1} ^{^^} , that is, ^ = where ^ ∈ ^{ 0,1 ^{}^^} . The fault tolerant preparation of polar code states according to the embodiments of the present invention is an important resource for exploiting quantum polar codes for fault tolerant systems of quantum computation. Fig. 10 schematically represents a computing system according to a preferred embodiment of the present invention. The computing system 11 comprises a classical computing system 13, a classical- quantum interface 15 and a quantum computing system 1 configured to prepare quantum polar codes according to the above embodiments of the present invention. The quantum computing system 1 is coupled to the classical computing system 13 via the classical-quantum interface 15. The classical-quantum interface 15 comprises a syndrome extractor 17 configured to extract a syndrome out of quantum measurements implemented by the quantum computing system 1. The classical computing system 13 comprises a classical decoder 19 configured to decode the quantum polar codes prepared according to the above methods by implementing successive cancelation decoding.