1. Modify the Fixpt function we built in class to use the Wegstein method. Upload the Wegstien_[your last name].m file to canvas.

function [x1,ea,iter]=fixpt(g,x0,es,maxit) % fixpt: fixed point root locator % This function deterimes the root of x=g(x) using fixed point iteration % The method is repeated until either the percent relative error (ea) is % equal to or less than es (default 1e-6) or the number of iterations % exceeds maxit (default 50). % % Inputs: % g: function for g(x) % x0: intial guess for x % es: stopping criterion for relative error % maxit: max number of iterations % Outputs % x1: calculated root % ea: relative error for x1 % iter: number of iterations calculated % Example: g=@(x) exp(-x); x0=0; [xr,ea,iter]=fixpt(g,x0,1e-6); if nargin = maxit) break end x0 = x1; end X = ['The root = ', num2str(x1), ' (ea = ',num2str(ea), '% in ',num2str(iter), ' iteration)']; disp(X) end