One numerical method for calculating the cubic root of a number, 3 squareroot halley's method. the solution process starts by choosing a value x_1 as a first estimate of the solution. using this value, a second, more accurate value x_2 is calculated with, x_2 = x_1(x_1^3 + 2p)/(2x_1^3 + p) which is then used for calculating a third, still more accurate value x_3 and so on. the general equation for calculating the value of x_i + 1 from the value of x_i is x_i + 1 = x_i(x_i^3 + 2p)/(2x_i^3 + p). write a matlab program that calculates the cubic root of a number. in the program use x_1 = p for the first estimate of the solution. then, by using the general equation in a loop, calculate new, more accurate values. stop the looping when the estimated relative error e defined by e = |x_i + 1 - x_i/x_i) is smaller than 0.1. use the program to calculate: (a) 3√ 800 (b) 3√ 59071 (c) 3√ -31.55