Week 3 - Discretization basics

Objective: Aim:The code to plot the first order, second and fourth order error for a range of dxclear allclose allclcx=pi/3;dx=linspace(pi/5,pi/5000,20);% f(x)=sin(x)/x^3;% f'(x)=(cos(x)*x^3 - 3*x^2*sin(x))/x^6;analytical_derivative =(cos(x)*x^3 - 3*x^2*sin(x))/x^6;for i=1:length(dx) % first order approximation first_order_approximation(i)=(sin(x+dx(i))/(x+dx(i))^3-sin(x)/x^3)/dx(i);…