Trapezodial Method using while loop

Question

0 个投票

Hi for my numerics class i have an assignment, which is supposed to be quite easy but I can't seem to find my mistake. We have to solve the integral of sin(x) within 0 and pi using the trapezodial method and double the interval until the solution is close enough to the true value 2. So far my code looks like this, and i'm not supposed to use trapz.

clear
f=@(x) sin(x);
a=0;
b=pi;
n=1;
h=(b-a)/n;
tol=10^-6;
I=h*s
s=0.5*(f(a)+f(b))
i=0
while abs(2-I)>=tol
   n=n*2;
   a=b-(b/n);
   b=b/n;
   h(n)=((b*(n-1)/n)-(a/n))/n;
   s(n)=0.5*(f(a)+f(b));
   I=sum(h(n)*s(n));
   i=i+1; 
end
I

but it takes a really long time and gives me the following error:

Requested 1073741824x1 (8.0GB) array exceeds maximum

array size preference. Creation of arrays greater

than this limit may take a long time and cause MATLAB

to become unresponsive. See array size limit or

preference panel for more information.

Error in NumMeth3 (line 18)

h(n)=((b*(n-1)/n)-(a/n))/n;

I don't know why this gets so 'big'. I've tried a lot but nothing seems to work, so I'd love some help!

Thanks in advance.

5 个评论
显示 3更早的评论隐藏 3更早的评论

LauraB 2020-5-18

在 MATLAB Online 中打开

clear
f=@(x) sin(x);  %Function for integration
a=0;            %Lower limit 
b=pi;           %Upper limit 
n=1;            %Number of intervals
h=(b-a)/n;      %Size of each interval 
tol=10^-6;      %Tolerance for solution
s=0.5*(f(a)+f(b))%Second part of the equation for Trapezodial method 
I=h*s           %Putting both parts together by multiplying
i=0             %How often is loop repeated until criteria is met 
while abs(2-I)>=tol %repeat until the solution is within the tolerance 
   n=n*2;           %double the intervals for every time loop is repeated
   a=b-(b*(n-1)/n); %new lower limit 
   b=b-b/n;           %new upper limit
   h=(b-a)/n;       %size of the new interval  
   s=0.5*(f(a)+f(b));%Second part of equation for new limits
   I=sum(h*s);      %Summing the solutions for each interval
   i=i+1;  %Counting iterations
end
I

Maybe it makes more sense what I'm trying to do now? Also I'm aware that the upper limit b=pi is not included anymore, but I can't think of a way to change it properly.

I took away the (n) and now it's not giving me an error anymore but also no solution. If I run it it now gives me NaN as soution for I, s, h and a

Sorry but I'm really new to all this ..

LauraB 2020-5-18

Okay, i guess this needs a whole different approach then..

I wrote down the first few iterations on paper and then tried to 'translate' it into code, but it seems I haven't worked thoroughly enough. I'm gonna try it your way this time :)

Rik 2020-5-18

在 MATLAB Online 中打开

One hint that should help dividing the range 0 to pi in segments:

linspace(0,pi,n)

请先登录，再进行评论。

请先登录，再回答此问题。

Follow Question

Answer 1

David Hill 2020-5-18

在 MATLAB Online 中打开

0 个投票

f=@(x) sin(x);
tol=1e-6;
x=linspace(0,pi,2);
A=sum(movsum(f(x),2,'Endpoints','discard')/2.*diff(x));
n=2;
while abs(2-A)>=tol
    x=linspace(0,pi,n+1);
    A=sum(movsum(f(x),2,'Endpoints','discard')/2.*diff(x));
    n=2*n;
end

2 个评论
显示无隐藏无

Rik 2020-5-18

Since this is homework I wouldn't encourage posting complete working examples. Also, without comments it isn't immediately obvious that this is an implementation of the trapezoidal method.

LauraB 2020-5-18

I agree, but it sometimes helps to take in a different view to find your own mistakes. So thank you anyway!

请先登录，再进行评论。

Trapezodial Method using while loop

5 个评论
显示 3更早的评论隐藏 3更早的评论

回答（1 个）

2 个评论
显示无隐藏无

类别

标签

Community Treasure Hunt

Trapezodial Method using while loop

5 个评论 显示 3更早的评论 隐藏 3更早的评论

回答（1 个）

2 个评论 显示 无 隐藏 无

类别

标签

另请参阅

Community Treasure Hunt

5 个评论
显示 3更早的评论隐藏 3更早的评论

2 个评论
显示无隐藏无