gammazad := 1;
  t2 := 10;
  t1 := 0;
  n := 1000;
  h := (t2 - t1) / n;
  al := pi / 180;
  gamma[1] := 0;
  //gamma[1] := 0;
  gammaprev := gamma[1];
  wx := 0;    // omega x - угл скорость по крену
  wxl := 0;
  wxlpr := 0;
  //x1 := 0;
  x2 := -1;
  
  
  for i := 2 to 1000 do
  begin
    ksi := gammaprev-gammazad;
    wxpr := wx;
    alpr := al;
    //x1pr := x1;
    x2pr := x2;
    wxlpr := wxl;
    wxl := -b1*wx-b3*(mue*wx+ie*gammaprev+nue*ksi);
    wx := wxpr + h * wxl;  //
    x2 := x2pr + h * wx;
    //gamma[i] := x2 + gammazad;
    gamma[i] := gamma[i-1] + h * wx;
    //gamma[i] := x2;
    gammaprev := gamma[i];
    //gamma[i] := gamma[i]*180/pi;
  end;