... equations for D(k) and D2(k), as given respec-tively by Eqs. (4. 16) and (4.30), and with the initial conditions for the function and its derivative. The first elements for the y, D, and D2 arrays ... (2*a(k)*y(k-1)/dt+a(k)D(k-1)+u(k));D(k)=(2/dt)*(y(k)-y(k-1))-D(k-1);endplot(t,y,t,u,' ')In-Class ExercisePb. 4.37 Plot the amplitude of y, and its dephasing from u, as a function ofa for large t.Example 4.9Find the first-order ... LLCdt=(tfin-tin)/(N-1);u=sin(t);a=(1/(2*pi))*ones(1,N);b=ones(1,N);y(1)=0;D(1)=(1/a(1))*(u(1)-b(1)*y(1)); for k=2:Ny(k)=((2*a(k)/dt+b(k))^ (-1 ))* (2*a(k)*y(k-1)/dt+a(k)D(k-1)+u(k));D(k)=(2/dt)*(y(k)-y(k-1))-D(k-1);endplot(t,y,t,u,'...