for E design</span>
W = V*diag(lambdaE)*V';
angles = linspace(0,2*pi,noangles);
R = chol(W);  <span class="comment">% W = R'*R</span>
ellipsoid = sqrt(eta)*(R\[cos(angles); sin(angles)]);
d=plot(0,0,<span class="string">'ok'</span>,ellipsoid(1,:), ellipsoid(2,:), <span class="string">'-'</span>);
set(d,<span class="string">'Color'</span>,[0 0.5 0]);
text(ellipsoid(1,4000),ellipsoid(2,4000),<span class="string">'E'</span>);

<span class="comment">% draw 90 percent confidence ellipsoid  for uniform design</sp