Solve by MATLAB

18.17 Develop, debug, and test a program in either a high-level language or macro language of your choice to implement Newton's interpolating polynomial based on Fig. 18.7. SUBROUTINE NewtInt (, y, n, xi, yint, ea) LOCAL fddon DOFOR 1 = 0. fddo = y END DO DOFOR j = 1, n DOFOR 1 = 0. n-j fddij = (fdd +1.3-1 - fdd 1-1)/(X:+ 3 X:) END DO END DO xterm = 1 yint = fddoo DOFOR order = 1. n xterm = xterm* (xi Xorder-1 yint2 = yintorder-) + fddo, order * xterm eaarder1 = yint2 - yintorder- yintorder = yint2 END Order END Newt Int FIGURE 18.7 An algorithm for Newton's interpolating polynomial written in pseudocode