Can you give me matlab code for this. A drained triaxial compression test is conducted on a lightly over$-$consolidated soil with

critical state parameters M $=0\mathrm{.}97,\backslash$kappa $=0\mathrm{.}04,\backslash$lambda $=0\mathrm{.}2,$ N $=2\mathrm{.}5,$ and Poisson$$s ratio $=$

$0\mathrm{.}15.$ The preconsolidation pressure $($p$0)$ is $100$ kPa and it is sheared from $60$ kPa. Use

modified Cam clay model $($associated flow rule$)$ to predict axial and volumetric strain.

Plot yield surface, critical state line and drained stress path in one figure. Also, plot

the variation of axial strain$($a$)$ vs$.$ deviatoric stress $($q$)$ and volumetric strain $($v$)$ vs

mean effective stress $($p$)$ in seperate figures. Use stress increment $\backslash$delta p $=1$ kPa