2y[n] = y[n-1] + x[n] / y[n-1]
If we use as input x[n] to this system (algorithm) a step function of amplitude A (i.e. x[n] = A u[n]), then y[n] will converge after several iterations to the square root of A.
• Write a MATLAB program that implements the above recursion to compute the square root of 16, 4, 5, and 3.
• How many iterations does it take to converge to the true value starting at y[-1] =0.5?
• Is this algorithm sensitive to the initial condition?