Finally I got my own answer.
a=ones(500,500); %matrix size
I= 0.5 %Intesity inside the cirlce
for i=1:500
for j=1:500
if ((j-x)^2+(i-y)^2)<=R^2 %condition to make the cirlce
a(i,j)=I;
end
end
end
imshow(a)