Week 3 - Speed Programming in Matlab Octave for non computer science engineers

Objective: #commentThe code seems to be an exact representation of the reference code. There is also a glitch % main bodyclose allclear allclct = [0 20]yo = [0 3][t, y] = ode45('ode_pendulum',t,yo)subplot(5,4,[13,14,15,16]);plot(t,y(:,1),'color','r');hold on;plot(t,y(:,2),'color','b');hold off;subplot(5,4,[17,18,19,20]);plot(y(:,1),y(:,2));xlabel('position_Radians')ylabel('Angular_velocity');ct=1;for…