please solve using Matlab and post code.

In reference to Chapter 17 of the textbook Soil Mechanics by A. Verruijt, calculate using an explicit finite difference algorithm the dimensionless excess pore pressure W(Z,T) and the degree of consolidation U(T). The problem is governed by the following set of equations: = aw 22w Diffusion equation: T az? Initial conditions: W = 1,0 0 and Oy = 0, Z =0,7 > 0 aw az = = = Use AZ =1/20 and AT = AAZwith a =1/4. Plot the isochrones of W for T=0.1, 0.2, 0.5, 1, and 2 calculated numerically and compare to the analytical solution given by equation (16.11). The plot should be similar to Fig. 16.2 of the textbook. Plot the function U(T) computed numerically by integrating W and compare to the analytical solution given by equation (16.19) as well as with the small time response (square root of time). The plot should be similar to Fig. 16.3 of the textbook. In reference to Chapter 17 of the textbook Soil Mechanics by A. Verruijt, calculate using an explicit finite difference algorithm the dimensionless excess pore pressure W(Z,T) and the degree of consolidation U(T). The problem is governed by the following set of equations: = aw 22w Diffusion equation: T az? Initial conditions: W = 1,0 0 and Oy = 0, Z =0,7 > 0 aw az = = = Use AZ =1/20 and AT = AAZwith a =1/4. Plot the isochrones of W for T=0.1, 0.2, 0.5, 1, and 2 calculated numerically and compare to the analytical solution given by equation (16.11). The plot should be similar to Fig. 16.2 of the textbook. Plot the function U(T) computed numerically by integrating W and compare to the analytical solution given by equation (16.19) as well as with the small time response (square root of time). The plot should be similar to Fig. 16.3 of the textbook