Solving equations analytically and using Finite Difference method

Objective: clear all close all clc x=pi/3; dx= linspace(pi/4,pi/4000,50); analytical_solution=(x^3*(cos(x))-sin(x)*3*x^2)/x^6; for i=1:length(dx) central_differencing= ((sin(x+dx(i))/(x+dx(i))^3)-(sin(x-dx(i))/(x-dx(i))^3))/(2*dx(i)); second_order_error(i)=second_order_discretization(x,dx(i)); forward_differencing= ((sin(x+dx(i))/(x+dx(i))^3)-(sin(x)/x^3))/(dx(i));…