Question

numerical method lab

A5.2. Write a matlab program to find out the root of equation f (x)-x*-3x - 1, using false-position method. Use initial and upper guesses -1 and 1.5 False Position method: xr-Xu Explanation of False Position f(x)(x-xu) Hints change) xu far-feval (foxr) ; root-0.3294 Xr

Answer 1

fun.m:

clc;

f=@(x) x^4-3*x-1;

a=2,b=3;

for i=1:10

x0=a;

x1=b;

fprintf('\n Hence root lies between (%.4f, %.0f)',a,b);

x2(i)=x0(x1-x0)/(f(x1)-f(x0))*f(x0);

if(x2(i)>0

b=x2(i);

else

a=x2(i);

end

fprintf('\n Therefore, x2=%.4f \n Here , f(x2)=%.4f',x2(i),f(x2(i)))

p=x2(i);

end

for i = 1: 10

error(i) = p - x2(i);

end

Answer=p;

plot(error)

grid on;

title('plot of error');

xlabel('terations');

ylabel('Error');

Output:

dj How to install Django on Windo x ng Diango models . Diango Girls Tu X © A52. Write A Matlab program T Execute MATLAB/Octa e Online + T Not secure www.compileonline.com/execute_matlab online.php codingground Execute MATLAB/Octave Online (GNU Octave, v4.2.1) A Project ▼ Edit ▼ * Setting ▼ Execute demo.m STDIN .l Result clc; 2 3 1 soctave-qf no-window-system demo.m f-6(x) x^4-3*x-1; a-2; b-3; Hence root lies between(2.0800,3) Therefore, x2-1.8548 Here, f(x2)-5.2720 Hence root 1ies between (2.0800,2) Therefore, x2-1.6496 Here, f(x2)-1.4554 Hence root lies between (2.0800,2) Therefore, x2-1.5828 Here, f(x2)-0.5170 Hence root lies between (2.080,2) Therefore, x2-1.5565 Here, f(x2)-0.1996 Hence root lies between(2.0800,2) Therefore, x2-1.5464 Here, f(x2)-0.8795 Hence root 1ies between (2.0800,2) Therefore, x2-1.5424 Here, f(x2)-0.8321 Hence root lies between (2.0800,2) Therefore, x2-1.5487 Here, f(x2)-0.8130 Hence root lies between (2.080,2) Therefore, x2-1.5401 Here, f(x2)-0.8053 Hence root lies between(2.0800,2) Therefore, x2-1.5398 Here, f(x2)-0.8021 Hence root 1ies between (2.0800,2) Therefore, x2-1.5397 Here, f(x2)-0.8009Answer 1.5397 5 for i-1:18 xe-a; x1-b; fprintf('\n Hence root lies between (%.4f,S.0f)',a,b) x2 (i)-x8-(x1-x8) ((x1)-f(xe))*(xe): if f(x2(i))>e 10 b-x2(1) else a-x2(1); 12 13 fprintf('In Therefore, x2-%.4f In Here, f(x20-%.4f ,x2(i),f(x2(i))) 14 15 end 16 for i-1:1e 17 18 end 19 Answer-p 20 plot (error) 21 grid on; 22 title("Plot of error) 23 xlabel('iterations: 24 ylabel( Error) error(i)-p-x2(i); warning: function /demo.m shadows a core library function O Type here to search :30 PM 3/20/2019 ^恤瘺 x