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:(adsbygoogle = window.adsbygoogle || []).push({});

% 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.

**Physics Forums | Science Articles, Homework Help, Discussion**

The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

# MATLAB Model a circle using finite difference equation in matlab

**Physics Forums | Science Articles, Homework Help, Discussion**