源程序:
syms t xi yj
fuchang=int ( (t-sin(t))*xi+(1-cos(t))*yj ,t,0,2*pi),
运行得:
fuchang =
2*xi*pi^2+2*yj*pi
即结果为
2*pi^2 *i+2*pi *j.
我算了一下,你的程序是没错的呀,不过最后一句应该改成I=int(F*ds,t,0,pi),这样才是正向上半椭圆吧。结果是I
=-2/15*c*(-2*c^4+15*b^4)/a/b^4。
而且,如果把t的范围取成pi到2*pi,可以得到结果为I=2/15*c*(-2*c^4+15*b^4)/a/b^4。你把t的范围取成0到2*pi,那么肯定结果为零。