probIj = 0.5; varIj = probIj*(1-probIj); var = 10*varIj; disp(var, "Thus, the required variance is")