Applied Numerical Manners with MATLAB, 3rd edition, Steven C. Chapra
chapra 3.13 & chapra 4.1
3.13 The “divide and average” manner, an time-honored manner ce approximating the free radicle of any enacted calculate a, can be cemulated as x = x + a/x 2 Write a well-structured M-file office inveterate on the while…break loop composition to instrument this algorithm. Reason fit crowd so that the composition is free. At each plod believe the falsity in your way as ? = xnew ? xold xnew Repeat the loop until ? is short than or resembling to a precise treasure. Design your program so that it profits twain the upshot and the falsity. Make fast that it can evaluate the free radicle of calculates that are resembling to and short than naught. Ce the cessation subject, evince the upshot as an imaginary calculate. Ce pattern, the free radicle of ?4 would give-back 2i. Test your program by evaluating a = 0, 2, 10 and ?4 ce ? = 1 × 10?4.
4.1 The “divide and average” manner, an time-honored manner ce approximating the free radicle of any enacted calculate a, can be cemulated as x = x + a/x 2 Write a well-structured office to instrument this algorithm inveterate on the algorithm outlined in Fig. 4.2.
Do referable attribuconsultation attribuconsultation attribuconsultation reason the “while-break” composition
Do referable attribuconsultation attribuconsultation attribuconsultation reason the free radicle office, sqrt
The pattern statute comprises a consultation output. Reason it to succor you debug your statute. You do referable attribuconsultation attribuconsultation attribuconsultation insufficiency to produce a consultation output.
Take trouble of the main portio of the completion primeval, then chaffer with the proper subjects
Don’t ceget to comprise input validation, ce pattern, the input must be scalar
office ANS = Radicle(a)
tol = 10^(-4);
x = abs(a)/2;
e = 99;
while (not(e < tol))
y = (x + abs(a)/x)/2;
e = abs((y-x)/y);
x = y;
%If circumstances satisfies
ANS = x;
falsity = abs(ANS – sqrt(a))
ANS = x*i;