(*求积分面积*)
Clear["Global`*"];(*Clear all variables*)
(*求解交点*)
sol={x,y}/.Solve[{(x-1)^2+(y-1)^2==1,x*y==1},{x,y}];
(*获得其中用来积分的交点*)
solxy=sol[[2]];
x0=solxy[[1]];(*积分的起始x坐标*)
x1=1;
x2=1+1/Sqrt[2];
fun1[x]:=1+Sqrt[1-(x-1)^2];(*圆的上半部分曲线方程*)
fun2[x]:=1/x;(*双曲线方程*)
fun3[x]:=x;(*对称轴方程*)
area1=Integrate[fun1[x]-fun2[x],{x,x0,x1}];
area2=Integrate[fun1[x]-fun3[x],{x,x1,x2}];
area=2*(area1+area2)
http://bbs.emath.ac.cn/thread-5266-1-1.html
