s = [ 0.5*cisd(-120) 0.08*cisd(-26); 0.08*cisd(-26) 0.99*cisd(172) ]