Consider boundary value problem given below.

$y^{\text{'}\text{'}}+P(x)y^{\text{'}}+Q(x)y=f(x)$

which is defined for $x_{0}x{x}_f$ with boundary conditions $y(x_0)=y_0$ and $y(x_f)=y_f.$

Write a MATLAB code that allows you to define $P(x),Q(x)$ and $f(x),$ and finds numerical solutions of this differential equation using Finite Difference Method.

Your MATLAB code shall accept $x_0,x_f,y_0,y_f,$ and number of internals $n,$ and plots the results.

You have to submit your MATLAB source code, plus a report that shows the output of your code for a specific differential equation you choose, and three different $n$ values.