function f = achromat2(v) % achromatic lens: adjust shape factors global cv th rn dn; pwr = get_powers; cv(2:3)=lens_shape(pwr(1),v(1),rn(2)); cv(4:5)=lens_shape(pwr(2),v(2),rn(4)); yap = 50/8; uco = tan(2.0*pi/180); ya = parax([yap 0],cv,th,rn); yc = parax([0 uco],cv,th,rn); w = ford(ya,yc,cv,th,rn,0.006); f(1) = w(1); % spherical aberration f(2) = w(2); % coma f = f(:);