Question

Task 2: (1 Mark) - Scripts Write a MATLAB script which plots the absolute error between your mySqrt() function and the built-in sqrt() function for a variety of N values with: N=1:20 • n=323158 err-le-3. Your plot should clearly show that error decreases until "saturating" somewhere below err. Because the algorithm "converges” (increased in accuracy) very quickly the true accuracy will be several orders of magnitude better than the err value. You should use semilogy () instead of plot() to plot with a logarithmic y-axis. To do this replace plot () with semilogy(). Both functions take the same arguments. You may ask your demonstrator for help regaring the use and interpretation of semilogy(). Plot title and axis labels are not required.

Answer 1

Please find the answer below.
Please do comments in case of any issue. Also, don't forget to rate the question. Thank You.

mysqrt.m

function x = mysqrt( n,N,err )
xn = n./2;
for i=1:N
x=0.5 *(xn+(n./xn));
if abs(x-xn)<err
return
end
xn = x;
end

end

script.m

clc
clear
n = 323158;
err = 1e-3;
data = zeros(1,20);
for N=1:20
myRoot = mysqrt(n,N,err);
inbuiltRoot = sqrt(N);
data(N) = abs(myRoot-inbuiltRoot);
end

x = 1:20;
semilogy(x,data);

output

Figure 1 Х G U File Edit View Insert Tools Desktop Window Help akpo 105 АКРО ticci and Wind =_wo m 104 = (n. qrt yRC 103 102 0 4 6 8 10 12 14 16 18 20 N