Помогите написать программу на MATLAB

[Maik2013]

Добрые люди пожалуйста помогите решит задачу на MATLAB
Для решения задачи

$$ \dfrac{\partial U}{\partial \partial t}=\dfrac{\partial^{2}U}{\partial x^{2}}, \quad U(x,0)=\Phi(x) $$

рассмотреть разностную схему

$$ \dfrac{u_{k}^{n+1}-u_{k}^{n}}{\tau}=\dfrac{u_{k+1}^{n}-2u_{k}^{n}+u_{k-1}^{n}}{h^{2}}, \quad u_{k}^{0}=\varphi_{k}. $$

Исследовать ее аппроксимацию, сходимость (при $2\tau/h^{2}\leq 1$)
и устойчивость (на функции $\varepsilon(-1)^{k}$)

Буду очень благодарен.