# MATLAB Model a circle using finite difference equation in matlab

#### tomallan

hello. I have a matlab skeleton provided because i want to model a distribution with a circular geometry. all in all, i want the 3d graph of the code to be some type of cylinder. This is the code:

% flat step condition
for ii=1:nHi,
for jj=1:nHj,
if (X(ii)/R_P)<1 & (Y(ii)/R_P)<1,
h(ii,jj)=h_init+h_step;
else,
h(ii,jj)=h_init;
end
end
end

nHi is the maximum size of nodes in X direction and nHj is the maximum size in Y direction while R_P is the radius of the circle. I am struggling with this code because when I execute this, it only appears with a rectangular graph, instead of a circle-ish, when viewed from the top.

what the code is trying to say is,
if node is less than the diameter (x or y node),
thickness is initial+step
else,
thickness is initial.

Any tips will be greatly appreciated.

Related MATLAB, Maple, Mathematica, LaTeX, Etc News on Phys.org

#### da_nang

Just a hunch, have you tried replacing this:

if (X(ii)/R_P)<1 & (Y(ii)/R_P)<1

with this?

if ((X(ii) - Xm)^2 + (Y(ii) - Ym)^2 < R_P^2)

where Xm and Ym are the coordinates of the center of the circle.

#### Strum

Are you trying to embed a circle to a square domain or are you trying to make a circle domain? Why not do it in polar coordinates?

### Physics Forums Values

We Value Quality
• Topics based on mainstream science
• Proper English grammar and spelling
We Value Civility
• Positive and compassionate attitudes
• Patience while debating
We Value Productivity
• Disciplined to remain on-topic
• Recognition of own weaknesses
• Solo and co-op problem solving