This preview has intentionally blurred sections. Sign up to view the full version.
View Full Document
Unformatted text preview: CE 335 Solutions to Homework 2 13) function [x, err] = sqrt_apprx(a, err_desired) %[x, err] = sqrt_apprx(a, err_desired) %approximate the square root of a real number a %with fractional precision err_desired %using the divideandaverage method. %returns the approximate square root a %and the approximate fractional error err %take care of cases where a is zero or negative if a == 0 x = 0; err = 0; return elseif a < 0 a = a; mult = i; %will multiply the result for a by i else mult = 1; end %set initial value for the square root x = a; err = err_desired + 1; %ensure that iteration runs at least once %because the condition for the while loop %will start out as true %do iteration while err > err_desired x_old = x; x = 0.5 * (x + a/x); err = abs((x  x_old)/x); end %multiply result by i if input was negative x = mult * x; testing: > [x, err] = sqrt_apprx(0, 1E4)...
View
Full Document
 Spring '11
 Lybas
 Taylor Series, Err, fractional error err

Click to edit the document details