here is my code;
% 1) Compute the volume of a cylindrical shell with an inner radius
% r1=5,
% an outer radius
% r2=6, and length=10.
r1 = 5; %inner radius
r2 = 6; %outer radius
h = 10; %height
dr = 0.01; dz = 0.1; dphi = 0.01; % set vector of iteration
r = [r1:dr:r2];z = [0:dz:h];phi = [0:dphi:2*pi];noi_r = length(r); noi_phi = length(phi); noi_z = length(z);x = zeros(noi_r,noi_phi,noi_z);y = zeros(noi_r,noi_phi,noi_z);z = zeros(noi_r,noi_phi,noi_z);plot3(NaN,NaN,NaN,'.')hold on%initial volume
Volume =0; for ir = 1:noi_r for iphi = 1:noi_phi for iz = 1:noi_z dv = r(ir)*dphi*dz*dr; Volume = Volume +dv; x(ir,iphi,iz) = r(ir)*cos(phi(iphi)); y(ir,iphi,iz) = r(ir)*sin(phi(iphi)); z(ir,iphi,iz) = z(iz); plot3(x,y,z,'.') end end endhold offdisp('volume is = ')disp(Volume)
Best Answer