% Test run for parser

Dirac

D: |a> |b> |c> |d> (c1) (c2) [O1] [O2] [U][sig2]
M: a=[1;0]; b=[0;1];
c1=0.7071
c2=0.7071
|d> = c1 |a> - c2|b>
|c> = c1 |a> + c2|b>
|a>
c1=0.7071
c2=0.7071
c2(|a>)
c2 ( c1 )
O1 = |a><b|
c2 O1
c2 (O1) (|a>)
c2(|a>)
c2 (O1)
c2 O1
c1 |a><a|
O1
exp(i O1)
M: sig2=[0 1;1 0]
M: U=eye(2)
U cos(1) + i sig2 sin(1) - exp(i sig2)
c1 O2
c1(O1)
c1 O1 + c1 O1
O2=O1
c1( O1 + O2 )
c1(O1+O2)
ss=c1(O1+O2)
c1(|a>)
c1(c2|a>)
c1(O1|a>)
sin(1) |a>
sin(1)
c1(c2)
sin(c1)
sin(c1)
c1(c1)
c1(|a>)
c1(|a>)+ sin(2)|a>
sin(2)|a>+c1(|a>)
sin(2)|a>+c1(|a>)
c1(|a>)+ sin(2)|a>
s=c1(|a>)+ sin(2)|a>
Q=|s><s|
exp(i Q)
s=c1(|a>)+ sin(2)|a>
sin(1) O1
(|c><c|)| a>
|c><c| a>
( |c><c| )| a>
(|c><c|)| a>
|d> = c1 |a> - c2|b>
|c> = c2 |a> + c1|b>
(|c><c|)| a>
sin(1) O1
s=c1(|a>)+ sin(2)|a>
(|c><c|)| a>
|c><c| a>
<c|c>
(|c><c|)| a>
 |c>=(|c><c|)|b>
|c><c|
<c|a>
O1=|c><c|
O2=|d><d|
O2|a>
(O1+O2)|a>
M: U=eye(2)
(O1+O2)|a>
(O1+O2)|a>
c1 (O1+O2)|a>
(c1)(O1+O2)|a>
c1 (O1+O2)|a>
c1(c2)
<c|= c1 <a| O2 O2
c1(O1+O2)|a>
U=(O1+O2)|a><a|
U=(|a><a|)
%%%  crashes
% <c|=  c1 c2 <a|  
 ( |a><a| + |a><a| )  |c>
   |c> =     ( |a><a|) |c>
|c> = ( c1 |a><a| + c2|a><a| ) |c>