I've tried reducing the digits to 8 and to compute the integral directly, both solutions work. Thank you!

I've tried reducing the digits to 8 and to compute the integral directly, both solutions work. Thank you!

The original code was a little too complex so I was reluctant to post it, but here's a simplified version where I still have the same problem:

> Tproc := proc (t, k, b) int((1-(1-exp(-k*(t-s)))^b)*(1-exp(-160*s^2)), s = 0 .. .15, numeric = true, method = _d01ajc) end proc:

> Tproc(2, 10, 2.5)

1.207298877*10^(-9)

> fdiff(Tproc, [2], [2, 10, 2.5])

fdiff(Tproc, [2], [2, 10, 2.5])

I think the problem has something to do with the 'method = _d01ajc'. But if I don't specify this method, the problem is even worse: 

> Tproc := proc (t, k, b) int((1-(1-exp(-k*(t-s)))^b)*(1-exp(-160*s^2)), s = 0 .. .15, numeric = true) end proc:

> Tproc(2, 10, 2.5)

1.207298877*10^(-9)

> fdiff(Tproc, [2], [2, 10, 2.5])

(no output, Maple hangs and crashes after a few minutes)