254
for
(i=0; i<n; i=i+1){
f1=fz(t0,p0,z0);
tv=t0+h/2.0;
pv=p0+h/2.0*z0;
zv=z0+h/2.0*f1;
f2=fz(tv,pv,zv);
t=t0+h;
p=p0+h*zv;
z=z0+h*f2;
file3<<t<<
"\t"
<<p<<
"\t"
<<z<<endl;
t0=t; p0=p; z0=z;
}
file3.close();
}
void
main()
{
cout<<
"n="
<<n<<endl;
eiler_prost();
eiler_usred();
eiler_centr();
}
Результат выполнения программы представлен на рис. 8.14.
Рис. 8.14. Результат выполнения программы
Результаты
расчета
методов
Эйлера
(простого
и
модифицированных) представлены на рис. 8.15 – рис. 8.17.
Рис. 8.15. Файл eiler_pr.xls, результат расчета простого метода Эйлера
