>i:=0;
>for h from 0.0 by 0.01 while h <= 1.14 do
i:=i+1:
upx[i]:=fsolve(-sqrt(1-x^2)=sqrt(1-x^2)-h,x=-1.1..1.1,complex):
s[i]:=upx[i]*h+2*int(sqrt(1-x^2),x=upx[i]..1);
end do;
>a:=seq([j*3.14/316,s[j]],j=1..115);
>plot({[a]},style=point);
>i:=115;
>for h from 1.15 by 0.01 while h <= 2 do
i:=i+1;
upx[i]:=fsolve(sqrt(1-x^2)=-sqrt(1-x^2)
+3
.1415926-h,x=-1.1..1.1,complex);
dnx[i]:=fsolve(sqrt(1-x^2)-h=-sqrt(1-x^2));
s[i]:=int(sqrt(1-x^2)-(-sqrt(1-x^2)+3.1415926-h),x=-upx[i]..0)+dnx[i]*h+2*int(sqrt(1-x^2),x=dnx[i]..1)
end do;
>a:=seq([j*3.14/316,s[j]],j=115..201);
>plot({[a]},style=point); >i:=201;
为什么这个程序在maple里运行不出来？？？？
> for h from 2 by 0.01 while h <= 3.14 do
upx[i]:=fsolve(sqrt(1-x^2)=-sqrt(1-x^2)+3.1415926-h,x=-1.1..1.1,complex);
s[i]:=int(sqrt(1-x^2)-(-sqrt(1-x^2)+3.1415926-h),x=-upx[i]..0)+3.1415926/2;
end do;
>a:=seq([j*3.14/316,s[j]],j=201..316);
>plot({[a]},style=point);