For discrete dynamical system: x(t+1) = x(t)^2 + c with x(0)=0. Estimate system for values of c=a+b*i in complex plane for a and b in interval [-1.5,1] with interval size of 0.005. Create c in matrix form using meshgrid. Evolve system up to t=20 for all values of c in parallel. Utilize image() to display, in complex plane, the values of c for which |x(20)| < 2. Looking for steps on being able to write this out and view image in MATLAB. Prefer code now, and it must be able to run without any requirement to change it.