>
|
use plots, plottools in
module()
local p1,p2,p3,p4,F,f;
F:=1+2*r*cos(theta)+3*r*sin(theta);
p1:=plot([cos(2*t),t,t=-Pi..Pi],coords=polar,color=red,thickness=2):
p2:=coordplot(polar,color=[gray,gray]):
f:=transform((x,y)->[x,y,0]):
p3:=plot3d(F,r=0..cos(2*theta),theta=-Pi/4..Pi/4,coords=cylindrical,filled=true,lightmodel=none,glossiness=0):
p4:=display(f(p1),f(p2),p3,lightmodel=none,axes=frame,labels=[x,y,z],tickmarks=[3,3,3],orientation=[-100,70,0]);
print(p4);
end module:
end use:
|