Area and volume properties : v=splder(td,zd,t) - File Exchange

Code covered by the BSD License

Highlights from
Area and volume properties

[Area,Ax,Ay,Ixx,Iyy,Pxy]=spln...
[Area,xcent,ycent,Axx,Ayy,Axy... [Area,xcent,ycent,Axx,Ayy,Axy]=areaprop(pdat,sdat,titl,doplot)
[Area,xcent,ycent,Ixx,Iyy,Ixy... [Area,xcent,ycent,Ixx,Iyy,Ixy]=ellipsecheck(n) checks
[V,Rc,VRR,I]=volprop(pdat,sda...
[pdat,sdat]=areabook [pdat,sdat]=areabook creates data for
[pdat,sdat]=elipwholes [pdat,sdat]=elipwholes creates an ellipse with
[pdat,sdat]=torconedat [pdat,sdat]=torconedat creates data for
[pdat,sdat]=triacirc This function defines a semicircle on the
[pdat,sdat]=volbook [pdat,sdat]=volbook example from
ellipticbowl Example 8: Elliptic Bowl with a Semi-circular Rim
examples Ten examples are given which show use
plotdat(pdat,sdat,titl) plotdat(pdat,sdat,titl) plots a geometry
plotvol(pdat,sdat,t1,t2,titl,... plotvol(pdat,sdat,t1,t2,titl,nrot,nspl)
readme Computation of geometrical properties of plane areas including
u=chkdat(v)
v=cplx2mat(z)
v=splder(td,zd,t) v=splder(td,zd,t) cubic spline differentiation
v=xymom(type,xd,yd,n,m)
View all files

from Area and volume properties by Howard Wilson
Area and volume properties for shapes defined by piecewise linear and spline boundaries

v=splder(td,zd,t)

function v=splder(td,zd,t)
% v=splder(td,zd,t) cubic spline differentiation
%    HBW, 9/6/09
n=length(td);
if n>3
  [b,c]=unmkpp(spline(td,zd));
  c=[3*c(:,1),2*c(:,2),c(:,3)];
  v=ppval(mkpp(b,c),t);
elseif n==3
  c=[[1;1;1],td(:),td(:).^2]\zd(:);
  v=c(2)+2*c(3)*t;    
elseif n==2   
  v=(zd(2)-zd(1))/(td(2)-td(1))*ones(size(t));   
else
 v=nan;  
 error('More than 1 data point is required.')
end

Highlights from Area and volume properties

Highlights from
Area and volume properties