Please include script code and output (Use MATLAB, do NOT do by hand)

Using MATLAB R2017b:

Write a matlab program that implements Heron's method to find the square root of 133,579 using 20 iterations of the algorithm. Heron's Algorithm for calculating square roots problem: Find sqrt(n)

1) Make a guess at the solution. (x₁)

2) x₂ = (x₁ + (n / x₁)) / 2

3) Repeat for x_3, x_4, x_{5, ...}