Numerical Integration of scalar function

Hi Dear researchers. I want to calculate numerically the integral f(z)=exp(i*I*pi*z+cos(z)) where n eisai a matrix. i am using the above command
function out=I1mn(g,N,L,d,ai,si)
M=length(ai);
I=-N:N;
I=I(ones(2*N+1,1),:);
I=I.'-I;
fp=(@(z,n)exp(j*2*pi*n*z+(g)*(d*cos*2*pi*z/L))));
out=quadv(@(z)fp(z,I),ai,ai+si)
but when i execute my code i take the above message
Warning:Maximum function count exceeded; singularity likely.
>in quadv
in I1mn
How i can solve this problem to take good results
With best Regards DImitrios Gerontitis.

6 Kommentare

Jan
Jan am 27 Feb. 2019
Bearbeitet: Jan am 27 Feb. 2019
"n eisai a matrix" ? f(z)=exp(i*I*pi*z+cos(z)) does not contain an "n"?
What are the variables g,N,L,d,ai,si?
I do not see, why you implement:
exp(i*I*pi*z + cos(z))
as
exp(j*2*pi*n*z + (g)*(d*cos*2*pi*z/L)))
This is at least confusing.
Bob Thompson
Bob Thompson am 27 Feb. 2019
Does the code run to completion and give you a result? If so, how do you know that the result isn't 'good'? What qualfies as a 'good' result with respect to the answers you are getting?
Unfortunately i don't take any result. Do you know any numerical integration routine?
Walter Roberson
Walter Roberson am 27 Feb. 2019
Could you give us some sample input values?
Are you using an old MATLAB that does not have integrate() ?
Jan
Jan am 27 Feb. 2019
@Dimitrios: Did you see my comment?
i find the solution with the command integral(fun,0,1,'ArrayValue',true)

Melden Sie sich an, um zu kommentieren.

Antworten (0)

Kategorien

Mehr zu Mathematics finden Sie in Hilfe-Center und File Exchange

Tags

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by