xxxxxxxxxxr.eval("""X3<-Z3<-X4<-Z4<-Y5<-Z5<-seq(from=-5,to=5,by=.1)cylinder1<-function(p,q,X,Z){ for (x in X){for (z in Z){ y<-q*(1-x^2/p^2)^0.5; v<-c(x,y,z,x,-y,z) if (sum(is.na(v))==0){ write.table(matrix(v,nrow=2,ncol=3,byrow=TRUE), file ='cylinder1.csv', append=TRUE,quote=FALSE, col.names=FALSE,row.names=FALSE)}}}}cylinder2<-function(p,q,X,Z){ for (x in X){for (z in Z){ y<-q*(x^2/p^2-1)^0.5; v<-c(x,y,z,x,-y,z) if (sum(is.na(v))==0){ write.table(matrix(v,nrow=2,ncol=3,byrow=TRUE), file='cylinder2.csv', append=TRUE,quote=FALSE, col.names=FALSE,row.names=FALSE)}}}}cylinder3<-function(p,Y,Z){ for (y in Y){for (z in Z){ x<-y^2/2/p; v<-c(x,y,z) write.table(matrix(v,nrow=1,ncol=3,byrow=TRUE), file ='cylinder3.csv', append=TRUE,quote=FALSE, col.names=FALSE,row.names=FALSE)}}}cylinder1(3,4,X3,Z3); cylinder2(2,3,X4,Z4); cylinder3(2,Y5,Z5)""")xxxxxxxxxximport pylab,numpyfrom mpl_toolkits.mplot3d import Axes3Dfrom matplotlib.lines import Line2Dfiles=['cylinder1.csv','cylinder2.csv', 'cylinder3.csv']colors=['#36ff36','#36ffff','#ff36ff']labels=['$x^2/3^2+y^2/4^2=1$', '$x^2/2^2-y^2/3^2=1$', '$y^2=2*2*x$']XYZ=[numpy.array(numpy.loadtxt(files[i],\ delimiter=' ',unpack=True)) for i in range(3)]fig=pylab.figure(figsize=(6.5,6.5))ax=fig.add_subplot(111,projection='3d')[ax.scatter(XYZ[i][0],XYZ[i][1],XYZ[i][2],s=1,\ c=colors[i],marker='1') for i in range(3)]fake2Dlines=[Line2D([0],[0],linestyle='none',\ c=c,marker='1') for c in colors]ax.legend(fake2Dlines,labels,loc=9) pylab.tight_layout(); pylab.show()
No comments:
Post a Comment