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 Commenti

Jan
Jan il 27 Feb 2019
Modificato: Jan il 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.
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)

Accedi per commentare.

Risposte (0)

Categorie

Scopri di più su Mathematics in Centro assistenza e File Exchange

Tag

Community Treasure Hunt

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

Start Hunting!

Translated by