844
figure=plot.figure()
Kl1, Kl2=plot.plot(t, u, ‘b-‘, t, x0*cos(omega0*t), ‘r—‘)
figure.legend((Kl1, Kl2), (‘numerical’, ‘theoretical’), ‘upper left’)
plot.xlabel(‘t’)
plot.show()
plot.savefigure(‘tmp.png’)
The comparison between the numerical solutions performed in Figure 2 and the exact solutions
of the differential equation determined analytically is shown in one plane. Surprisingly, the incorrect
view of the numerical solution is due to a programming error or a problem with thechosen Euler
method?
Do'stlaringiz bilan baham: