Исправление

@sovershenstvo · Регистрация: 18.02.2012

Author24 — интернет-сервис помощи студентам

Парашютист прыгает с некоторой высоты и летит, не открывая парашюта; на какой высоте (или через какое время) ему следует открыть парашют, чтобы иметь к моменту приземления безопасную скорость (не большую 10 м/с)?

Помогите исправить программу,проверку нужно оформить в мат пакете.
Преподователь требует чтобы начальная высота была задана например 1000,а безопасная колеблилась от 1000 до 150 ,а на графике мы должны показать какая скорость при этих безопасных высотах

Pascal

uses crt;
const
    vbez = 10.0;    plot = 1.303; m = 90.0; mu = 0.0182;    r = 0.3; e=0.01;
type
    vartype = array [1..3] of real;
    func = function(vars : vartype): real;
var
    res,s: real;
  v0, v1: vartype;
  H0, i:integer;
  tt:text;
  //c парашютом
function f1(v: vartype):real;
begin
    f1:= v[3];
end;
function f2(v: vartype):real;
begin
    f2 := (0.5 *0.4 * 2.7 * plot * v[3] * v[3]) / m;
end;
//без параш
function f3(v: vartype):real;
begin
     f3:= v[3];
end;
function f4(v: vartype):real;
begin
     f4:= (-m*9.8 + 6 * 3.14*mu*r*v[3]+0.5*0.55*1.7*0.4*plot*v[3]*v[3]) / m; //1.7м-средний рост человека,0,4-полуобхват грудной клетки выберем в качестве характерного расстояния
end;
function rez(var h1: vartype; var h2: vartype ): real;
begin
    while (abs(h2[2]) > abs(h1[2])) do begin
        h1[2] := h1[2] + e*f1(h1);
        h1[3] := h1[3] + e*f2(h1);
        h1[1] := h1[1] + e;
        h2[2] := h2[2] + e*f3(h2);
        h2[3] := h2[3] + e*f4(h2);
        h2[1] := h2[1] + e;
    end;
rez := h2[2];
end;
begin
write('Высота прыжка:');
readln(H0);s:=0.1;
assign(tt,'парашют1.txt');
rewrite(tt);
 
 
while (s < 3) do begin
 
    v0[1] := 0;
    v0[2] := 0;
    v0[3] := vbez;
    v1[1] := 0;
    v1[2] := H0;
    v1[3] := 0;
    res := rez(v0, v1);
    writeln(tt, H0,'   ',s,'   ',res:6:2,'   ',abs(v1[3]):6:2);
  s:= s+0.1;
  end;
close(tt);
end.

	04.12.2014, 18:31