dam.m

2005MCMA题程序

clc,clear
[T,Z]=ode45('damfail',[0 2.5],52.9)
figure(1)
plot(T,Z,'linewidth',2)
axis([0 2.5 0 60]);
grid on;
xlabel('Time-------t/hrs')
ylabel('Vertical Distance from Saluda River----Z/m')