function t=torus(a, b)
v=(1/4)*(pi^2)*(a+b)*((b-a)^2);
c=(pi^2)*(b^2-a^2);
%put v and c into a vector
t=[v,c]
end