欧拉法解常微分方程组

用ode45求解器求解上述方程组，得到了理想中的曲线，

现在想用欧拉格式直接解出，得到的曲线却不是之前的曲线，能否帮忙看看问题出在哪里？下面是欧拉法的求解代码和运行结果
clear
clc
format long
yr1(1,:)=0;
dyr1(1,:)=0;
ddyr1(1,:)=0;
h=1;
a1=-64;
a2=-48;
a3=-12;
h=0.0001;%Step size
t=0:10;
for i=2:11
yr1(i,:)=dyr1(i-1,:)h+yr1(i-1,:);
dyr1(i,:)=ddyr1(i-1,:)h+dyr1(i-1,:);
if rem(i,12)>=0 & rem(i,12)<5
ddyr1(i,:)=(a1*(yr1(i-1,:)-0.1)+a2*(dyr1(i-1,:)-0)+a3*(ddyr1(i-1,:)-0))*h+ddyr1(i-1,:);
elseif rem(i,12)>=5 & rem(i,12)<12
ddyr1(i,:)=(a1(yr1(i-1,:)+0)+a2(dyr1(i-1,:)-0)+a3*(ddyr1(i-1,:)-0))*h+ddyr1(i-1,:);
end
end

figure(1)
plot(t,yr1(:,1),'b--','LineWidth',2),grid on
xlabel('Time (s)')
ylabel('position(mm)')
legend('$y_d$')