登录  
 加关注
   显示下一条  |  关闭
温馨提示!由于新浪微博认证机制调整,您的新浪微博帐号绑定已过期,请重新绑定!立即重新绑定新浪微博》  |  关闭

徐小湛的博客

高等数学课件、函数图形动画、数学软件应用

 
 
 

日志

 
 

两张相交的球面  

2013-04-26 18:14:44|  分类: 高等数学图形动画 |  标签: |举报 |字号 订阅

  下载LOFTER 我的照片书  |

两张相交的球面 - calculus - 徐小湛的博客两张相交的球面 - calculus - 徐小湛的博客

with(plots):R:=3:r:=2:
qumian:=plot3d([r*sin(phi)*cos(theta),r*sin(phi)*sin(theta),R+r*cos(phi)],phi=Pi/2+arcsin(r/R/2)..Pi,theta=0..2*Pi,color=green):
qumian1:=plot3d([r*sin(phi)*cos(theta),r*sin(phi)*sin(theta),R+r*cos(phi)],phi=0..Pi/2+arcsin(r/R/2),theta=0..2*Pi,color=yellow):
qumian2:=plot3d([R*sin(phi)*cos(theta),R*sin(phi)*sin(theta),R*cos(phi)],phi=0..Pi,theta=0..2*Pi,color=blue,style=wireframe):
x_axis:=plot3d([u,0,0],u=-R-1..R+1,v=0..0.01,thickness=3):
y_axis:=plot3d([0,u,0],u=-R-1..R+1,v=0..0.01,thickness=3):
z_axis:=plot3d([0,0,u],u=-R-1..R+r+1,v=0..0.01,thickness=3):
display(qumian,qumian1,qumian2,x_axis,y_axis,z_axis,scaling=constrained); 两张相交的球面 - calculus - 徐小湛的博客
 
两张相交的球面 - calculus - 徐小湛的博客
 
两张相交的球面 - calculus - 徐小湛的博客  两张相交的球面 - calculus - 徐小湛的博客
with(plots):R:=3:r:=2:
qumian:=plot3d([r*sin(phi)*cos(theta),r*sin(phi)*sin(theta),R+r*cos(phi)],phi=Pi/2+arcsin(r/R/2)..Pi,theta=0..2*Pi,color=red):
qumian1:=plot3d([r*sin(phi)*cos(theta),r*sin(phi)*sin(theta),R+r*cos(phi)],phi=0..Pi/2+arcsin(r/R/2),theta=0..2*Pi,color=green,style=wireframe):
qumian2:=plot3d([R*sin(phi)*cos(theta),R*sin(phi)*sin(theta),R*cos(phi)],phi=0..Pi,theta=0..2*Pi,color=blue,style=wireframe):
x_axis:=plot3d([u,0,0],u=-R-1..R+1,v=0..0.01,thickness=3):
y_axis:=plot3d([0,u,0],u=-R-1..R+1,v=0..0.01,thickness=3):
z_axis:=plot3d([0,0,u],u=-R-1..R+r+1,v=0..0.01,thickness=3):
display(qumian,qumian1,qumian2,x_axis,y_axis,z_axis,scaling=constrained);
 
     
两张相交的球面 - calculus - 徐小湛的博客

面积最大的情形 r:=R*4/3:
 
两张相交的球面 - calculus - 徐小湛的博客两张相交的球面 - calculus - 徐小湛的博客
with(plots):R:=3:r:=R*4/3:
qumian:=plot3d([r*sin(phi)*cos(theta),r*sin(phi)*sin(theta),R+r*cos(phi)],phi=Pi/2+arcsin(r/R/2)..Pi,theta=0..2*Pi,color=red):
qumian1:=plot3d([r*sin(phi)*cos(theta),r*sin(phi)*sin(theta),R+r*cos(phi)],phi=0..Pi/2+arcsin(r/R/2),theta=0..2*Pi,color=green,style=wireframe):
qumian2:=plot3d([R*sin(phi)*cos(theta),R*sin(phi)*sin(theta),R*cos(phi)],phi=0..Pi,theta=0..2*Pi,color=blue,style=wireframe):
x_axis:=plot3d([u,0,0],u=-R-1..R+1,v=0..0.01,thickness=3):
y_axis:=plot3d([0,u,0],u=-R-1..R+1,v=0..0.01,thickness=3):
z_axis:=plot3d([0,0,u],u=-R-1..R+r+.31,v=0..0.01,thickness=3):
display(qumian,qumian1,qumian2,x_axis,y_axis,z_axis,scaling=constrained);
  评论这张
 
阅读(2221)| 评论(2)

历史上的今天

评论

<#--最新日志,群博日志--> <#--推荐日志--> <#--引用记录--> <#--博主推荐--> <#--随机阅读--> <#--首页推荐--> <#--历史上的今天--> <#--被推荐日志--> <#--上一篇,下一篇--> <#-- 热度 --> <#-- 网易新闻广告 --> <#--右边模块结构--> <#--评论模块结构--> <#--引用模块结构--> <#--博主发起的投票-->
 
 
 
 
 
 
 
 
 
 
 
 
 
 

页脚

网易公司版权所有 ©1997-2018