Series5.Add(p1[i]);

//Заполняютаблицу

StringGrid1.Cells[1,i+1]:= FloatToStrF(x[i],ffFixed,3,3);

StringGrid1.Cells[2,i+1]:= FloatToStrF(r[i],ffFixed,3,3);

StringGrid1.Cells[3,i+1]:= FloatToStrF(m[i],ffFixed,3,3);

StringGrid1.Cells[4,i+1]:= FloatToStrF(p[i],ffFixed,3,3);

StringGrid1.Cells[5,i+1]:= FloatToStrF(m1[i],ffFixed,3,3);

StringGrid1.Cells[6,i+1]:= FloatToStrF(p1[i],ffFixed,3,3);

except

end;

end;

end;

end.

Скришоты работы программы в сравнении с графиками MathCad:

Вывод

Проведенный в программе расчет на первых трех участках совпадает с данными, полученными при моделировании в MathCAD (с точностью до 10³ ), поэтому принимается, что программа расчета дает справедливые результаты для остальных участков обтекаемого тела.

Список литературы

1. Кузнецов В.И., Макаров В.В. Сверхзвуковое обтекание заостренных тел вращения: Учебное пособие.- Омск: Изд-во ОмГТУ, 1999г.- 80с.

2. Кузнецов В.И., Макаров В.В. Аэрогазодинамика в вопросах и задачах: практикум.- Омск: Изд-во ОмГТУ, 2002г.- 80с.