% truncated cone diagram % Modified for PGF/TikZ def O (0,0,0) def I [1,0,0] def J [0,1,0] def K [0,0,1] def p0 (1,2) def p1 (1.5,0) def N 8 def seg_rot rotate(360 / N, [J]) def dx 2 <> 2.3 def dy 2 <> 3.3 def dz dx def basic_cone { % draw the cone; this is the easy part! sweep[cull=false,fill=blue!20,fill opacity=0.8] { N, [[seg_rot]] } line(p0)(p1) % draw the axes def ax (dx,0,0) def ay (0,dy,0) def az (0,0,dz) line[arrows=<->,linewidth=.4pt](ax)(O)(ay) line[arrows=->,linewidth=.4pt](O)(az) % repeat dotted as an overlay to hint at the hidden lines %line[lay=over,linestyle=dotted,linewidth=.4pt](ax)(O)(ay) % line[lay=over,linestyle=dotted,linewidth=.4pt](O)(az) % label special|\path #1 node[below] {$x$} #2 node[above] {$y$} #3 node[left] {$z$};| (ax)(ay)(az) % height measurement mark takes too much code! def c0 (p0) then scale([J]) def hdim_ref unit((p1) - (O)) then [[seg_rot]]^2 def h00 (c0) + 1.1 * [hdim_ref] def h01 (c0) + 1.9 * [hdim_ref] def h02 (c0) + 1.8 * [hdim_ref] line(h00)(h01) def h10 (O) + 1.6 * [hdim_ref] def h11 (O) + 1.9 * [hdim_ref] def h12 (O) + 1.8 * [hdim_ref] line(h10)(h11) line[arrows=<->](h02)(h12) def hm2 ((h02)-(O)+(h12)-(O)) / 2 + (O) special|\node[ann] at #1 {$h$};|(hm2) % radius measurement marks def gap [0,.2,0] % first r1 def up1 [0,3.1,0] def r1 ((p1) then [[seg_rot]]^-2) + [up1] def r1c (r1) then scale([J]) def r1t (r1) + [gap] def r1b ((r1t) then scale([1,0,1])) + [gap] line[arrows=<->](r1c)(r1) line(r1b)(r1t) def r1m ((r1) - (O) + (r1c) - (O)) / 2 + (O) special |\node[ann] at #1 {$r_1$};|(r1m) % same drill for r0, but must project down first def up0 [0,2.7,0] def r0 ((p0) then scale([1,0,1]) then [[seg_rot]]^-2) + [up0] def r0c (r0) then scale([J]) def r0t (r0) + [gap] def r0b ((p0) then [[seg_rot]]^-2) + [gap] line[arrows=<->](r0c)(r0) line(r0b)(r0t) def r0m ((r0) - (O) + (r0c) - (O)) / 2 + (O) special |\node[ann] at #1 {$r_0$};|(r0m) } def labeled_cone { % the "ghost" of the entire cone sweep[draw=lightgray,cull=false] { N-1, [[seg_rot]] } line(p0)(p1) % for the highlighted face, we need explicit points def p00 (p0) then [[seg_rot]]^-1 def p10 (p1) then [[seg_rot]]^-1 def p01 (p0) def p11 (p1) %polygon[showpoints=true](p00)(p10)(p11)(p01) polygon[fillcolor=red,dotsep=semitransparent](p00)(p10)(p11)(p01) % TikZ does not have a showpoints option. Use dots. dots(p00)(p10)(p11)(p01) % TikZ special for labels. special|\fill[black,font=\footnotesize] #1 node [above] {$P_{00}$} #2 node [below] {$P_{10}$} #3 node [above] {$P_{01}$} #4 node [below] {$P_{11}$};| (p00)(p10)(p01)(p11) def mid ((p00)-(O)+(p10)-(O)+(p11)-(O)+(p01)-(O))/4+(O) % The TikZ arc operation starts at the current point. We therefore % need to shift it to get mid to be the center of the arc special|\draw #1+(-60:.25) [yscale=1.3,->] arc(-60:240:.25);| [lay=over](mid) def mid_left ((p00)-(O)+(p10)-(O))/2+(O) def mid_right ((p11)-(O)+(p01)-(O))/2+(O) special|\path[font=\footnotesize] #1 node[left] {$j$} #2 node[right] {$j\hbox{$+$}1$};| (mid_left)(mid_right) def top_lbl (p01) then [[seg_rot]]^2 def bot_lbl (p11) then [[seg_rot]]^2 special|\path[font=\footnotesize] #1 node[right] {$i\hbox{$=$}0$} #2 node[right] {$i\hbox{$=$}1$};| (top_lbl)(bot_lbl) } def cone {labeled_cone} <> {basic_cone} put { view((10,4,2)) } {cone} % Cool trick: lay = under forces this to be output first in the tikz picture block. special |\tikzstyle{ann} = [fill=white,font=\footnotesize,inner sep=1pt]|[lay=under] global { language tikz }