$(30$pts$)($Computer Simulation$)$ Simulate the Lorenz Attractor chaotic system using MATLAB. Run for $50$ seconds with all initial states equal to $0\mathrm{.}1,=10,r=28,b=\frac{8}{3}.$ Plot states versus time and make $3-$D plot of $x_1,x_2,x_3$ using plot $).($Include your MATLAB code$)$

$x_1^{}=-(x_1-x_2)$

$x_2^{}=m{x}_1-x_2-x_1{x}_3$

$x_3^{}=-b{x}_3+x_1{x}_2$