Modify the following code to solve the diffusion equation with Neumann boundary conditions using finite differences.

For the code "[D,x]=cheb(N);", cheb() is one of the codng functions within chebfun, an open-source package for computing with functions to about 15-digit accuracy. In order to access the package,you will need to go to chebfun.org to download it.

Need help!! Please provide a step-by-step solution through the use of MALAB code along with a detailed explanation!!!! It is alright if you handwrite the step-by-step solution but I would prefer the solution to be shown in MATLAB code. This is very URGENT!!! PLEASE reply as soon as possible!! Thank you!!

"-t = u xx-ISK