Discontinuities when computing integration of error functions using integral function
조회 수: 5 (최근 30일)
이전 댓글 표시
I am trying to integrate a function over a region in different time intervals. The integration looks something like this.
fun_uz = @(u)1./sqrt(u).*exp(-Z.^2./(2.*u));
fun_Y = @(u)(erf((Y+B)./sqrt(2.*u))-erf((Y-B)./sqrt(2.*u)));
fun_Z = @(u)(erf((X+L+u)./sqrt(2.*u))-erf((X-L+u)./sqrt(2.*u)));
fun = @(u)inc.*fun_uz(u).*fun_Y(u).*fun_Z(u);
fint = integral(fun,0,upperl);
The variable 'upperl' is the upper limit of the integral function. I have to perform this integration over different X,Y, and Z regions and different 'upperl' values. I am getting profiles which are discontinuous for different 'upperl' values. I have shown here profiles at few different 'upperl' values.
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/278091/image.jpeg)
I am not able to understand why the discontinuity are occuring, any help is greatly appreciated. Thanks.
댓글 수: 6
darova
2020년 3월 19일
time = 1e-2; % The variable which is changed to generate different contour plots
xs = (-500:10:10)./1000;
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/278262/image.png)
채택된 답변
darova
2020년 3월 19일
The function you are trying to integrate looks like following
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/278274/image.png)
Put these lines inside for loops
ezplot(fun,[0 upperl])
pause(0.01)
When time > 0.1 upperl is big. When you call integral you don't know how many points it takes. Maybe it misses something
댓글 수: 5
추가 답변 (1개)
Walter Roberson
2020년 3월 20일
Change the integral to
fint = integral(fun,0,upperl, 'waypoints', L-X);
You have two erf that only have an input near 0 (and so a measurable output) near-ish -(X+L) to -(X-L) . Some of your integral() calls just happened to evaluate near there, and some of them did not happen to evaluate near there and predicted that there was nothing interesting in that area. The above forces evaluation near that area.
참고 항목
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!