- 积分
- 0
- 注册时间
- 2007-5-13
- 仿真币
-
- 最后登录
- 1970-1-1
|
发表于 2009-12-7 17:24:09
|
显示全部楼层
来自 上海
尝试了一下,设a=1,代码:
> restart;
> equ1 := cos(u)*(cos((1/2)*u)*sin(v)-sin((1/2)*u)*sin(2*v)+3);
> equ2 := sin(u)*(cos((1/2)*u)*sin(v)-sin((1/2)*u)*sin(2*v)+3);
> equ3 := sin((1/2)*u)*sin(v)+cos((1/2)*u)*sin(2*v);
> with(plots);
> equ00 := plot3d([equ1, equ2, equ3], u = -10 .. 10, v = -10 .. 10, numpoints = 5000);
> equ11 := plot3d([equ1, equ2, 7], u = -10 .. 10, v = -10 .. 10, numpoints = 5000);
> equ22 := plot3d([equ1, 7, equ3], u = -10 .. 10, v = -10 .. 10, numpoints = 5000);
> equ33 := plot3d([7, equ2, equ3], u = -10 .. 10, v = -10 .. 10, numpoints = 5000);
> display(equ00, equ11, equ22, equ33);
|
本帖子中包含更多资源
您需要 登录 才可以下载或查看,没有账号?注册
×
评分
-
1
查看全部评分
-
|