先把积分化简:
问题即转化为解多项式方程
代码:
p = [1/3 1 0 -3];
roots(p)
结果有三个根:
ans =
-2.2130 + 1.1891i
-2.2130 - 1.1891i
1.4260 + 0.0000i
其中前两个为复数根,第三个为实数根。
参考:How can I solve this integral equation?
syms a x y
y=int(x^2+a,x,0,a)
然后运行
solve('1/3*a^3+a^2=3','a');
vpa(a)
这结果出来了,三个。自己根据需要选择合适解。一般是第一个,后面俩个是复数解