Trying to find the error in my modified false position method code

Question

Lily 2024-4-16

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2107771-trying-to-find-the-error-in-my-modified-false-position-method-code

评论： Lily 2024-4-25

I am getting this error (Index exceeds the number of array elements. Index must not exceed 1) for my code and unsure why. Any help would be appreciated.

clc;
clear all;
close all;
%Mini Project
%Modified False Position Method
syms f(x);
f(x)=(-4.868e-7)*(4*x^3)+(1.496e-4)*(3*x^2)+(-0.0094)*(2*x)+(-0.3324);
i=1;
xl(i)=150;
xu(i)=200;
iu(i)=0;
il(i)=0;
fu(i)=f(xu(i));
fl(i)=f(xl(i));
xr(i)=xu(i)-((fu(i)*(xl(i)-xu(i))/(fl(i)-fu(i))));
fr(i)=f(xr(i));
ea = zeros(1,i);
if f(xr(i))~=0
    while i<100
        i=i+1;
        if fl(i-1)*fr(i-1)<0
            xu(i)=xr(i-1);
            xl(i)=xl(i-1);
            iu(i)=iu(i-1);
            il(i)=il(i-1)+1;
            if il(i)>=2
                fl(i)=fl(i-1)/2;
                fu(i)=fr(i-1);
            end
        else
            xl(i)=xr(i-1);
            xu(i)=xu(i-1);
            il(i)=il(i-1);
            iu(i)=iu(i-1)+1;
            if iu(i)>=2
                fu(i)=fu(i-1)/2;
                fl(i)=fr(i-1);
            end
        end
        xr(i)=xu(i)-((fu(i)*(xl(i)-xu(i))/(fl(i)-fu(i))));
        fr(i)=f(xr(i));
            ea(i)=abs((xr(i)-xr(i-1))/xr(i))*100;
            if ea(i)<0.05
                break;
            end
    end
end
Index exceeds the number of array elements. Index must not exceed 1.

Error in indexing (line 962)
                R_tilde = builtin('subsref',L_tilde,Idx);
figure(1)
plot(ea)
hold on
grid on
xlabel('Iteration')
ylabel('Percent Approximation Error')
y=-200:10:200;
figure(2)
hold on
grid on
box on
for i=1:1:length(y)
    plot(y,f(y),'LineWidth',2)
end
plot(xr,f(xr),'*','MarkerSize',10);
xlabel('Temperature (K)')
ylabel('Derivative of f(T)')
yline(0);
xr=double(xr);
r_mfpm=(-4.868e-7)*(xr^4)+(1.496e-4)*(xr^3)+(-0.0094)*(xr^2)+(-0.3324*xr)+14.55;   
et=abs((166.222229010761-xr)/-166.222229010761)*100;
         
 

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

James Tursa 2024-4-16

在 MATLAB Online 中打开

To debug this, type the following at the command prompt:

dbstop if error

Then run your code. When the error occurs, the program will pause at the line causing the error with all variables intact. Then you can examine the variables on that line to figure out the cause of the error. In your case, some of the variables are scalars but you are using an index of 2, hence the error. Backtrack in your code to figure out why these variables are not the size you expected.

请先登录，再进行评论。

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

Answer 1

Aman 2024-4-22

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2107771-trying-to-find-the-error-in-my-modified-false-position-method-code#answer_1445511

在 MATLAB Online 中打开

Hi Lily,

You are getting the below error from the below line:

xr(i)=xu(i)-((fu(i)*(xl(i)-xu(i))/(fl(i)-fu(i))));

Here "fl" and "fu" are 1x1 syms, and the code is trying to access the second index, due to which it is erroring out. Please reiterate your code once and make appropriate changes, like accessing the proper index or assigning the proper type to the "fl" and "fu.".

I hope this helps!