bonjour tout le monde. je dois créer un cercle avec des bords. bref, dans un carré de dimension 0.61m*0.61m, j'aimerai :
-un cercle centré en zero de rayon 0.56m (donc rien à faire pour cette partie)
-diviser ce carré en carrés de dimension 0.05m*0.05m
-pour tous les carrés, vérifier si oui ou non leur centre est à l'interieur du cercle, si oui ne rien faire, si non, ecrire dans un fichier .dat "#condition = xmin xmax ymin ymax" avec x y correspondant aux coins des petits carrés. cad xmax=xmin+0.05 et ymax=ymin+0.05
 
voila ou jen suis pour le moment:
 
%oxysol ertlab creation cellule
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 
clear all;
close all;
clc;
 
% DESCRIPTION Cercle
 
title('oxysol')
xlabel('m')
ylabel('m')
 
% Cercle
disp('Veuillez rentrer la valeur du rayon R:')
R= input('Valeur du rayon R: ')
disp('Veuillez rentrer l''abscisse du centre du cercle:')
Xc= input('Abscisse :')
disp('Veuillez rentrer l''ordonnée du centre du cercle:')
Yc= input('Ordonnée:')
VThetaDeg = 0:1:360;
VTheta = VThetaDeg *pi / 180;
XCercle = Xc + R * cos(VTheta);
YCercle = Yc + R * sin(VTheta);
%figure
plot(XCercle, YCercle)
axis([Xc-R-1,Xc+R+1,Yc-R-1,Yc+R+1])
hold on
grid on
 
% Carrée
for i=-Xc-R-1:Xc+R+1
   for j=-Yc-R-1:Yc+R+1
rectangle('Position',[i,j,1,1],...
          'Curvature',[0,0],...
         'LineWidth',1,'LineStyle','--')
daspect([1,1,1]);grid on
plot(i+1/2,j+1/2,'rx')
   end
end
 
 
 
merci pour votre aide