Ooh the joys of matlab!
lam11=[];
lam22=[];
x11=[];
x22=[];
for i=[100:350]
q=3+i/(100/(11/3-3));
r1=-58.*(219.*(q-11/3)-(732381.*(q-11/3).^2+713168.*(q-11/3)).^(1./2))./(178292+171105.*(q-11/3))./cos(theta);
alp1 = 18.*(4.*(q-11/3)+5.*r1.^2.*cos(theta).^2+30.*(q-11/3).*r1.^2.*cos(theta).^2)./(116+45.*r1.*cos(theta)+690.*r1.^2.*cos(theta).^2);
x1= r1.*cos(theta)+(1-1./(alp1+3));
lam1=alp1+3;
r2=-58.*(219.*(q-11/3)+(732381.*(q-11/3).^2+713168.*(q-11/3)).^(1./2))./(178292+171105.*(q-11/3))./cos(theta);
alp2 = 18.*(4.*(q-11/3)+5.*r2.^2.*cos(theta).^2+30.*(q-11/3).*r2.^2.*cos(theta).^2)./(116+45.*r2.*cos(theta)+690.*r2.^2.*cos(theta).^2);
x2 = r2.*cos(theta)+(1-1./(alp2+3));
lam2=alp2+3;
x11=union(x11,x1);
lam11=union(lam11,lam1);
x22=union(x22,x2);
lam22=union(lam22,lam2);
end
x22=sort(x22,'descend');
plot(lam11,x11)
hold on
plot(lam22,x22)