[go: up one dir, main page]
Scribd
Upload
57 views4 pages

Program As

This document contains MATLAB code for three numerical methods for solving nonlinear equations: the bisection method, Newton's method, and the false position method. The bisection method iteratively finds the root of a function f(x) = 0 by narrowing the range between two values a and b where f(a) and f(b) have opposite signs. Newton's method iteratively finds the root using the function value and derivative. The false position method iteratively finds the root using the function values at two points to estimate where the function crosses zero.

Uploaded by

keila
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
57 views4 pages

Program As

This document contains MATLAB code for three numerical methods for solving nonlinear equations: the bisection method, Newton's method, and the false position method. The bisection method iteratively finds the root of a function f(x) = 0 by narrowing the range between two values a and b where f(a) and f(b) have opposite signs. Newton's method iteratively finds the root using the function value and derivative. The false position method iteratively finds the root using the function values at two points to estimate where the function crosses zero.

Uploaded by

keila
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
You are on page 1/ 4

function bisect(f,a,b,tol,n)

% Bisection method for solving the nonlinear equation


f(x)=0.
a0=a; b0=b;
iter=0;
u=feval(f,a);
v=feval(f,b);
c=(a+b)*0.5;
err=abs(b-a)*0.5;
disp('_______________________________________________')
disp(' iter a b c f(c) |b-a|/2 ')
disp('_______________________________________________')
fprintf('\n')
if (u*v<=0)
while (err>tol)&(iter<=n)
w=feval(f,c);
fprintf('%2.0f %10.4f %10.4f %12.6f %10.6f %10.6f\n',
iter,a, b, c, w, err)
if (w*u<0)
b=c;v=w;
end;
if (w*u>0)
a=c;u=w;
end;
iter=iter+1;
c=(a+b)*0.5;
err=abs(b-a)*0.5;
end;
if (iter>n)
disp(' Method failed to converge')
end;
else
disp(' The method cannot be applied f(a)f(b)>0')
end;
% Plot f(x) in the interval [a,b].
fplot(f, [a0 b0])
xlabel('x'); ylabel('f(x)'); grid
function newton(f,df,x0,tol,n)
% Newton's method for solving the
nonlinear equation f(x)=0.
iter=0;
u=feval(f,x0);
v=feval(df,x0);
err=abs(u/v);
disp('__________________________________
____________')
disp(' iter x f(x) df(x) |xn+1-xn| ')
disp('__________________________________
____________')
fprintf('%2.0f %12.6f %12.6f %12.6f\n',
iter, x0, u, v)
while (err>tol)&(iter<=n)&(v~=0)
x1=x0-u/v;
err=abs(x1-x0);
x0=x1;
u=feval(f,x0);
v=feval(df,x0);
iter=iter+1;
fprintf('%2.0f %12.6f %12.6f %12.6f
%12.6f\n',iter,x0,u,v,err)
end
if (v==0)
disp(' division by zero')
end
if (iter>n)
disp(' Method failed to converge')
end
function falsep(f,a,b,tol,n )
% False position method for solving the
nonlinear equation f(x)=0.
a0=a; b0=b;
iter=0;
u=feval(f,a);
v=feval(f,b);
c=(v*a-u*b)/(v-u);
w=feval(f,c);
disp('__________________________________
____________________________')
disp(' iter a b
c f(c) |b-a| ')
disp('__________________________________
____________________________')
fprintf('\n')
if (u*v<=0)
while (abs(w)>tol)&(abs(b-
a)>tol)&(iter<=n)&((v-u)~=0)
w=feval(f,c);
fprintf('%2.0f %12.4f %12.4f %12.6f
%10.6f %10.6f\n', iter, a,b, c, w,
abs(b-a))
if (w*u<0)
b=c;v=w;
end;
if (w*u>0)
a=c;u=w;
end;
iter=iter+1;
c=(v*a-u*b)/(v-u);
end;
if (iter>n)
disp(' Method failed to converge')
end;
if (v-u==0)
disp(' Division by zero')
end;
else
disp(' The method cannot be applied
f(a)f(b)>0')
end;
fplot(f,[a0 b0])
xlabel('x');ylabel('f(x)'); grid

You might also like