I need a MATLAB CODE to solve this PDE using the crank nicolson method with the constraints h $=0\mathrm{.}1$; $\%$ mesh size for x

k $=0\mathrm{.}05$; $\%$ mesh size for t

L $=1$; $\%$ spatial domain is $[0,$L$]$

T $=0\mathrm{.}95$; $\%$ time domain is $[0,$T$]$

m $=$ L$/$h; $\%$ number of subintervals in x

n $=$ T$/$k; $\%$ number of subintervals in t

x $=[0$:h:L$]$; $\%$ vectors containing grid points

t $=[0$:k:T$]$; $\%$ vectors containing grid points

$\%$ initial condition

u$\_0=$ x$.*(1-$ x$)$;$.$ Please create a MATLAB CODE for the solution!