请求积分问题

weilinhy

有一个积分关于x的积分 主要是积分函数含有积分下限
H=0.68; s=1; t=30;
syms x;
f=int(power(x,H-0.5)*power(x-s,H-1.5),x,s,t);

运行结果
Warning: Explicit integral could not be found.
> In sym.int at 58

但是个积分在Mathematica里面可以得到具体数值
很奇怪

如果用
double(int(power(x,H-1.5)*power(x-s,H-0.5),x,s,t));
也不行 有waring

请位高手帮忙看看了
谢谢

rocwoods

1. 
   
2. H = 0.68;
   
3. t = 30;
   
4. IntResult = @(s) quadl(@(x) power(x,H-0.5).*power(x-s,H-1.5),s,t);
   
5. %当s=1时
   
6. s = 1;
   
7. IntResult(s)
   

复制代码

手头没有MATLAB，楼主运行看看

weilinhy

谢谢 recwoods
运行有结果 但是
Warning: Minimum step size reached; singularity possible.
> In quadl at 95
  In @(s) quadl(@(x) power(x,H-0.5).*power(x-s,H-1.5),s,t)
请问可以改进么
多谢了

ChaChing

ff=@(x) power(x,H-0.5).*power(x-s,H-1.5)
ezplot(ff)
的确在s=1处singularity

weilinhy

谢谢chachinq  论坛高手好多 拜谢