--------------------------------------------
%% Coordinate
l=1;
T=5;
Nx=100;
Ny=100;
u=linspace(-T/2,T/2,Nx);
v=linspace(-T/2,T/2,Ny);
[x,y]=meshgrid(u,v);
r=(x.^2+y.^2).^(1/2);
phi=angle(x+1i*y);

%% Calculation
E1=exp(-(r.^2)).*exp(-1i*phi*l).*(2*r.^2).^(l/2);
E2=exp(-1i*2*pi*4*x);
Etot=E1+E2;
I=abs(Etot).^2;

% Output
save EH;
imagesc(I);
colormap(gray);
-------------------------------------------
大家好，我是Matlab新手

以上是我的Matlab code

如果我要加入截取如下

Ｉ>=0.5 and r<=1, 則 I=1
 else                            I=0 

請問要如何加入與怎麼寫這程式呢？

我有試著寫條件指令，可是都失敗

請大家教導了。謝謝大家！