Trying to make a loop for a definite integral
조회 수: 2 (최근 30일)
이전 댓글 표시
I am trying to create a loop for the following definite integral for y=0 to y=20 with step 0.01 in order to make a plot.any ideas are welcome
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/149104/image.png)
댓글 수: 0
답변 (1개)
Greig
2015년 9월 16일
Is something like this what you are after?
clear vars *;
% Define the function - Note I simplify the denominator in the last term
% Double check the numbers are typed correctly!
F =@(x) 2.906663106 .* x .* ( 1 ./(1 + 1.38 .* x.^4)).^0.4311594203...
- ( (3.458929096 .* x.^5) ./ (1 + 1.38 .* x.^4).^1.4311594203 );
% define teh Y values and the steps
Y = 0:0.01:20;
nSteps = length(Y);
% preallocate the integral values for speed
Ival = NaN(nSteps-1, 1);
% Do the loop
for ii = 1:nSteps-1
Ival(ii) = integral(F, Y(ii), Y(ii+1));
end
% Check it all adds up...
% The whole integral
integral(F, 0, 20)
% The sum of the parts
sum(Ival)
% Then add some plotting function here
댓글 수: 0
참고 항목
카테고리
Help Center 및 File Exchange에서 Loops and Conditional Statements에 대해 자세히 알아보기
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!