xxxxxxxxxx
def _(N=[5,4,..,1],pp=[5,10,..,30]):
f1(x,y,z,t)=x^2+y^2-.5*z*t^2
f2(x,y,z,t)=x^2+y^2+.5*z*t^2
p1=sum([implicit_plot3d(
f1(x,y,z,i),(x,-i,i),(y,-i,i),(z,-2,2),
color=hue(sin(.24*i)),opacity=1-.15*i,
plot_points=pp,frame=False) for i in [1..N]])
p2=sum([implicit_plot3d(
f2(x,y,z,i),(x,-i,i),(y,-i,i),(z,-1,0),
color=hue(sin(.24*i)),opacity=.3,
plot_points=pp,frame=False) for i in [1..N]])
p3=plot3d(lambda x,y:0,(-N,N),(-N,N),
color='steelblue',opacity=.5)
(p1+p2+p3).show()
No comments:
Post a Comment