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 Comments

"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.
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?
Could you give us some sample input values?
Are you using an old MATLAB that does not have integrate() ?
@Dimitrios: Did you see my comment?
i find the solution with the command integral(fun,0,1,'ArrayValue',true)

Sign in to comment.

Answers (0)

Categories

Find more on MATLAB in Help Center and File Exchange

Tags

Community Treasure Hunt

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

Start Hunting!