% Return computational basis states |0>,|1>.

function [z0,z1] = qstates()
z0 = [1;0];
z1 = [0;1];
end