Items tagged with integral integral Tagged Items Feed

Why this simplifes:

z1:=n*(Int(cos(x)^(n-2), x))-(Int(cos(x)^(n-2), x));


But when adding an extra term to z1, it no longer simplfies the above any more:

z2 := cos(x)^(n-1)*sin(x)+n*(Int(cos(x)^(n-2), x))-(Int(cos(x)^(n-2), x));


You can see the second term, which is z1, was not simplfied any more.

Why? And how would one go about simplifying z2 such that the second term gets simplfies as with z1, but while using z2 expression. It seems simplify stopped at first term and did not look ahead any more?

Maple 18.02, windows.

I am trying to find the root of an equation that involves a procedure and a definite integral (solved numerically). Of course, I don't need the root to be found symbolically, but numerically would be fine. The problem is, I keep getting the error

"Error, (in fsolve) Can't handle expressions with typed procedures"

whenever I try to solve it. Anyone have any ideas? My worksheet is here: 


Excuse me,


  I have the following input




I would like to obtain

1/(k+1)* (x2^(k+1)), k<>-1


but maple gives me just the same integral. Is there any way to get the expression?


Thank you!




hi everyone , i need your point of view in my question,any help would be appreciated in advance .

we have a discrete function named g(t) and a continous function f(t) in in convolution integral just like this :
int(f(t-x)*g(x),x=0..t) ; 
we have just g(x) in some special points int the interval (0..t) , thus i need to convert this integral to a series.
how should i do this ? can anyone help or any idea ? i need at first a mathematical solution or idea about how to do this and then, how to do this in a software ?
tnx again.

I'd like to pay attention to the recent article "The Misfortunes of a Trio of Mathematicians Using Computer Algebra Systems. Can We Trust in Them?"

In particular, the authors consider the integral


stating "Both Mathematica and Maple return zero as the answer to this calculation. Yet this cannot be correct, because the integrand is clearly positive and nonzero in the indicated region". Unfortunately, they give only the Mathematica command to this end.

Of course, the integral under consideration is complicated so the the simple-minded trials

int(evalc(abs(exp((2*Pi*I)*x)+exp((2*Pi*I)*y))), [x = 0 .. 1, y = 0 .. 1]);


VectorCalculus:-int(evalc(abs(exp((2*Pi*I)*x)+exp((2*Pi*I)*y))), [x,y]=Rectangle( 0 .. 1, 0 .. 1));

fail. However,this can be found with Maple (I think with Mathematica too.) in such a way.


A := evalc(abs(exp((2*Pi*I)*x)+exp((2*Pi*I)*y)))




B := simplify(A, trig)






combine(op(B)[1], x)



C := eval(B, op(B)[1] = combine(op(B)[1], x))



int(C, [x = 0 .. 1, y = 0 .. 1])









I wanted to get the following function integrated with Maple 17:

I used the following command to get a numerical solution for my choise of limits:

Unfortunately, I just get the Integral itself back. Also, taking pi/2 as the upper limit for x does improve anything. I also tried the AllSolutions option but it does not help as well as first trying to integrate only over y and leaving x as it is in order to get at least one integration done.

I hope someone can help. I do not necessarily need to solve it for myself, so if someone gets the result and posts it here, I'm okay with that. Anyway, being able to calculate it myself would be even better.

1) Consider the followin Integral:


How can I tell maple to give the results in terms of hyperbolic and harmonic funtions?
The maple returns the answer in term of exponentials with imaginary powers.
Also, I could not use the "simplify" command to turn the answer returned by maple into an answer in terms of what I want! Please see the code below.

2) Is there a way that I can update the maple integration database manaully?

Below is the function that I have.


f := (t-1)^(1/3)


b[n] := 2/p*(Int(f*sin(2*Pi*n*t/p), t = 0 .. p))


I also included a picture below to show what it is doing. Some help would be greatly appreciated. All I need to know is why maple doesn't want to evaluate bn?


Maple Code


Hi everyone!

I wander whether maple can solve the integral of trigonometric series with parametal N, the number of sereis, and how. The formation is showed as below. N is a  variable and 'm' belongs to 'k', 'n' belongs to 'l'.

the intergral of series and the orthogonality conditions

A := int(int(sum(sum(cos(2*k*Pi*x/a)*(1-cos(2*l*Pi*y/b))*(1-cos(2*m*Pi*x/a))*(1-cos(2*n*Pi*y/b)), k = 1 ..N), l = 1 .. N), x = 0 .. a), y = 0 .. b)

orthogonality codition 1:

OrthCondition1 := int(sum(cos(2*k*Pi*x/a)*cos(2*m*Pi*x/a), k = 1 .. N), x = 0 .. a) = (1/2)*a

orthogonality codition 2:

OrthCondition2 := int(sum(cos(2*l*Pi*x/b)*cos(2*m*Pi*x/b), l = 1 .. N), x = 0 .. b) = (1/2)*b;



The alpha and beta are 'randomly' chosen, both >0, to produce some values, hence a value for the target function to optimize.

I am having trouble to evaluate the expression in the middle of the calculation.

Sometimes, the 'chosen' alpha and beta works fine with 'method = _d01amc' for numerical integration. and such way is really "fast".


But sometimes, the the 'chosen' alpha and beta will fail. Even when both of them are perfectly defined. and the integral can be easily evaluated using Int() and then evalf().


So what's be best way to proceed?




I have a line integral and i want to calculate and graph this line integral.

The line integral comes from a physical problem, the integrand is good in the interval [0, 2*Pi] and the value of the integrand at zero = the value of the integrand at 2*Pi.

The problem is when i use the trapizidal rule to calculate the line integral.

The results is so bad, i.e.  the value of the line integral at zero not equal to the value of the line integral at 2* Pi.

where is the problem.

Any one could help me.

Attached is my sheet.

On my system if just execute the whole worksheet maple return the last statement's integral un evaluated, however at this point if re-execute the statements from the definition of L2, it looks that maple has evaluated the last statement's integral.


What am I missing here, any help?


I am working on fopurier analysis and below is what I have. The method I used worked fine when I had a piecewise continuous function, but now the bottom int function below is not calculating the integral.


f := 2*sin(Pi-Pi*exp(-t))


A[0] := (1/3)*(Int(f, t = 0 .. 3))

A[0] := (1/3)*(int(f, t = 0 .. 3))   I think this one is evaluating it correctly. The bold and italicized below is the answer from this calulation

A[0] := (2/3)*Si(Pi)-(2/3)*Si(exp(-3)*Pi)


a[n] := (2/3)*(Int(f*cos(2*Pi*n*t*(1/3)), t = 0 .. 3))  this prints the integral like I want it to.

a[n] := (2/3)*(int(f*cos(2*Pi*n*t*(1/3)), t = 0 .. 3))  this just reporints the integral, I want it to evaluate it. Am I doing something wrong? Even with the assume(n, integer); it still just prints the integral.



Is there any help that y'all can provide? My professor is not ver proficent in the Math softwares out there.


 While calculating an integral including complex numbers, I have encountered with the output "undefined if a+ib>0". What does this mean?A complex number bigger than zero???


I read this thread, but it didnt really help. (i noticed an error in the expression pj, now fixed!)

See this file above.

The function needs to be maximized is "loglik", with constraints,

f0>0, sigma>0 and mu is FREE. Only three variables.


Using some calculus and by changing variables, I was able to use Gauss–Hermite quadrature, and do some approximation, to get the values:





I am wondering if Maple can do this straight away. By that, I mean, without doing much further,  how could I get the values? (most efficiently and quickly)





1 2 3 4 5 6 7 Last Page 1 of 32