Experimental and modeling studies of desensitization of P2X3 receptors (Sokolova et al. 2006)

 Download zip file 
Help downloading and running models
Accession:117691
"The function of ATP-activated P2X3 receptors involved in pain sensation is modulated by desensitization, a phenomenon poorly understood. The present study used patch-clamp recording from cultured rat or mouse sensory neurons and kinetic modeling to clarify the properties of P2X3 receptor desensitization. ... Desensitization properties were well accounted for by a cyclic model in which receptors could be desensitized from either open or closed states. Recovery was assumed to be a multistate process with distinct kinetics dependent on the agonist-dependent dissociation rate from desensitized receptors. ... By using subthreshold concentrations of an HAD (high-affinity desensitization)-potent agonist, it might be possible to generate sustained inhibition of P2X3 receptors for controlling chronic pain."
Reference:
1 . Sokolova E, Skorinkin A, Moiseev I, Agrachev A, Nistri A, Giniatullin R (2006) Experimental and modeling studies of desensitization of P2X3 receptors. Mol Pharmacol 70:373-82 [PubMed]
Model Information (Click on a link to find other models with that property)
Model Type: Channel/Receptor;
Brain Region(s)/Organism:
Cell Type(s):
Channel(s):
Gap Junctions:
Receptor(s): Sensory Receptors;
Gene(s):
Transmitter(s):
Simulation Environment: Pascal/Delphi;
Model Concept(s): Ion Channel Kinetics; Nociception;
Implementer(s): Skorinkin, Andrei [askorink at yandex.ru];
Search NeuronDB for information about:  Sensory Receptors;
unit DS_especial_pas;

interface

uses
  Windows, Messages, SysUtils, Classes, Graphics, Controls, Forms, Dialogs,
  StdCtrls, ExtCtrls, Math, TeEngine, Series, TeeProcs, Chart, ComCtrls;

type
  TForm1 = class(TForm)
    Image1: TImage;
    Label1: TLabel;
    Image2: TImage;
    Image4: TImage;
    Label2: TLabel;
    Label4: TLabel;
    Chart1: TChart;
    Edit1: TEdit;
    Edit2: TEdit;
    Edit3: TEdit;
    Edit4: TEdit;
    Edit5: TEdit;
    Edit6: TEdit;
    Edit7: TEdit;
    Edit8: TEdit;
    Label6: TLabel;
    Edit9: TEdit;
    Label7: TLabel;
    Edit10: TEdit;
    Label8: TLabel;
    Edit11: TEdit;
    Label10: TLabel;
    Button1: TButton;
    Button2: TButton;
    Button3: TButton;
    Button4: TButton;
    Edit13: TEdit;
    Edit14: TEdit;
    Edit15: TEdit;
    Label11: TLabel;
    Label12: TLabel;
    Label19: TLabel;
    Edit24: TEdit;
    Label20: TLabel;
    Edit25: TEdit;
    Image5: TImage;
    Edit26: TEdit;
    Edit27: TEdit;
    Edit28: TEdit;
    Edit29: TEdit;
    Edit37: TEdit;
    Edit38: TEdit;
    Edit39: TEdit;
    Edit40: TEdit;
    CheckBox1: TCheckBox;
    Label24: TLabel;
    Edit32: TEdit;
    Image8: TImage;
    Label5: TLabel;
    Edit12: TEdit;
    Edit17: TEdit;
    Edit19: TEdit;
    Edit20: TEdit;
    Edit21: TEdit;
    Edit22: TEdit;
    Label30: TLabel;
    Edit31: TEdit;
    Label31: TLabel;
    Edit42: TEdit;
    Button5: TButton;
    Button6: TButton;
    Label13: TLabel;
    Edit18: TEdit;
    Image13: TImage;
    Label33: TLabel;
    Edit44: TEdit;
    ProgressBar1: TProgressBar;
    Label34: TLabel;
    Edit46: TEdit;
    Label35: TLabel;
    Edit47: TEdit;
    Edit48: TEdit;
    Button7: TButton;
    Button8: TButton;
    Series3: TPointSeries;
    Series1: TPointSeries;
    Series2: TPointSeries;
    OpenDialog1: TOpenDialog;
    SaveDialog1: TSaveDialog;
    Series4: TPointSeries;
    Series5: TPointSeries;
    Label36: TLabel;
    Edit49: TEdit;
    Edit50: TEdit;
    Label39: TLabel;
    Label40: TLabel;
    Label41: TLabel;
    Label42: TLabel;
    Label43: TLabel;
    Label44: TLabel;
    Label45: TLabel;
    Label46: TLabel;
    Label47: TLabel;
    Label48: TLabel;
    Image14: TImage;
    Label49: TLabel;
    Label50: TLabel;
    Image15: TImage;
    Label51: TLabel;
    Image16: TImage;
    Label52: TLabel;
    Label53: TLabel;
    Image17: TImage;
    Label54: TLabel;
    Image18: TImage;
    Label55: TLabel;
    Label56: TLabel;
    Image19: TImage;
    Label57: TLabel;
    Image20: TImage;
    Label58: TLabel;
    Edit51: TEdit;
    Label3: TLabel;
    Label9: TLabel;
    Edit53: TEdit;
    Edit54: TEdit;
    Image21: TImage;
    Label14: TLabel;
    Label15: TLabel;
    Edit55: TEdit;
    Edit56: TEdit;
    Edit57: TEdit;
    Label16: TLabel;
    Image6: TImage;
    Label17: TLabel;
    Edit58: TEdit;
    Edit59: TEdit;
    Label18: TLabel;
    Image22: TImage;
    Label21: TLabel;
    Edit60: TEdit;
    Edit61: TEdit;
    Label22: TLabel;
    Image23: TImage;
    Label23: TLabel;
    Edit62: TEdit;
    Image24: TImage;
    Label26: TLabel;
    Label27: TLabel;
    Edit63: TEdit;
    Edit64: TEdit;
    Image25: TImage;
    Label28: TLabel;
    Label29: TLabel;
    Edit65: TEdit;
    Edit66: TEdit;
    Edit67: TEdit;
    Label32: TLabel;
    Image26: TImage;
    Label37: TLabel;
    Edit68: TEdit;
    Edit69: TEdit;
    Label38: TLabel;
    Image27: TImage;
    Label59: TLabel;
    Edit70: TEdit;
    Edit71: TEdit;
    Label60: TLabel;
    Image28: TImage;
    Label61: TLabel;
    Edit72: TEdit;
    Label62: TLabel;
    Label63: TLabel;
    Label64: TLabel;
    Label65: TLabel;
    Edit73: TEdit;
    Edit74: TEdit;
    Edit75: TEdit;
    Edit76: TEdit;
    Label25: TLabel;
    Edit16: TEdit;
    Edit23: TEdit;
    Edit30: TEdit;
    Label66: TLabel;
    Edit33: TEdit;
    Edit34: TEdit;
    Edit35: TEdit;
    Edit36: TEdit;
    Label67: TLabel;
    Edit41: TEdit;
    Edit43: TEdit;
    Edit45: TEdit;
    Edit52: TEdit;
    Label68: TLabel;
    Edit77: TEdit;
    Edit78: TEdit;
    Edit79: TEdit;
    Label69: TLabel;
    Edit80: TEdit;
    Edit81: TEdit;
    Edit82: TEdit;
    Label70: TLabel;
    Edit83: TEdit;
    Edit84: TEdit;
    Edit85: TEdit;
    Edit86: TEdit;
    Label71: TLabel;
    Edit87: TEdit;
    Edit88: TEdit;
    Edit89: TEdit;
    Edit90: TEdit;
    Label72: TLabel;
    Label73: TLabel;
    Label74: TLabel;
    Label75: TLabel;
    Edit91: TEdit;
    Label76: TLabel;
    Edit92: TEdit;
    Label77: TLabel;
    Edit93: TEdit;
    Label78: TLabel;
    Image3: TImage;
    Edit94: TEdit;
    Edit95: TEdit;
    Image7: TImage;
    Label79: TLabel;
    Edit96: TEdit;
    Edit97: TEdit;
    Image9: TImage;
    Label80: TLabel;
    Edit98: TEdit;
    Edit99: TEdit;
    Image10: TImage;
    Label81: TLabel;
    Edit100: TEdit;
    Edit101: TEdit;
    procedure Button3Click(Sender: TObject);
    procedure Button1Click(Sender: TObject);
    procedure Button4Click(Sender: TObject);
    procedure Button2Click(Sender: TObject);
    procedure Button5Click(Sender: TObject);
    procedure Button6Click(Sender: TObject);
    procedure Button8Click(Sender: TObject);
    procedure Button7Click(Sender: TObject);
  private
    { Private declarations }
  public
    { Public declarations }
  end;

var
  Form1: TForm1;
  N,NJ:integer;
  yh:array[0..27] of real;
  ar,a2r,a3r,a3ro,tx,
  rr,a3do,a3d,a2d,ad,d,
  bd,br,abd,abr,b2d,b2r,
  a2bd,a2br,ab2d,ab2r,b3d,b3r,
  a2bdo,a2bro,ab2do,ab2ro,b3do,b3ro,
  sumo,sumr,sumd:array[0..20000] of real;
  am,ta1,ta2,rt,tt0,tt1,tt2,agon,ag,tb0,tb1,tb2,bgon,bg,
  rish,arish,a2rish,a3rish,a3roish,
  dish,adish,a2dish,a3dish,a3doish,
  bdish,brish,abdish,abrish,b2dish,b2rish,
  a2bdish,a2brish,ab2dish,ab2rish,b3dish,b3rish,
  a2bdoish,a2broish,ab2doish,ab2roish,b3doish,b3roish,
  k01,k10,k02,k20,k13,k31,k14,k41,
  k24,k42,k25,k52,k36,k63,k37,k73,k47,k74,
  k48,k84,k58,k85,k59,k95,k610,k106,
  k711,k117,k812,k128,k913,k139,
  k1415,k1514,k1416,k1614,k1517,k1715,k1518,k1815,
  k1618,k1816,k1619,k1916,k1720,k2017,k1721,k2117,k1821,k2118,
  k1822,k2218,k1922,k2219,k1923,k2319,k2024,k2420,
  k2125,k2521,k2226,k2622,k2327,k2723,
  k014,k140,k115,k151,k216,k162,k317,k173,k418,k184,
  k519,k195,k620,k206,k721,k217,k822,k228,k923,k239,
  k1024,k2410,k1125,k2511,k1226,k2612,k1327,k2713,
  tauf,taus:real;

implementation

{$R *.DFM}


function Sign(xs:real):real;
begin
   if xs=0 then Sign:=0 else
   if xs>0 then Sign:=1 else Sign:=-1;
end;

function Vspl(xt,x0,x1,xa:real):real;
var xf1,xf2,taf,tas,pok:real;
begin
   pok:=1.0/6.0;
   taf:=tauf/Ln((1-Power(0.1,pok))/(1-Power(0.9,pok)));
   tas:=taus/Ln((1-Power(0.1,pok))/(1-Power(0.9,pok)));
   if xt<x0 then xf1:=0.0
   else xf1:=Power(1.0-exp((x0-xt)/taf),6.0);
   if xt<(x0+x1) then xf2:=1.0
   else xf2:=Power(1.0-(1.0-exp((x0+x1-xt)/tas)),6.0);
   Vspl:=xf1*xf2*xa;
end;

function so1(so1so0,so1so1,so1so3,so1so4,so1so15:real):real;
begin
so1:=k01*ag*so1so0-k10*so1so1-k13*ag*so1so1+k31*so1so3-k14*bg*so1so1+k41*so1so4
-k115*so1so1+k151*so1so15;
end;


function so2(so2so0,so2so2,so2so4,so2so5,so2so16:real):real;
begin
so2:=k02*bg*so2so0-k20*so2so2-k24*ag*so2so2+k42*so2so4-k25*bg*so2so2+k52*so2so5
-k216*so2so2+k162*so2so16;
end;


function so3(so3so1,so3so3,so3so6,so3so7,so3so17:real):real;
begin
so3:=k13*ag*so3so1-k31*so3so3-k36*ag*so3so3+k63*so3so6-k37*bg*so3so3+k73*so3so7
-k317*so3so3+k173*so3so17;
end;

function so4(so4so1,so4so2,so4so4,so4so7,so4so8,so4so18:real):real;
begin
so4:=k14*bg*so4so1-k41*so4so4+k24*ag*so4so2-k42*so4so4-k47*ag*so4so4+k74*so4so7
-k48*bg*so4so4+k84*so4so8-k418*so4so4+k184*so4so18;
end;

function so5(so5so2,so5so5,so5so8,so5so9,so5so19:real):real;
begin
so5:=k25*bg*so5so2-k52*so5so5-k58*ag*so5so5+k85*so5so8-k59*bg*so5so5+k95*so5so9
-k519*so5so5+k195*so5so19;
end;

function so6(so6so3,so6so6,so6so10,so6so20:real):real;
begin
so6:=k36*ag*so6so3-k63*so6so6-k610*so6so6+k106*so6so10
-k620*so6so6+k206*so6so20;
end;

function so7(so7so3,so7so4,so7so7,so7so11,so7so21:real):real;
begin
so7:=k37*bg*so7so3-k73*so7so7+k47*ag*so7so4-k74*so7so7-k711*so7so7+k117*so7so11
-k721*so7so7+k217*so7so21;
end;

function so8(so8so4,so8so5,so8so8,so8so12,so8so22:real):real;
begin
so8:=k48*bg*so8so4-k84*so8so8+k58*ag*so8so5-k85*so8so8-k812*so8so8+k128*so8so12
-k822*so8so8+k228*so8so22;
end;

function so9(so9so5,so9so9,so9so13,so9so23:real):real;
begin
so9:=k59*bg*so9so5-k95*so9so9-k913*so9so9+k139*so9so13
-k923*so9so9+k239*so9so23;
end;

function so10(so10so6,so10so10,so10so24:real):real;
begin
so10:=k610*so10so6-k106*so10so10-k1024*so10so10+k2410*so10so24;
end;

function so11(so11so7,so11so11,so11so25:real):real;
begin
so11:=k711*so11so7-k117*so11so11-k1125*so11so11+k2511*so11so25;
end;

function so12(so12so8,so12so12,so12so26:real):real;
begin
so12:=k812*so12so8-k128*so12so12-k1226*so12so12+k2612*so12so26;
end;

function so13(so13so9,so13so13,so13so27:real):real;
begin
so13:=k913*so13so9-k139*so13so13-k1327*so13so13+k2713*so13so27;
end;

function so14(so14so0,so14so14,so14so15,so14so16:real):real;
begin
so14:=k014*so14so0-k140*so14so14-k1415*ag*so14so14+k1514*so14so15
-k1416*bg*so14so14+k1614*so14so16;
end;

function so15(so15so14,so15so15,so15so17,so15so18,so15so1:real):real;
begin
so15:=k1415*ag*so15so14-k1514*so15so15-k1517*ag*so15so15+k1715*so15so17
-k1518*bg*so15so15+k1815*so15so18-k151*so15so15+k115*so15so1;
end;


function so16(so16so14,so16so16,so16so18,so16so19,so16so2:real):real;
begin
so16:=k1416*bg*so16so14-k1614*so16so16-k1618*ag*so16so16+k1816*so16so18
-k1619*bg*so16so16+k1916*so16so19-k162*so16so16+k216*so16so2;
end;


function so17(so17so15,so17so17,so17so20,so17so21,so17so3:real):real;
begin
so17:=k1517*ag*so17so15-k1715*so17so17-k1720*ag*so17so17+k2017*so17so20
-k1721*bg*so17so17+k2117*so17so21-k173*so17so17+k317*so17so3;
end;

function so18(so18so15,so18so16,so18so18,so18so21,so18so22,so18so4:real):real;
begin
so18:=k1518*bg*so18so15-k1815*so18so18+k1618*ag*so18so16-k1816*so18so18
-k1821*ag*so18so18+k2118*so18so21-k1822*bg*so18so18+k2218*so18so22
-k184*so18so18+k418*so18so4;
end;

function so19(so19so16,so19so19,so19so22,so19so23,so19so5:real):real;
begin
so19:=k1619*bg*so19so16-k1916*so19so19-k1922*ag*so19so19+k2219*so19so22
-k1923*bg*so19so19+k2319*so19so23-k195*so19so19+k519*so19so5;
end;

function so20(so20so17,so20so20,so20so24,so20so6:real):real;
begin
so20:=k1720*ag*so20so17-k2017*so20so20-k2024*so20so20+k2420*so20so24
-k206*so20so20+k620*so20so6;
end;

function so21(so21so17,so21so18,so21so21,so21so25,so21so7:real):real;
begin
so21:=k1721*bg*so21so17-k2117*so21so21+k1821*ag*so21so18-k2118*so21so21
-k2125*so21so21+k2521*so21so25-k217*so21so21+k721*so21so7;
end;

function so22(so22so18,so22so19,so22so22,so22so26,so22so8:real):real;
begin
so22:=k1822*bg*so22so18-k2218*so22so22+k1922*ag*so22so19-k2219*so22so22
-k2226*so22so22+k2622*so22so26-k228*so22so22+k822*so22so8;
end;

function so23(so23so19,so23so23,so23so27,so23so9:real):real;
begin
so23:=k1923*bg*so23so19-k2319*so23so23-k2327*so23so23+k2723*so23so27
-k239*so23so23+k923*so23so9;
end;

function so24(so24so20,so24so24,so24so10:real):real;
begin
so24:=k2024*so24so20-k2420*so24so24-k2410*so24so24+k1024*so24so10;
end;

function so25(so25so21,so25so25,so25so11:real):real;
begin
so25:=k2125*so25so21-k2521*so25so25-k2511*so25so25+k1125*so25so11;
end;

function so26(so26so22,so26so26,so26so12:real):real;
begin
so26:=k2226*so26so22-k2622*so26so26-k2612*so26so26+k1226*so26so12;
end;

function so27(so27so23,so27so27,so27so13:real):real;
begin
so27:=k2327*so27so23-k2723*so27so27-k2713*so27so27+k1327*so27so13;
end;


procedure TForm1.Button3Click(Sender: TObject);
var iii,jjj:integer;
agu,bgu:real;
fou:TextFile;
ti,si:string;
begin
ti:=TimeToStr(Time);
jjj:=Length(ti)-1;
for iii:=0 to jjj do if ti[iii]=':' then ti[iii]:='_';
si:='DS_especial_'+ti+'_r.dat';
AssignFile(fou,si);
Rewrite(fou);
Writeln(fou,'    t       r      ar    br    a2r    abr   b2r   a3r    a2br     ab2r     b3r');
for iii:=1 to NJ do
Writeln(fou,tx[iii]:9:3,rr[iii]:7:3,ar[iii]:7:3,br[iii]:7:3,a2r[iii]:7:3,
abr[iii]:7:3,b2r[iii]:7:3,a3r[iii]:7:3,a2br[iii]:7:3,ab2r[iii]:7:3,b3r[iii]:7:3);
Writeln(fou);
if CheckBox1.Checked=false then
begin
Writeln(fou,'ag = ',agon:9:6,'   t0 = ',tt0:9:6,'   t1 = ',tt1:9:6,'   t2 = ',tt2:9:6);
Writeln(fou,'bg = ',bgon:9:6,'   t0 = ',tb0:9:6,'   t1 = ',tb1:9:6,'   t2 = ',tb2:9:6);
Writeln(fou,'am = ',am:9:6,'  tau1 = ',ta1:12:6,'  tau2 = ',ta2:12:6,'  rt = ',rt:12:6);
Writeln(fou,'tauf = ',tauf:9:6,'  taus = ',taus:9:6,'  Parametr file: ',Edit48.Text);
end;
Writeln(fou);
Writeln(fou);
Writeln(fou,'k01 = ',k01:13:5,'      k10 = ',k10:13:5);
Writeln(fou,'k02 = ',k02:13:5,'      k20 = ',k20:13:5);
Writeln(fou,'k13 = ',k13:13:5,'      k31 = ',k31:13:5);
Writeln(fou,'k14 = ',k14:13:5,'      k41 = ',k41:13:5);
Writeln(fou,'k24 = ',k24:13:5,'      k42 = ',k42:13:5);
Writeln(fou,'k25 = ',k25:13:5,'      k52 = ',k52:13:5);
Writeln(fou,'k36 = ',k36:13:5,'      k63 = ',k63:13:5);
Writeln(fou,'k37 = ',k37:13:5,'      k73 = ',k73:13:5);
Writeln(fou,'k47 = ',k47:13:5,'      k74 = ',k74:13:5);
Writeln(fou,'k48 = ',k48:13:5,'      k84 = ',k84:13:5);
Writeln(fou,'k58 = ',k58:13:5,'      k85 = ',k85:13:5);
Writeln(fou,'k59 = ',k59:13:5,'      k95 = ',k95:13:5);
Writeln(fou);
Writeln(fou,'k610 = ',k610:13:5,'      k106 = ',k106:13:5);
Writeln(fou,'k711 = ',k711:13:5,'      k117 = ',k117:13:5);
Writeln(fou,'k812 = ',k812:13:5,'      k128 = ',k128:13:5);
Writeln(fou,'k913 = ',k913:13:5,'      k139 = ',k139:13:5);
Writeln(fou);
Writeln(fou);
Writeln(fou,'k1415 = ',k1415:13:5,'      k1514 = ',k1514:13:5);
Writeln(fou,'k1416 = ',k1416:13:5,'      k1614 = ',k1614:13:5);
Writeln(fou,'k1517 = ',k1517:13:5,'      k1715 = ',k1715:13:5);
Writeln(fou,'k1518 = ',k1518:13:5,'      k1815 = ',k1815:13:5);
Writeln(fou,'k1618 = ',k1618:13:5,'      k1816 = ',k1816:13:5);
Writeln(fou,'k1619 = ',k1619:13:5,'      k1916 = ',k1916:13:5);
Writeln(fou,'k1720 = ',k1720:13:5,'      k2017 = ',k2017:13:5);
Writeln(fou,'k1721 = ',k1721:13:5,'      k2117 = ',k2117:13:5);
Writeln(fou,'k1821 = ',k1821:13:5,'      k2118 = ',k2118:13:5);
Writeln(fou,'k1822 = ',k1822:13:5,'      k2218 = ',k2218:13:5);
Writeln(fou,'k1922 = ',k1922:13:5,'      k2219 = ',k2219:13:5);
Writeln(fou,'k1923 = ',k1923:13:5,'      k2319 = ',k2319:13:5);
Writeln(fou);
Writeln(fou,'k2024 = ',k2024:13:5,'      k2420 = ',k2420:13:5);
Writeln(fou,'k2125 = ',k2125:13:5,'      k2521 = ',k2521:13:5);
Writeln(fou,'k2226 = ',k2226:13:5,'      k2622 = ',k2622:13:5);
Writeln(fou,'k2327 = ',k2327:13:5,'      k2723 = ',k2723:13:5);
Writeln(fou);
Writeln(fou);
Writeln(fou,'k014 = ',k014:13:5,'      k140 = ',k140:13:5);
Writeln(fou,'k115 = ',k115:13:5,'      k151 = ',k151:13:5);
Writeln(fou,'k216 = ',k216:13:5,'      k162 = ',k162:13:5);
Writeln(fou,'k317 = ',k317:13:5,'      k173 = ',k173:13:5);
Writeln(fou,'k418 = ',k418:13:5,'      k184 = ',k184:13:5);
Writeln(fou,'k519 = ',k519:13:5,'      k195 = ',k195:13:5);
Writeln(fou,'k620 = ',k620:13:5,'      k206 = ',k206:13:5);
Writeln(fou,'k721 = ',k721:13:5,'      k217 = ',k217:13:5);
Writeln(fou,'k822 = ',k822:13:5,'      k228 = ',k228:13:5);
Writeln(fou,'k923 = ',k923:13:5,'      k239 = ',k239:13:5);
Writeln(fou);
Writeln(fou,'k1024 = ',k1024:13:5,'      k2410 = ',k2410:13:5);
Writeln(fou,'k1125 = ',k1125:13:5,'      k2511 = ',k2511:13:5);
Writeln(fou,'k1226 = ',k1226:13:5,'      k2612 = ',k2612:13:5);
Writeln(fou,'k1327 = ',k1327:13:5,'      k2713 = ',k2713:13:5);
Writeln(fou);
CloseFile(fou);
si:='DS_especial_'+ti+'_d.dat';
AssignFile(fou,si);
Rewrite(fou);
Writeln(fou,'    t       d      ad    bd    a2d    abd   b2d   a3d    a2bd     ab2d     b3d');
for iii:=1 to NJ do
Writeln(fou,tx[iii]:9:3,d[iii]:7:3,ad[iii]:7:3,bd[iii]:7:3,a2d[iii]:7:3,
abd[iii]:7:3,b2d[iii]:7:3,a3d[iii]:7:3,a2bd[iii]:7:3,ab2d[iii]:7:3,b3d[iii]:7:3);
Writeln(fou);
if CheckBox1.Checked=false then
begin
Writeln(fou,'ag = ',agon:9:6,'   t0 = ',tt0:9:6,'   t1 = ',tt1:9:6,'   t2 = ',tt2:9:6);
Writeln(fou,'bg = ',bgon:9:6,'   t0 = ',tb0:9:6,'   t1 = ',tb1:9:6,'   t2 = ',tb2:9:6);
Writeln(fou,'am = ',am:9:6,'  tau1 = ',ta1:12:6,'  tau2 = ',ta2:12:6,'  rt = ',rt:12:6);
Writeln(fou,'tauf = ',tauf:9:6,'  taus = ',taus:9:6,'  Parametr file: ',Edit48.Text);
end;
Writeln(fou);
Writeln(fou);
Writeln(fou,'k01 = ',k01:13:5,'      k10 = ',k10:13:5);
Writeln(fou,'k02 = ',k02:13:5,'      k20 = ',k20:13:5);
Writeln(fou,'k13 = ',k13:13:5,'      k31 = ',k31:13:5);
Writeln(fou,'k14 = ',k14:13:5,'      k41 = ',k41:13:5);
Writeln(fou,'k24 = ',k24:13:5,'      k42 = ',k42:13:5);
Writeln(fou,'k25 = ',k25:13:5,'      k52 = ',k52:13:5);
Writeln(fou,'k36 = ',k36:13:5,'      k63 = ',k63:13:5);
Writeln(fou,'k37 = ',k37:13:5,'      k73 = ',k73:13:5);
Writeln(fou,'k47 = ',k47:13:5,'      k74 = ',k74:13:5);
Writeln(fou,'k48 = ',k48:13:5,'      k84 = ',k84:13:5);
Writeln(fou,'k58 = ',k58:13:5,'      k85 = ',k85:13:5);
Writeln(fou,'k59 = ',k59:13:5,'      k95 = ',k95:13:5);
Writeln(fou);
Writeln(fou,'k610 = ',k610:13:5,'      k106 = ',k106:13:5);
Writeln(fou,'k711 = ',k711:13:5,'      k117 = ',k117:13:5);
Writeln(fou,'k812 = ',k812:13:5,'      k128 = ',k128:13:5);
Writeln(fou,'k913 = ',k913:13:5,'      k139 = ',k139:13:5);
Writeln(fou);
Writeln(fou);
Writeln(fou,'k1415 = ',k1415:13:5,'      k1514 = ',k1514:13:5);
Writeln(fou,'k1416 = ',k1416:13:5,'      k1614 = ',k1614:13:5);
Writeln(fou,'k1517 = ',k1517:13:5,'      k1715 = ',k1715:13:5);
Writeln(fou,'k1518 = ',k1518:13:5,'      k1815 = ',k1815:13:5);
Writeln(fou,'k1618 = ',k1618:13:5,'      k1816 = ',k1816:13:5);
Writeln(fou,'k1619 = ',k1619:13:5,'      k1916 = ',k1916:13:5);
Writeln(fou,'k1720 = ',k1720:13:5,'      k2017 = ',k2017:13:5);
Writeln(fou,'k1721 = ',k1721:13:5,'      k2117 = ',k2117:13:5);
Writeln(fou,'k1821 = ',k1821:13:5,'      k2118 = ',k2118:13:5);
Writeln(fou,'k1822 = ',k1822:13:5,'      k2218 = ',k2218:13:5);
Writeln(fou,'k1922 = ',k1922:13:5,'      k2219 = ',k2219:13:5);
Writeln(fou,'k1923 = ',k1923:13:5,'      k2319 = ',k2319:13:5);
Writeln(fou);
Writeln(fou,'k2024 = ',k2024:13:5,'      k2420 = ',k2420:13:5);
Writeln(fou,'k2125 = ',k2125:13:5,'      k2521 = ',k2521:13:5);
Writeln(fou,'k2226 = ',k2226:13:5,'      k2622 = ',k2622:13:5);
Writeln(fou,'k2327 = ',k2327:13:5,'      k2723 = ',k2723:13:5);
Writeln(fou);
Writeln(fou);
Writeln(fou,'k014 = ',k014:13:5,'      k140 = ',k140:13:5);
Writeln(fou,'k115 = ',k115:13:5,'      k151 = ',k151:13:5);
Writeln(fou,'k216 = ',k216:13:5,'      k162 = ',k162:13:5);
Writeln(fou,'k317 = ',k317:13:5,'      k173 = ',k173:13:5);
Writeln(fou,'k418 = ',k418:13:5,'      k184 = ',k184:13:5);
Writeln(fou,'k519 = ',k519:13:5,'      k195 = ',k195:13:5);
Writeln(fou,'k620 = ',k620:13:5,'      k206 = ',k206:13:5);
Writeln(fou,'k721 = ',k721:13:5,'      k217 = ',k217:13:5);
Writeln(fou,'k822 = ',k822:13:5,'      k228 = ',k228:13:5);
Writeln(fou,'k923 = ',k923:13:5,'      k239 = ',k239:13:5);
Writeln(fou);
Writeln(fou,'k1024 = ',k1024:13:5,'      k2410 = ',k2410:13:5);
Writeln(fou,'k1125 = ',k1125:13:5,'      k2511 = ',k2511:13:5);
Writeln(fou,'k1226 = ',k1226:13:5,'      k2612 = ',k2612:13:5);
Writeln(fou,'k1327 = ',k1327:13:5,'      k2713 = ',k2713:13:5);
Writeln(fou);
CloseFile(fou);
si:='DS_especial_'+ti+'_o.dat';
AssignFile(fou,si);
Rewrite(fou);
Writeln(fou,'    t    a3ro   a2bro   ab2ro   b3ro    a3do   a2bdo   ab2do   b3do');
for iii:=1 to NJ do
Writeln(fou,tx[iii]:9:3,a3ro[iii]:7:3,a2bro[iii]:7:3,ab2ro[iii]:7:3,b3ro[iii]:7:3,
a3do[iii]:7:3,a2bdo[iii]:7:3,ab2do[iii]:7:3,b3do[iii]:7:3);
Writeln(fou);
if CheckBox1.Checked=false then
begin
Writeln(fou,'ag = ',agon:9:6,'   t0 = ',tt0:9:6,'   t1 = ',tt1:9:6,'   t2 = ',tt2:9:6);
Writeln(fou,'bg = ',bgon:9:6,'   t0 = ',tb0:9:6,'   t1 = ',tb1:9:6,'   t2 = ',tb2:9:6);
Writeln(fou,'am = ',am:9:6,'  tau1 = ',ta1:12:6,'  tau2 = ',ta2:12:6,'  rt = ',rt:12:6);
Writeln(fou,'tauf = ',tauf:9:6,'  taus = ',taus:9:6,'  Parameter file: ',Edit48.Text);
end;
Writeln(fou);
Writeln(fou);
Writeln(fou,'k01 = ',k01:13:5,'      k10 = ',k10:13:5);
Writeln(fou,'k02 = ',k02:13:5,'      k20 = ',k20:13:5);
Writeln(fou,'k13 = ',k13:13:5,'      k31 = ',k31:13:5);
Writeln(fou,'k14 = ',k14:13:5,'      k41 = ',k41:13:5);
Writeln(fou,'k24 = ',k24:13:5,'      k42 = ',k42:13:5);
Writeln(fou,'k25 = ',k25:13:5,'      k52 = ',k52:13:5);
Writeln(fou,'k36 = ',k36:13:5,'      k63 = ',k63:13:5);
Writeln(fou,'k37 = ',k37:13:5,'      k73 = ',k73:13:5);
Writeln(fou,'k47 = ',k47:13:5,'      k74 = ',k74:13:5);
Writeln(fou,'k48 = ',k48:13:5,'      k84 = ',k84:13:5);
Writeln(fou,'k58 = ',k58:13:5,'      k85 = ',k85:13:5);
Writeln(fou,'k59 = ',k59:13:5,'      k95 = ',k95:13:5);
Writeln(fou);
Writeln(fou,'k610 = ',k610:13:5,'      k106 = ',k106:13:5);
Writeln(fou,'k711 = ',k711:13:5,'      k117 = ',k117:13:5);
Writeln(fou,'k812 = ',k812:13:5,'      k128 = ',k128:13:5);
Writeln(fou,'k913 = ',k913:13:5,'      k139 = ',k139:13:5);
Writeln(fou);
Writeln(fou);
Writeln(fou,'k1415 = ',k1415:13:5,'      k1514 = ',k1514:13:5);
Writeln(fou,'k1416 = ',k1416:13:5,'      k1614 = ',k1614:13:5);
Writeln(fou,'k1517 = ',k1517:13:5,'      k1715 = ',k1715:13:5);
Writeln(fou,'k1518 = ',k1518:13:5,'      k1815 = ',k1815:13:5);
Writeln(fou,'k1618 = ',k1618:13:5,'      k1816 = ',k1816:13:5);
Writeln(fou,'k1619 = ',k1619:13:5,'      k1916 = ',k1916:13:5);
Writeln(fou,'k1720 = ',k1720:13:5,'      k2017 = ',k2017:13:5);
Writeln(fou,'k1721 = ',k1721:13:5,'      k2117 = ',k2117:13:5);
Writeln(fou,'k1821 = ',k1821:13:5,'      k2118 = ',k2118:13:5);
Writeln(fou,'k1822 = ',k1822:13:5,'      k2218 = ',k2218:13:5);
Writeln(fou,'k1922 = ',k1922:13:5,'      k2219 = ',k2219:13:5);
Writeln(fou,'k1923 = ',k1923:13:5,'      k2319 = ',k2319:13:5);
Writeln(fou);
Writeln(fou,'k2024 = ',k2024:13:5,'      k2420 = ',k2420:13:5);
Writeln(fou,'k2125 = ',k2125:13:5,'      k2521 = ',k2521:13:5);
Writeln(fou,'k2226 = ',k2226:13:5,'      k2622 = ',k2622:13:5);
Writeln(fou,'k2327 = ',k2327:13:5,'      k2723 = ',k2723:13:5);
Writeln(fou);
Writeln(fou);
Writeln(fou,'k014 = ',k014:13:5,'      k140 = ',k140:13:5);
Writeln(fou,'k115 = ',k115:13:5,'      k151 = ',k151:13:5);
Writeln(fou,'k216 = ',k216:13:5,'      k162 = ',k162:13:5);
Writeln(fou,'k317 = ',k317:13:5,'      k173 = ',k173:13:5);
Writeln(fou,'k418 = ',k418:13:5,'      k184 = ',k184:13:5);
Writeln(fou,'k519 = ',k519:13:5,'      k195 = ',k195:13:5);
Writeln(fou,'k620 = ',k620:13:5,'      k206 = ',k206:13:5);
Writeln(fou,'k721 = ',k721:13:5,'      k217 = ',k217:13:5);
Writeln(fou,'k822 = ',k822:13:5,'      k228 = ',k228:13:5);
Writeln(fou,'k923 = ',k923:13:5,'      k239 = ',k239:13:5);
Writeln(fou);
Writeln(fou,'k1024 = ',k1024:13:5,'      k2410 = ',k2410:13:5);
Writeln(fou,'k1125 = ',k1125:13:5,'      k2511 = ',k2511:13:5);
Writeln(fou,'k1226 = ',k1226:13:5,'      k2612 = ',k2612:13:5);
Writeln(fou,'k1327 = ',k1327:13:5,'      k2713 = ',k2713:13:5);
Writeln(fou);
CloseFile(fou);
si:='DS_especial_'+ti+'_s.dat';
AssignFile(fou,si);
Rewrite(fou);
Writeln(fou,'    t       r       o       d       ag       bg');
for iii:=1 to NJ do
begin
agu:=vspl(tx[iii],tt0,tt1,agon);
bgu:=vspl(tx[iii],tb0,tb1,bgon);
Writeln(fou,tx[iii]:9:3,sumr[iii]:7:3,sumo[iii]:7:3,sumd[iii]:7:3,
agu:10:6,bgu:10:6);
end;
Writeln(fou);
if CheckBox1.Checked=false then
begin
Writeln(fou,'ag = ',agon:9:6,'   t0 = ',tt0:9:6,'   t1 = ',tt1:9:6,'   t2 = ',tt2:9:6);
Writeln(fou,'bg = ',bgon:9:6,'   t0 = ',tb0:9:6,'   t1 = ',tb1:9:6,'   t2 = ',tb2:9:6);
Writeln(fou,'am = ',am:9:6,'  tau1 = ',ta1:12:6,'  tau2 = ',ta2:12:6,'  rt = ',rt:12:6);
Writeln(fou,'tauf = ',tauf:9:6,'  taus = ',taus:9:6,'  Parametr file: ',Edit48.Text);
end;
Writeln(fou);
Writeln(fou);
Writeln(fou,'k01 = ',k01:13:5,'      k10 = ',k10:13:5);
Writeln(fou,'k02 = ',k02:13:5,'      k20 = ',k20:13:5);
Writeln(fou,'k13 = ',k13:13:5,'      k31 = ',k31:13:5);
Writeln(fou,'k14 = ',k14:13:5,'      k41 = ',k41:13:5);
Writeln(fou,'k24 = ',k24:13:5,'      k42 = ',k42:13:5);
Writeln(fou,'k25 = ',k25:13:5,'      k52 = ',k52:13:5);
Writeln(fou,'k36 = ',k36:13:5,'      k63 = ',k63:13:5);
Writeln(fou,'k37 = ',k37:13:5,'      k73 = ',k73:13:5);
Writeln(fou,'k47 = ',k47:13:5,'      k74 = ',k74:13:5);
Writeln(fou,'k48 = ',k48:13:5,'      k84 = ',k84:13:5);
Writeln(fou,'k58 = ',k58:13:5,'      k85 = ',k85:13:5);
Writeln(fou,'k59 = ',k59:13:5,'      k95 = ',k95:13:5);
Writeln(fou);
Writeln(fou,'k610 = ',k610:13:5,'      k106 = ',k106:13:5);
Writeln(fou,'k711 = ',k711:13:5,'      k117 = ',k117:13:5);
Writeln(fou,'k812 = ',k812:13:5,'      k128 = ',k128:13:5);
Writeln(fou,'k913 = ',k913:13:5,'      k139 = ',k139:13:5);
Writeln(fou);
Writeln(fou);
Writeln(fou,'k1415 = ',k1415:13:5,'      k1514 = ',k1514:13:5);
Writeln(fou,'k1416 = ',k1416:13:5,'      k1614 = ',k1614:13:5);
Writeln(fou,'k1517 = ',k1517:13:5,'      k1715 = ',k1715:13:5);
Writeln(fou,'k1518 = ',k1518:13:5,'      k1815 = ',k1815:13:5);
Writeln(fou,'k1618 = ',k1618:13:5,'      k1816 = ',k1816:13:5);
Writeln(fou,'k1619 = ',k1619:13:5,'      k1916 = ',k1916:13:5);
Writeln(fou,'k1720 = ',k1720:13:5,'      k2017 = ',k2017:13:5);
Writeln(fou,'k1721 = ',k1721:13:5,'      k2117 = ',k2117:13:5);
Writeln(fou,'k1821 = ',k1821:13:5,'      k2118 = ',k2118:13:5);
Writeln(fou,'k1822 = ',k1822:13:5,'      k2218 = ',k2218:13:5);
Writeln(fou,'k1922 = ',k1922:13:5,'      k2219 = ',k2219:13:5);
Writeln(fou,'k1923 = ',k1923:13:5,'      k2319 = ',k2319:13:5);
Writeln(fou);
Writeln(fou,'k2024 = ',k2024:13:5,'      k2420 = ',k2420:13:5);
Writeln(fou,'k2125 = ',k2125:13:5,'      k2521 = ',k2521:13:5);
Writeln(fou,'k2226 = ',k2226:13:5,'      k2622 = ',k2622:13:5);
Writeln(fou,'k2327 = ',k2327:13:5,'      k2723 = ',k2723:13:5);
Writeln(fou);
Writeln(fou);
Writeln(fou,'k014 = ',k014:13:5,'      k140 = ',k140:13:5);
Writeln(fou,'k115 = ',k115:13:5,'      k151 = ',k151:13:5);
Writeln(fou,'k216 = ',k216:13:5,'      k162 = ',k162:13:5);
Writeln(fou,'k317 = ',k317:13:5,'      k173 = ',k173:13:5);
Writeln(fou,'k418 = ',k418:13:5,'      k184 = ',k184:13:5);
Writeln(fou,'k519 = ',k519:13:5,'      k195 = ',k195:13:5);
Writeln(fou,'k620 = ',k620:13:5,'      k206 = ',k206:13:5);
Writeln(fou,'k721 = ',k721:13:5,'      k217 = ',k217:13:5);
Writeln(fou,'k822 = ',k822:13:5,'      k228 = ',k228:13:5);
Writeln(fou,'k923 = ',k923:13:5,'      k239 = ',k239:13:5);
Writeln(fou);
Writeln(fou,'k1024 = ',k1024:13:5,'      k2410 = ',k2410:13:5);
Writeln(fou,'k1125 = ',k1125:13:5,'      k2511 = ',k2511:13:5);
Writeln(fou,'k1226 = ',k1226:13:5,'      k2612 = ',k2612:13:5);
Writeln(fou,'k1327 = ',k1327:13:5,'      k2713 = ',k2713:13:5);
Writeln(fou);
CloseFile(fou);
end;

procedure TForm1.Button1Click(Sender: TObject);
var fop:TextFile;
par:real;
ps,psn:string;
pj,pk,pm:integer;
begin
ForceCurrentDirectory:=True;
OpenDialog1.DefaultExt:='tpt';
OpenDialog1.Filter:='Parameter files(*.tpt)|*.tpt';
OpenDialog1.FilterIndex := 1;
if OpenDialog1.Execute then
begin
ps:=OpenDialog1.FileName;
pk:=Length(ps);
pm:=1;
for pj:=1 to pk do if ps[pj]='\' then pm:=pj;
psn:=ps;
Delete(psn,1,pm);
AssignFile(fop,ps);
Reset(fop);
readln(fop,par);
Edit1.Text:=FloatToStrF(par,ffFixed,9,1);
readln(fop,par);
Edit2.Text:=FloatToStrF(par,ffFixed,9,1);
readln(fop,par);
Edit3.Text:=FloatToStrF(par,ffFixed,9,1);
readln(fop,par);
Edit4.Text:=FloatToStrF(par,ffFixed,9,1);
readln(fop,par);
Edit5.Text:=FloatToStrF(par,ffFixed,9,1);
readln(fop,par);
Edit6.Text:=FloatToStrF(par,ffFixed,9,5);
readln(fop,par);
Edit7.Text:=FloatToStrF(par,ffFixed,9,1);
readln(fop,par);
Edit8.Text:=FloatToStrF(par,ffFixed,9,5);
readln(fop,par);
Edit19.Text:=FloatToStrF(par,ffFixed,9,2);
readln(fop,par);
Edit20.Text:=FloatToStrF(par,ffFixed,9,4);
readln(fop,par);
Edit21.Text:=FloatToStrF(par,ffFixed,9,2);
readln(fop,par);
Edit22.Text:=FloatToStrF(par,ffFixed,9,4);
readln(fop,par);
Edit26.Text:=FloatToStrF(par,ffFixed,9,2);
readln(fop,par);
Edit27.Text:=FloatToStrF(par,ffFixed,9,4);
readln(fop,par);
Edit28.Text:=FloatToStrF(par,ffFixed,9,2);
readln(fop,par);
Edit29.Text:=FloatToStrF(par,ffFixed,9,4);
readln(fop,par);
Edit37.Text:=FloatToStrF(par,ffFixed,9,2);
readln(fop,par);
Edit38.Text:=FloatToStrF(par,ffFixed,9,4);
readln(fop,par);
Edit39.Text:=FloatToStrF(par,ffFixed,9,2);
readln(fop,par);
Edit40.Text:=FloatToStrF(par,ffFixed,9,4);
readln(fop,par);
Edit51.Text:=FloatToStrF(par,ffFixed,9,2);
readln(fop,par);
Edit12.Text:=FloatToStrF(par,ffFixed,9,4);
readln(fop,par);
Edit17.Text:=FloatToStrF(par,ffFixed,9,2);
readln(fop,par);
Edit18.Text:=FloatToStrF(par,ffFixed,9,4);
readln(fop,par);
Edit53.Text:=FloatToStrF(par,ffFixed,9,5);
readln(fop,par);
Edit54.Text:=FloatToStrF(par,ffFixed,9,5);
readln(fop,par);
Edit55.Text:=FloatToStrF(par,ffFixed,9,5);
readln(fop,par);
Edit56.Text:=FloatToStrF(par,ffFixed,9,5);
readln(fop,par);
Edit57.Text:=FloatToStrF(par,ffFixed,9,5);
readln(fop,par);
Edit58.Text:=FloatToStrF(par,ffFixed,9,5);
readln(fop,par);
Edit59.Text:=FloatToStrF(par,ffFixed,9,5);
readln(fop,par);
Edit60.Text:=FloatToStrF(par,ffFixed,9,5);
readln(fop,par);
Edit61.Text:=FloatToStrF(par,ffFixed,9,5);
readln(fop,par);
Edit62.Text:=FloatToStrF(par,ffFixed,9,5);
readln(fop,par);
Edit63.Text:=FloatToStrF(par,ffFixed,9,5);
readln(fop,par);
Edit64.Text:=FloatToStrF(par,ffFixed,9,5);
readln(fop,par);
Edit65.Text:=FloatToStrF(par,ffFixed,9,5);
readln(fop,par);
Edit66.Text:=FloatToStrF(par,ffFixed,9,5);
readln(fop,par);
Edit67.Text:=FloatToStrF(par,ffFixed,9,5);
readln(fop,par);
Edit68.Text:=FloatToStrF(par,ffFixed,9,5);
readln(fop,par);
Edit69.Text:=FloatToStrF(par,ffFixed,9,5);
readln(fop,par);
Edit70.Text:=FloatToStrF(par,ffFixed,9,5);
readln(fop,par);
Edit71.Text:=FloatToStrF(par,ffFixed,9,5);
readln(fop,par);
Edit72.Text:=FloatToStrF(par,ffFixed,9,5);
readln(fop,par);
Edit94.Text:=FloatToStrF(par,ffFixed,9,5);
readln(fop,par);
Edit95.Text:=FloatToStrF(par,ffFixed,9,5);
readln(fop,par);
Edit96.Text:=FloatToStrF(par,ffFixed,9,5);
readln(fop,par);
Edit97.Text:=FloatToStrF(par,ffFixed,9,5);
readln(fop,par);
Edit98.Text:=FloatToStrF(par,ffFixed,9,5);
readln(fop,par);
Edit99.Text:=FloatToStrF(par,ffFixed,9,5);
readln(fop,par);
Edit100.Text:=FloatToStrF(par,ffFixed,9,5);
readln(fop,par);
Edit101.Text:=FloatToStrF(par,ffFixed,9,5);
readln(fop,par);
Edit73.Text:=FloatToStrF(par,ffFixed,9,7);
readln(fop,par);
Edit74.Text:=FloatToStrF(par,ffFixed,7,2);
readln(fop,par);
Edit75.Text:=FloatToStrF(par,ffFixed,7,2);
readln(fop,par);
Edit76.Text:=FloatToStrF(par,ffFixed,7,2);
readln(fop,par);
Edit9.Text:=FloatToStrF(par,ffFixed,9,7);
readln(fop,par);
Edit10.Text:=FloatToStrF(par,ffFixed,7,2);
readln(fop,par);
Edit31.Text:=FloatToStrF(par,ffFixed,7,2);
readln(fop,par);
Edit42.Text:=FloatToStrF(par,ffFixed,7,2);
readln(fop,par);
Edit46.Text:=FloatToStrF(par,ffFixed,8,4);
readln(fop,par);
Edit47.Text:=FloatToStrF(par,ffFixed,8,4);
Edit48.Text:=psn;
CloseFile(fop);
end
else
begin
Edit1.Text:='2300.0';
Edit2.Text:='12.0';
Edit3.Text:='40000.0';
Edit4.Text:='12.0';
Edit5.Text:='10.0';
Edit6.Text:='0.2';
Edit7.Text:='10.0';
Edit8.Text:='0.02';
Edit19.Text:='140.0';
Edit20.Text:='1.0';
Edit21.Text:='7.0';
Edit22.Text:='0.1';
Edit26.Text:='105.0';
Edit27.Text:='1.0';
Edit28.Text:='8.7';
Edit29.Text:='0.1';
Edit37.Text:='70.0';
Edit38.Text:='1.0';
Edit39.Text:='10.4';
Edit40.Text:='0.1';
Edit51.Text:='35.0';
Edit12.Text:='1.0';
Edit17.Text:='12.0';
Edit18.Text:='0.1';
Edit53.Text:='0.00001';
Edit54.Text:='0.25';
Edit55.Text:='1.0';
Edit56.Text:='0.00001';
Edit57.Text:='12.0';
Edit58.Text:='0.00001';
Edit59.Text:='0.8';
Edit60.Text:='0.00001';
Edit61.Text:='9.0';
Edit62.Text:='0.00001';
Edit63.Text:='0.5';
Edit64.Text:='0.00001';
Edit65.Text:='6.0';
Edit66.Text:='0.00001';
Edit67.Text:='0.4';
Edit68.Text:='0.00001';
Edit69.Text:='0.67';
Edit70.Text:='0.00001';
Edit71.Text:='0.53';
Edit72.Text:='0.00001';

Edit94.Text:='1.0';
Edit95.Text:='0.12';
Edit96.Text:='0.7';
Edit97.Text:='0.12';
Edit98.Text:='0.4';
Edit99.Text:='0.12';
Edit100.Text:='0.1';
Edit101.Text:='0.12';

Edit73.Text:='0.01';
Edit74.Text:='90.0';
Edit75.Text:='2.0';
Edit76.Text:='0.0';
Edit9.Text:='0.000002';
Edit10.Text:='0.0';
Edit31.Text:='90.0';
Edit42.Text:='0.0';
Edit46.Text:='0.04';
Edit47.Text:='0.04';
Edit48.Text:='control.tpt';
end;
Edit11.Text:='0.001';
Edit13.Text:='1.0';
rish:=1.0;
Edit14.Text:='0.0';
Edit15.Text:='0.0';
Edit16.Text:='0.0';
Edit23.Text:='0.0';
Edit30.Text:='0.0';
Edit33.Text:='0.0';
Edit34.Text:='0.0';
Edit35.Text:='0.0';
Edit36.Text:='0.0';
Edit41.Text:='0.0';
Edit43.Text:='0.0';
Edit45.Text:='0.0';
Edit52.Text:='0.0';
Edit77.Text:='0.0';
Edit78.Text:='0.0';
Edit79.Text:='0.0';
Edit80.Text:='0.0';
Edit81.Text:='0.0';
Edit82.Text:='0.0';
Edit83.Text:='0.0';
Edit84.Text:='0.0';
Edit85.Text:='0.0';
Edit86.Text:='0.0';
Edit87.Text:='0.0';
Edit88.Text:='0.0';
Edit89.Text:='0.0';
Edit90.Text:='0.0';
Edit44.Text:='???.?';
Edit24.Text:='???.?';
Edit49.Text:='1.0';
Edit50.Text:='???.?';
Edit25.Text:='???.?';
Edit32.Text:='???.?';
Edit91.Text:='1.0';
Edit92.Text:='0.0';
Edit93.Text:='0.0';
end;

procedure TForm1.Button4Click(Sender: TObject);
begin
Series3.Clear;
Series2.Clear;
Series1.Clear;
Series4.Clear;
Series5.Clear;
end;

procedure TForm1.Button2Click(Sender: TObject);
label 1,2,51,79,99;
const
   Nmax=10000000;
   uround:real=1.73e-18;
   c2:real=1./18. ;
   c3:real=1./12. ;
   c4:real=1./8.  ;
   c5:real=5./16. ;
   c6:real=3./8.   ;
   c7:real=59./400.;
   c8:real=93./200.;
   c9:real=5490023248./9719169821.;
   c10:real=13./20. ;
   c11:real=1201146811./1299019798.;
   c12:real=1.      ;
   c13:real=1.      ;
   a21:real=1./18.  ;
   a31:real=1./48.   ;
   a32:real=1./16.   ;
   a41:real=1./32.   ;
   a43:real=3./32.    ;
   a51:real=5./16.    ;
   a53:real= -75./64. ;
   a54:real= 75./64.    ;
   a61:real=3./80.    ;
   a64:real=3./16.    ;
   a65:real=3./20.     ;
   a71:real=29443841./614563906.   ;
   a74:real=77736538./692538347.   ;
   a75:real= -28693883./1125.e6    ;
   a76:real=23124283./18.e8        ;
   a81:real=16016141./946692911.    ;
   a84:real=61564180./158732637.     ;
   a85:real=22789713./633445777.     ;
   a86:real=545815736./2771057229.   ;
   a87:real= -180193667./1043307555. ;
   a91:real=39632708./573591083.      ;
   a94:real= -433636366./683701615.   ;
   a95:real= -421739975./2616292301.  ;
   a96:real=100302831./723423059.     ;
   a97:real=790204164./839813087.     ;
   a98:real=800635310./3783071287.     ;
   a101:real=246121993./1340847787.     ;
   a104:real= -37695042795./15268766246. ;
   a105:real= -309121744./1061227803.    ;
   a106:real= -12992083./490766935.      ;
   a107:real=6005943493./2108947869.     ;
   a108:real=393006217./1396673457.      ;
   a109:real=123872331./1001029789.      ;
   a111:real= -1028468189./846180014.    ;
   a114:real=8478235783./508512852.      ;
   a115:real=1311729495./1432422823.      ;
   a116:real= -10304129995./1701304382.    ;
   a117:real= -48777925059./3047939560.   ;
   a118:real=15336726248./1032824649.       ;
   a119:real= -45442868181./3398467696.     ;
   a1110:real=3065993473./597172653.        ;
   a121:real=185892177./718116043.          ;
   a124:real= -3185094517./667107341.       ;
   a125:real= -477755414./1098053517.       ;
   a126:real= -703635378./230739211.        ;
   a127:real=5731566787./1027545527.        ;
   a128:real=5232866602./850066563.         ;
   a129:real= -4093664535./808688257.       ;
   a1210:real=3962137247./1805957418.       ;
   a1211:real=65686358./487910083.          ;
   a131:real=403863854./491063109.           ;
   a134:real= -5068492393./434740067.         ;
   a135:real= -411421997./543043805.           ;
   a136:real=652783627./914296604.              ;
   a137:real=11173962825./925320556.      ;
   a138:real= -13158990841./6184727034.   ;
   a139:real=3936647629./1978049680.      ;
   a1310:real= -160528059./685178525.     ;
   a1311:real=248638103./1413531060.      ;
   b1:real=14005451./335480064.           ;
   b6:real= -59238493./1068277825.        ;
   b7:real=181606767./758867731.          ;
   b8:real=561292985./797845732.          ;
   b9:real= -1041891430./1371343529.      ;
   b10:real=760417239./1151165299.        ;
   b11:real=118820643./751138087.         ;
   b12:real= -528747749./2220607170.      ;
   b13:real=1./4.                         ;
   bh1:real=13451932./455176623.     ;
   bh6:real= -808719846./976000145.  ;
   bh7:real=1757004468./5645159321.  ;
   bh8:real=656045339./265891186.      ;
   bh9:real= -3867574721./1518517206.  ;
   bh10:real=465885868./322736535.     ;
   bh11:real=53011238./667516719.      ;
   bh12:real=2./45. ;
var j,jj,ii,im,i,nstep,imi:integer;
dp,x,xend,xold,eps,hmax,h,su,r1,ddp,
t1,t9,t2,t8,a2,a8,r,mam,
posneg,y11s,y12s,hnew,xph,err,denom,fac,
tt,tta,ttb,xx,a5,ami,amo,t5,aka,bka:real;
k1,k2,k3,k4,k5,k6,k7,y1,y:array[1..27] of real;
reject,quit:boolean;
{fff:TextFile;   }
begin
{ Assignfile(fff,'sig.dat');
  Rewrite(fff);}
mam:=StrToFloat(Edit49.Text);
if mam<=0.0 then mam:=1.0;
agon:=StrToFloat(Edit73.Text);
tt0:=StrToFloat(Edit74.Text);
tt1:=StrToFloat(Edit75.Text);
tt2:=StrToFloat(Edit76.Text);
bgon:=StrToFloat(Edit9.Text);
tb0:=StrToFloat(Edit10.Text);
tb1:=StrToFloat(Edit31.Text);
tb2:=StrToFloat(Edit42.Text);
tta:=tt0+tt1+tt2;
ttb:=tb0+tb1+tb2;
tt:=Max(tta,ttb);
tauf:=StrToFloat(Edit46.Text);
taus:=StrToFloat(Edit47.Text);
dp:=StrToFloat(Edit11.Text);
if tt<2.0 then dp:=dp*tt/2.0;

  if CheckBox1.Checked=false then
  begin
if agon>=bgon then
begin
  if tt0<=2.0 then jj:=Trunc(tt0/(10.0*dp))+1 else
  jj:=Trunc((2.0/(10.0*dp))+((tt0-2.0)/(100.0*dp)))+1;
  if tt1>0.0 then if tt1<=6.0 then
  jj:=jj+Trunc(tt1/dp)+1 else
  jj:=jj+Trunc((6.0/dp)+((tt1-6.0)/(10.0*dp)))+1;
  if tt2>0.0 then
  jj:=jj+Trunc(tt2/(100.0*dp))+1;
  if jj<3 then jj:=3;
end
else
begin
  if tb0<=2.0 then jj:=Trunc(tb0/(10.0*dp))+1 else
  jj:=Trunc((2.0/(10.0*dp))+((tb0-2.0)/(100.0*dp)))+1;
  if tb1>0.0 then if tb1<=6.0 then
  jj:=jj+Trunc(tb1/dp)+1 else
  jj:=jj+Trunc((6.0/dp)+((tb1-6.0)/(10.0*dp)))+1;
  if tb2>0.0 then
  jj:=jj+Trunc(tb2/(100.0*dp))+1;
  if jj<3 then jj:=3;
end
  end
  else
  jj:=Trunc(tt/(500.0*dp))+1;

  ProgressBar1.Min:=0;
  ProgressBar1.Max:=jj;
  ProgressBar1.Position:=0;

k36:=StrToFloat(Edit1.Text);
k47:=k36; k58:=k36; k13:=2*k36; k24:=k13; k01:=3*k36;
k10:=StrToFloat(Edit2.Text);
k31:=2*k10; k42:=k31; k63:=3*k10; k74:=k63; k85:=k63;
k37:=StrToFloat(Edit3.Text);
k48:=k37; k59:=k37; k14:=2*k37; k25:=k14; k02:=3*k37;
k20:=StrToFloat(Edit4.Text);
k52:=2*k20; k41:=k52; k73:=3*k20; k84:=k73; k95:=k73;

k1720:=StrToFloat(Edit5.Text);
k1821:=k1720; k1922:=k1720; k1517:=2*k1720; k1618:=k1517; k1415:=3*k1720;
k1514:=StrToFloat(Edit6.Text);
k1715:=2*k1514; k1816:=k1715; k2017:=3*k1514; k2118:=k2017; k2219:=k2017;
k1721:=StrToFloat(Edit7.Text);
k1822:=k1721; k1923:=k1721; k1518:=2*k1721; k1619:=k1518; k1416:=3*k1721;
k1614:=StrToFloat(Edit8.Text);
k1916:=2*k1614; k1815:=k1916; k2117:=3*k1614; k2218:=k2117; k2319:=k2117;

k2024:=StrToFloat(Edit95.Text);
k2420:=StrToFloat(Edit94.Text);
k2125:=StrToFloat(Edit97.Text);
k2521:=StrToFloat(Edit96.Text);
k2226:=StrToFloat(Edit99.Text);
k2622:=StrToFloat(Edit98.Text);
k2327:=StrToFloat(Edit101.Text);
k2723:=StrToFloat(Edit100.Text);

k610:=StrToFloat(Edit19.Text);
k106:=StrToFloat(Edit20.Text);
k1024:=StrToFloat(Edit21.Text);
k2410:=StrToFloat(Edit22.Text);
k711:=StrToFloat(Edit26.Text);
k117:=StrToFloat(Edit27.Text);
k1125:=StrToFloat(Edit28.Text);
k2511:=StrToFloat(Edit29.Text);

k812:=StrToFloat(Edit37.Text);
k128:=StrToFloat(Edit38.Text);
k1226:=StrToFloat(Edit39.Text);
k2612:=StrToFloat(Edit40.Text);
k913:=StrToFloat(Edit51.Text);
k139:=StrToFloat(Edit12.Text);
k1327:=StrToFloat(Edit17.Text);
k2713:=StrToFloat(Edit18.Text);
k014:=StrToFloat(Edit53.Text);
k140:=StrToFloat(Edit54.Text);

k115:=StrToFloat(Edit55.Text);
k151:=StrToFloat(Edit56.Text);
k317:=StrToFloat(Edit57.Text);
k173:=StrToFloat(Edit58.Text);
k620:=StrToFloat(Edit59.Text);
k206:=StrToFloat(Edit60.Text);
k418:=StrToFloat(Edit61.Text);
k184:=StrToFloat(Edit62.Text);

k216:=StrToFloat(Edit63.Text);
k162:=StrToFloat(Edit64.Text);
k519:=StrToFloat(Edit65.Text);
k195:=StrToFloat(Edit66.Text);
k923:=StrToFloat(Edit67.Text);
k239:=StrToFloat(Edit68.Text);
k721:=StrToFloat(Edit69.Text);
k217:=StrToFloat(Edit70.Text);
k822:=StrToFloat(Edit71.Text);
k228:=StrToFloat(Edit72.Text);

rr[0]:=StrToFloat(Edit13.Text);
ar[0]:=StrToFloat(Edit14.Text);
br[0]:=StrToFloat(Edit15.Text);
a2r[0]:=StrToFloat(Edit16.Text);
abr[0]:=StrToFloat(Edit23.Text);
b2r[0]:=StrToFloat(Edit30.Text);
a3r[0]:=StrToFloat(Edit33.Text);
a2br[0]:=StrToFloat(Edit34.Text);
ab2r[0]:=StrToFloat(Edit35.Text);
b3r[0]:=StrToFloat(Edit36.Text);
a3ro[0]:=StrToFloat(Edit41.Text);
a2bro[0]:=StrToFloat(Edit43.Text);
ab2ro[0]:=StrToFloat(Edit45.Text);
b3ro[0]:=StrToFloat(Edit52.Text);
d[0]:=StrToFloat(Edit77.Text);
ad[0]:=StrToFloat(Edit78.Text);
bd[0]:=StrToFloat(Edit79.Text);
a2d[0]:=StrToFloat(Edit80.Text);
abd[0]:=StrToFloat(Edit81.Text);
b2d[0]:=StrToFloat(Edit82.Text);
a3d[0]:=StrToFloat(Edit83.Text);
a2bd[0]:=StrToFloat(Edit84.Text);
ab2d[0]:=StrToFloat(Edit85.Text);
b3d[0]:=StrToFloat(Edit86.Text);
a3do[0]:=StrToFloat(Edit87.Text);
a2bdo[0]:=StrToFloat(Edit88.Text);
ab2do[0]:=StrToFloat(Edit89.Text);
b3do[0]:=StrToFloat(Edit90.Text);
tx[0]:=0.0;
sumr[0]:=rr[0]+ar[0]+br[0]+a2r[0]+abr[0]+b2r[0]+a3r[0]+a2br[0]+ab2r[0]+b3r[0];
sumo[0]:=a3ro[0]+a2bro[0]+ab2ro[0]+b3ro[0];
sumd[0]:=d[0]+ad[0]+bd[0]+a2d[0]+abd[0]+b2d[0]+a3d[0]+a2bd[0]+ab2d[0]+b3d[0]
+a3do[0]+a2bdo[0]+ab2do[0]+b3do[0];
su:=sumr[0]+sumo[0]+sumd[0];
if (su>1.001) OR (su<0.999) then
begin
rr[0]:=rr[0]/su; ar[0]:=ar[0]/su; br[0]:=br[0]/su;
a2r[0]:=a2r[0]/su; abr[0]:=abr[0]/su; b2r[0]:=b2r[0]/su;
a3r[0]:=a3r[0]/su; a2br[0]:=a2br[0]/su; ab2r[0]:=ab2r[0]/su;
b3r[0]:=b3r[0]/su; a3ro[0]:=a3ro[0]/su; a2bro[0]:=a2bro[0]/su;
ab2ro[0]:=ab2ro[0]/su; b3ro[0]:=b3ro[0]/su;
d[0]:=d[0]/su; ad[0]:=ad[0]/su; bd[0]:=bd[0]/su;
a2d[0]:=a2d[0]/su; abd[0]:=abd[0]/su; b2d[0]:=b2d[0]/su;
a3d[0]:=a3d[0]/su; a2bd[0]:=a2bd[0]/su; ab2d[0]:=ab2d[0]/su;
b3d[0]:=b3d[0]/su; a3do[0]:=a3do[0]/su; a2bdo[0]:=a2bdo[0]/su;
ab2do[0]:=ab2do[0]/su; b3do[0]:=b3do[0]/su;
sumr[0]:=rr[0]+ar[0]+br[0]+a2r[0]+abr[0]+b2r[0]+a3r[0]+a2br[0]+ab2r[0]+b3r[0];
sumo[0]:=a3ro[0]+a2bro[0]+ab2ro[0]+b3ro[0];
sumd[0]:=d[0]+ad[0]+bd[0]+a2d[0]+abd[0]+b2d[0]+a3d[0]+a2bd[0]+ab2d[0]+b3d[0]
+a3do[0]+a2bdo[0]+ab2do[0]+b3do[0];
end;
Edit91.Text:=FloatToStrF(sumr[jj],ffFixed,7,4);
Edit92.Text:=FloatToStrF(sumo[jj],ffFixed,7,4);
Edit93.Text:=FloatToStrF(sumd[jj],ffFixed,7,4);

xend:=tt; eps:=0.05; h:=0.0001; hmax:=0.5*dp;
y[1]:=ar[0]; y[2]:=br[0]; y[3]:=a2r[0]; y[4]:=abr[0]; y[5]:=b2r[0];
y[6]:=a3r[0]; y[7]:=a2br[0]; y[8]:=ab2r[0]; y[9]:=b3r[0];
y[10]:=a3ro[0]; y[11]:=a2bro[0]; y[12]:=ab2ro[0]; y[13]:=b3ro[0];
y[14]:=d[0];
y[15]:=ad[0]; y[16]:=bd[0]; y[17]:=a2d[0]; y[18]:=abd[0]; y[19]:=b2d[0];
y[20]:=a3d[0]; y[21]:=a2bd[0]; y[22]:=ab2d[0]; y[23]:=b3d[0];
y[24]:=a3do[0]; y[25]:=a2bdo[0]; y[26]:=ab2do[0]; y[27]:=b3do[0];

jj:=0; x:=0.0; xold:=x;
        {begin of Runge-Kutta calc}

  quit:=false;
  posneg:=Sign(xend-x);
  hmax:=abs(hmax);
  h:=min(max(1.0e-10,abs(h)),hmax);
  h:=h*posneg;
  reject:=false; nstep:=0;
  {naccpt:=0; nreject:=0; nfcn:=0;}
  Series3.AddXY(x,sumo[0],'',clBlue);
  Series1.AddXY(x,sumr[0],'',clRed);
  Series4.AddXY(x,sumd[0],'',clLime);
  Series2.AddXY(x,0.0,'',clYellow);
  Series5.AddXY(x,0.0,'',clAqua);
{  writeln(fff,j:5,jj:8,y[3]:8:5);}
  { Basic Integration Step }
1:
  if {quit or} (nstep>nmax) or ((x+0.03*h)=x) then goto 79;
  if (x-xend)*posneg>0 then goto 79;        { !!!!!!!!!!! }
  if (x+h-xend)*posneg>0 then h:=xend-x;
  if CheckBox1.Checked=false then
  if x>tt0 then ag:=vspl(x,tt0,tt1,agon) else ag:=0.0
  else ag:=0.0;
  r:=1.0-y[1]-y[2]-y[3]-y[4]-y[5]-y[6]-y[7]-y[8]-y[9]
  -y[10]-y[11]-y[12]-y[13]-y[14]-y[15]-y[16]-y[17]-y[18]
  -y[19]-y[20]-y[21]-y[22]-y[23]-y[24]-y[25]-y[26]-y[27];
  k1[1]:=so1(r,y[1],y[3],y[4],y[15]);
  k1[2]:=so2(r,y[2],y[4],y[5],y[16]);
  k1[3]:=so3(y[1],y[3],y[6],y[7],y[17]);
  k1[4]:=so4(y[1],y[2],y[4],y[7],y[8],y[18]);
  k1[5]:=so5(y[2],y[5],y[8],y[9],y[19]);
  k1[6]:=so6(y[3],y[6],y[10],y[20]);
  k1[7]:=so7(y[3],y[4],y[7],y[11],y[21]);
  k1[8]:=so8(y[4],y[5],y[8],y[12],y[22]);
  k1[9]:=so9(y[5],y[9],y[13],y[23]);
  k1[10]:=so10(y[6],y[10],y[24]);
  k1[11]:=so11(y[7],y[11],y[25]);
  k1[12]:=so12(y[8],y[12],y[26]);
  k1[13]:=so13(y[9],y[13],y[27]);
  k1[14]:=so14(r,y[14],y[15],y[16]);
  k1[15]:=so15(y[14],y[15],y[17],y[18],y[1]);
  k1[16]:=so16(y[14],y[16],y[18],y[19],y[2]);
  k1[17]:=so17(y[15],y[17],y[20],y[21],y[3]);
  k1[18]:=so18(y[15],y[16],y[18],y[21],y[22],y[4]);
  k1[19]:=so19(y[16],y[19],y[22],y[23],y[5]);
  k1[20]:=so20(y[17],y[20],y[24],y[6]);
  k1[21]:=so21(y[17],y[18],y[21],y[25],y[7]);
  k1[22]:=so22(y[18],y[19],y[22],y[26],y[8]);
  k1[23]:=so23(y[19],y[23],y[27],y[9]);
  k1[24]:=so24(y[20],y[24],y[10]);
  k1[25]:=so25(y[21],y[25],y[11]);
  k1[26]:=so26(y[22],y[26],y[12]);
  k1[27]:=so27(y[23],y[27],y[13]);
2:
{  inc(nstep);}
{ the First 9 stages }
  for i:=1 to 27 do
  y1[i]:=y[i]+h*a21*k1[i];
  r:=1.0-y1[1]-y1[2]-y1[3]-y1[4]-y1[5]-y1[6]-y1[7]-y1[8]-y1[9]
  -y1[10]-y1[11]-y1[12]-y1[13]-y1[14]-y1[15]-y1[16]-y1[17]-y1[18]
  -y1[19]-y1[20]-y1[21]-y1[22]-y1[23]-y1[24]-y1[25]-y1[26]-y1[27];
  k2[1]:=so1(r,y1[1],y1[3],y1[4],y1[15]);
  k2[2]:=so2(r,y1[2],y1[4],y1[5],y1[16]);
  k2[3]:=so3(y1[1],y1[3],y1[6],y1[7],y1[17]);
  k2[4]:=so4(y1[1],y1[2],y1[4],y1[7],y1[8],y1[18]);
  k2[5]:=so5(y1[2],y1[5],y1[8],y1[9],y1[19]);
  k2[6]:=so6(y1[3],y1[6],y1[10],y1[20]);
  k2[7]:=so7(y1[3],y1[4],y1[7],y1[11],y1[21]);
  k2[8]:=so8(y1[4],y1[5],y1[8],y1[12],y1[22]);
  k2[9]:=so9(y1[5],y1[9],y1[13],y1[23]);
  k2[10]:=so10(y1[6],y1[10],y1[24]);
  k2[11]:=so11(y1[7],y1[11],y1[25]);
  k2[12]:=so12(y1[8],y1[12],y1[26]);
  k2[13]:=so13(y1[9],y1[13],y1[27]);
  k2[14]:=so14(r,y1[14],y1[15],y1[16]);
  k2[15]:=so15(y1[14],y1[15],y1[17],y1[18],y1[1]);
  k2[16]:=so16(y1[14],y1[16],y1[18],y1[19],y1[2]);
  k2[17]:=so17(y1[15],y1[17],y1[20],y1[21],y1[3]);
  k2[18]:=so18(y1[15],y1[16],y1[18],y1[21],y1[22],y1[4]);
  k2[19]:=so19(y1[16],y1[19],y1[22],y1[23],y1[5]);
  k2[20]:=so20(y1[17],y1[20],y1[24],y1[6]);
  k2[21]:=so21(y1[17],y1[18],y1[21],y1[25],y1[7]);
  k2[22]:=so22(y1[18],y1[19],y1[22],y1[26],y1[8]);
  k2[23]:=so23(y1[19],y1[23],y1[27],y1[9]);
  k2[24]:=so24(y1[20],y1[24],y1[10]);
  k2[25]:=so25(y1[21],y1[25],y1[11]);
  k2[26]:=so26(y1[22],y1[26],y1[12]);
  k2[27]:=so27(y1[23],y1[27],y1[13]);
  for i:=1 to 27 do
  y1[i]:=y[i]+h*(a31*k1[i]+a32*k2[i]);
  r1:=1.0-y1[1]-y1[2]-y1[3]-y1[4]-y1[5]-y1[6]-y1[7]-y1[8]-y1[9]
  -y1[10]-y1[11]-y1[12]-y1[13]-y1[14]-y1[15]-y1[16]-y1[17]-y1[18]
  -y1[19]-y1[20]-y1[21]-y1[22]-y1[23]-y1[24]-y1[25]-y1[26]-y1[27];
  xx:=x+c3*h;
  if CheckBox1.Checked=false then
  if xx>tt0 then ag:=vspl(xx,tt0,tt1,agon) else ag:=0.0
  else ag:=0.0;
  if CheckBox1.Checked=false then
  if xx>tb0 then bg:=vspl(xx,tb0,tb1,bgon) else bg:=0.0
  else bg:=0.0;
  k3[1]:=so1(r1,y1[1],y1[3],y1[4],y1[15]);
  k3[2]:=so2(r1,y1[2],y1[4],y1[5],y1[16]);
  k3[3]:=so3(y1[1],y1[3],y1[6],y1[7],y1[17]);
  k3[4]:=so4(y1[1],y1[2],y1[4],y1[7],y1[8],y1[18]);
  k3[5]:=so5(y1[2],y1[5],y1[8],y1[9],y1[19]);
  k3[6]:=so6(y1[3],y1[6],y1[10],y1[20]);
  k3[7]:=so7(y1[3],y1[4],y1[7],y1[11],y1[21]);
  k3[8]:=so8(y1[4],y1[5],y1[8],y1[12],y1[22]);
  k3[9]:=so9(y1[5],y1[9],y1[13],y1[23]);
  k3[10]:=so10(y1[6],y1[10],y1[24]);
  k3[11]:=so11(y1[7],y1[11],y1[25]);
  k3[12]:=so12(y1[8],y1[12],y1[26]);
  k3[13]:=so13(y1[9],y1[13],y1[27]);
  k3[14]:=so14(r1,y1[14],y1[15],y1[16]);
  k3[15]:=so15(y1[14],y1[15],y1[17],y1[18],y1[1]);
  k3[16]:=so16(y1[14],y1[16],y1[18],y1[19],y1[2]);
  k3[17]:=so17(y1[15],y1[17],y1[20],y1[21],y1[3]);
  k3[18]:=so18(y1[15],y1[16],y1[18],y1[21],y1[22],y1[4]);
  k3[19]:=so19(y1[16],y1[19],y1[22],y1[23],y1[5]);
  k3[20]:=so20(y1[17],y1[20],y1[24],y1[6]);
  k3[21]:=so21(y1[17],y1[18],y1[21],y1[25],y1[7]);
  k3[22]:=so22(y1[18],y1[19],y1[22],y1[26],y1[8]);
  k3[23]:=so23(y1[19],y1[23],y1[27],y1[9]);
  k3[24]:=so24(y1[20],y1[24],y1[10]);
  k3[25]:=so25(y1[21],y1[25],y1[11]);
  k3[26]:=so26(y1[22],y1[26],y1[12]);
  k3[27]:=so27(y1[23],y1[27],y1[13]);
  for i:=1 to 27 do
  y1[i]:=y[i]+h*(a41*k1[i]+a43*k3[i]);
  r1:=1.0-y1[1]-y1[2]-y1[3]-y1[4]-y1[5]-y1[6]-y1[7]-y1[8]-y1[9]
  -y1[10]-y1[11]-y1[12]-y1[13]-y1[14]-y1[15]-y1[16]-y1[17]-y1[18]
  -y1[19]-y1[20]-y1[21]-y1[22]-y1[23]-y1[24]-y1[25]-y1[26]-y1[27];
  xx:=x+c4*h;
  if CheckBox1.Checked=false then
  if xx>tt0 then ag:=vspl(xx,tt0,tt1,agon) else ag:=0.0
  else ag:=0.0;
  if CheckBox1.Checked=false then
  if xx>tb0 then bg:=vspl(xx,tb0,tb1,bgon) else bg:=0.0
  else bg:=0.0;
  k4[1]:=so1(r1,y1[1],y1[3],y1[4],y1[15]);
  k4[2]:=so2(r1,y1[2],y1[4],y1[5],y1[16]);
  k4[3]:=so3(y1[1],y1[3],y1[6],y1[7],y1[17]);
  k4[4]:=so4(y1[1],y1[2],y1[4],y1[7],y1[8],y1[18]);
  k4[5]:=so5(y1[2],y1[5],y1[8],y1[9],y1[19]);
  k4[6]:=so6(y1[3],y1[6],y1[10],y1[20]);
  k4[7]:=so7(y1[3],y1[4],y1[7],y1[11],y1[21]);
  k4[8]:=so8(y1[4],y1[5],y1[8],y1[12],y1[22]);
  k4[9]:=so9(y1[5],y1[9],y1[13],y1[23]);
  k4[10]:=so10(y1[6],y1[10],y1[24]);
  k4[11]:=so11(y1[7],y1[11],y1[25]);
  k4[12]:=so12(y1[8],y1[12],y1[26]);
  k4[13]:=so13(y1[9],y1[13],y1[27]);
  k4[14]:=so14(r1,y1[14],y1[15],y1[16]);
  k4[15]:=so15(y1[14],y1[15],y1[17],y1[18],y1[1]);
  k4[16]:=so16(y1[14],y1[16],y1[18],y1[19],y1[2]);
  k4[17]:=so17(y1[15],y1[17],y1[20],y1[21],y1[3]);
  k4[18]:=so18(y1[15],y1[16],y1[18],y1[21],y1[22],y1[4]);
  k4[19]:=so19(y1[16],y1[19],y1[22],y1[23],y1[5]);
  k4[20]:=so20(y1[17],y1[20],y1[24],y1[6]);
  k4[21]:=so21(y1[17],y1[18],y1[21],y1[25],y1[7]);
  k4[22]:=so22(y1[18],y1[19],y1[22],y1[26],y1[8]);
  k4[23]:=so23(y1[19],y1[23],y1[27],y1[9]);
  k4[24]:=so24(y1[20],y1[24],y1[10]);
  k4[25]:=so25(y1[21],y1[25],y1[11]);
  k4[26]:=so26(y1[22],y1[26],y1[12]);
  k4[27]:=so27(y1[23],y1[27],y1[13]);
  for i:=1 to 27 do
  y1[i]:=y[i]+h*(a51*k1[i]+a53*k3[i]+a54*k4[i]);
  r1:=1.0-y1[1]-y1[2]-y1[3]-y1[4]-y1[5]-y1[6]-y1[7]-y1[8]-y1[9]
  -y1[10]-y1[11]-y1[12]-y1[13]-y1[14]-y1[15]-y1[16]-y1[17]-y1[18]
  -y1[19]-y1[20]-y1[21]-y1[22]-y1[23]-y1[24]-y1[25]-y1[26]-y1[27];
  xx:=x+c5*h;
  if CheckBox1.Checked=false then
  if xx>tt0 then ag:=vspl(xx,tt0,tt1,agon) else ag:=0.0
  else ag:=0.0;
  if CheckBox1.Checked=false then
  if xx>tb0 then bg:=vspl(xx,tb0,tb1,bgon) else bg:=0.0
  else bg:=0.0;
  k5[1]:=so1(r1,y1[1],y1[3],y1[4],y1[15]);
  k5[2]:=so2(r1,y1[2],y1[4],y1[5],y1[16]);
  k5[3]:=so3(y1[1],y1[3],y1[6],y1[7],y1[17]);
  k5[4]:=so4(y1[1],y1[2],y1[4],y1[7],y1[8],y1[18]);
  k5[5]:=so5(y1[2],y1[5],y1[8],y1[9],y1[19]);
  k5[6]:=so6(y1[3],y1[6],y1[10],y1[20]);
  k5[7]:=so7(y1[3],y1[4],y1[7],y1[11],y1[21]);
  k5[8]:=so8(y1[4],y1[5],y1[8],y1[12],y1[22]);
  k5[9]:=so9(y1[5],y1[9],y1[13],y1[23]);
  k5[10]:=so10(y1[6],y1[10],y1[24]);
  k5[11]:=so11(y1[7],y1[11],y1[25]);
  k5[12]:=so12(y1[8],y1[12],y1[26]);
  k5[13]:=so13(y1[9],y1[13],y1[27]);
  k5[14]:=so14(r1,y1[14],y1[15],y1[16]);
  k5[15]:=so15(y1[14],y1[15],y1[17],y1[18],y1[1]);
  k5[16]:=so16(y1[14],y1[16],y1[18],y1[19],y1[2]);
  k5[17]:=so17(y1[15],y1[17],y1[20],y1[21],y1[3]);
  k5[18]:=so18(y1[15],y1[16],y1[18],y1[21],y1[22],y1[4]);
  k5[19]:=so19(y1[16],y1[19],y1[22],y1[23],y1[5]);
  k5[20]:=so20(y1[17],y1[20],y1[24],y1[6]);
  k5[21]:=so21(y1[17],y1[18],y1[21],y1[25],y1[7]);
  k5[22]:=so22(y1[18],y1[19],y1[22],y1[26],y1[8]);
  k5[23]:=so23(y1[19],y1[23],y1[27],y1[9]);
  k5[24]:=so24(y1[20],y1[24],y1[10]);
  k5[25]:=so25(y1[21],y1[25],y1[11]);
  k5[26]:=so26(y1[22],y1[26],y1[12]);
  k5[27]:=so27(y1[23],y1[27],y1[13]);
  for i:=1 to 27 do
  y1[i]:=y[i]+h*(a61*k1[i]+a64*k4[i]+a65*k5[i]);
  r1:=1.0-y1[1]-y1[2]-y1[3]-y1[4]-y1[5]-y1[6]-y1[7]-y1[8]-y1[9]
  -y1[10]-y1[11]-y1[12]-y1[13]-y1[14]-y1[15]-y1[16]-y1[17]-y1[18]
  -y1[19]-y1[20]-y1[21]-y1[22]-y1[23]-y1[24]-y1[25]-y1[26]-y1[27];
  xx:=x+c6*h;
  if CheckBox1.Checked=false then
  if xx>tt0 then ag:=vspl(xx,tt0,tt1,agon) else ag:=0.0
  else ag:=0.0;
  if CheckBox1.Checked=false then
  if xx>tb0 then bg:=vspl(xx,tb0,tb1,bgon) else bg:=0.0
  else bg:=0.0;
  k6[1]:=so1(r1,y1[1],y1[3],y1[4],y1[15]);
  k6[2]:=so2(r1,y1[2],y1[4],y1[5],y1[16]);
  k6[3]:=so3(y1[1],y1[3],y1[6],y1[7],y1[17]);
  k6[4]:=so4(y1[1],y1[2],y1[4],y1[7],y1[8],y1[18]);
  k6[5]:=so5(y1[2],y1[5],y1[8],y1[9],y1[19]);
  k6[6]:=so6(y1[3],y1[6],y1[10],y1[20]);
  k6[7]:=so7(y1[3],y1[4],y1[7],y1[11],y1[21]);
  k6[8]:=so8(y1[4],y1[5],y1[8],y1[12],y1[22]);
  k6[9]:=so9(y1[5],y1[9],y1[13],y1[23]);
  k6[10]:=so10(y1[6],y1[10],y1[24]);
  k6[11]:=so11(y1[7],y1[11],y1[25]);
  k6[12]:=so12(y1[8],y1[12],y1[26]);
  k6[13]:=so13(y1[9],y1[13],y1[27]);
  k6[14]:=so14(r1,y1[14],y1[15],y1[16]);
  k6[15]:=so15(y1[14],y1[15],y1[17],y1[18],y1[1]);
  k6[16]:=so16(y1[14],y1[16],y1[18],y1[19],y1[2]);
  k6[17]:=so17(y1[15],y1[17],y1[20],y1[21],y1[3]);
  k6[18]:=so18(y1[15],y1[16],y1[18],y1[21],y1[22],y1[4]);
  k6[19]:=so19(y1[16],y1[19],y1[22],y1[23],y1[5]);
  k6[20]:=so20(y1[17],y1[20],y1[24],y1[6]);
  k6[21]:=so21(y1[17],y1[18],y1[21],y1[25],y1[7]);
  k6[22]:=so22(y1[18],y1[19],y1[22],y1[26],y1[8]);
  k6[23]:=so23(y1[19],y1[23],y1[27],y1[9]);
  k6[24]:=so24(y1[20],y1[24],y1[10]);
  k6[25]:=so25(y1[21],y1[25],y1[11]);
  k6[26]:=so26(y1[22],y1[26],y1[12]);
  k6[27]:=so27(y1[23],y1[27],y1[13]);
  for i:=1 to 27 do
  y1[i]:=y[i]+h*(a71*k1[i]+a74*k4[i]+a75*k5[i]+a76*k6[i]);
  r1:=1.0-y1[1]-y1[2]-y1[3]-y1[4]-y1[5]-y1[6]-y1[7]-y1[8]-y1[9]
  -y1[10]-y1[11]-y1[12]-y1[13]-y1[14]-y1[15]-y1[16]-y1[17]-y1[18]
  -y1[19]-y1[20]-y1[21]-y1[22]-y1[23]-y1[24]-y1[25]-y1[26]-y1[27];
  xx:=x+c7*h;
  if CheckBox1.Checked=false then
  if xx>tt0 then ag:=vspl(xx,tt0,tt1,agon) else ag:=0.0
  else ag:=0.0;
  if CheckBox1.Checked=false then
  if xx>tb0 then bg:=vspl(xx,tb0,tb1,bgon) else bg:=0.0
  else bg:=0.0;
  k7[1]:=so1(r1,y1[1],y1[3],y1[4],y1[15]);
  k7[2]:=so2(r1,y1[2],y1[4],y1[5],y1[16]);
  k7[3]:=so3(y1[1],y1[3],y1[6],y1[7],y1[17]);
  k7[4]:=so4(y1[1],y1[2],y1[4],y1[7],y1[8],y1[18]);
  k7[5]:=so5(y1[2],y1[5],y1[8],y1[9],y1[19]);
  k7[6]:=so6(y1[3],y1[6],y1[10],y1[20]);
  k7[7]:=so7(y1[3],y1[4],y1[7],y1[11],y1[21]);
  k7[8]:=so8(y1[4],y1[5],y1[8],y1[12],y1[22]);
  k7[9]:=so9(y1[5],y1[9],y1[13],y1[23]);
  k7[10]:=so10(y1[6],y1[10],y1[24]);
  k7[11]:=so11(y1[7],y1[11],y1[25]);
  k7[12]:=so12(y1[8],y1[12],y1[26]);
  k7[13]:=so13(y1[9],y1[13],y1[27]);
  k7[14]:=so14(r1,y1[14],y1[15],y1[16]);
  k7[15]:=so15(y1[14],y1[15],y1[17],y1[18],y1[1]);
  k7[16]:=so16(y1[14],y1[16],y1[18],y1[19],y1[2]);
  k7[17]:=so17(y1[15],y1[17],y1[20],y1[21],y1[3]);
  k7[18]:=so18(y1[15],y1[16],y1[18],y1[21],y1[22],y1[4]);
  k7[19]:=so19(y1[16],y1[19],y1[22],y1[23],y1[5]);
  k7[20]:=so20(y1[17],y1[20],y1[24],y1[6]);
  k7[21]:=so21(y1[17],y1[18],y1[21],y1[25],y1[7]);
  k7[22]:=so22(y1[18],y1[19],y1[22],y1[26],y1[8]);
  k7[23]:=so23(y1[19],y1[23],y1[27],y1[9]);
  k7[24]:=so24(y1[20],y1[24],y1[10]);
  k7[25]:=so25(y1[21],y1[25],y1[11]);
  k7[26]:=so26(y1[22],y1[26],y1[12]);
  k7[27]:=so27(y1[23],y1[27],y1[13]);
  for i:=1 to 27 do
  y1[i]:=y[i]+h*(a81*k1[i]+a84*k4[i]+a85*k5[i]+a86*k6[i]+a87*k7[i]);
  r1:=1.0-y1[1]-y1[2]-y1[3]-y1[4]-y1[5]-y1[6]-y1[7]-y1[8]-y1[9]
  -y1[10]-y1[11]-y1[12]-y1[13]-y1[14]-y1[15]-y1[16]-y1[17]-y1[18]
  -y1[19]-y1[20]-y1[21]-y1[22]-y1[23]-y1[24]-y1[25]-y1[26]-y1[27];
  xx:=x+c8*h;
  if CheckBox1.Checked=false then
  if xx>tt0 then ag:=vspl(xx,tt0,tt1,agon) else ag:=0.0
  else ag:=0.0;
  if CheckBox1.Checked=false then
  if xx>tb0 then bg:=vspl(xx,tb0,tb1,bgon) else bg:=0.0
  else bg:=0.0;
  k2[1]:=so1(r1,y1[1],y1[3],y1[4],y1[15]);
  k2[2]:=so2(r1,y1[2],y1[4],y1[5],y1[16]);
  k2[3]:=so3(y1[1],y1[3],y1[6],y1[7],y1[17]);
  k2[4]:=so4(y1[1],y1[2],y1[4],y1[7],y1[8],y1[18]);
  k2[5]:=so5(y1[2],y1[5],y1[8],y1[9],y1[19]);
  k2[6]:=so6(y1[3],y1[6],y1[10],y1[20]);
  k2[7]:=so7(y1[3],y1[4],y1[7],y1[11],y1[21]);
  k2[8]:=so8(y1[4],y1[5],y1[8],y1[12],y1[22]);
  k2[9]:=so9(y1[5],y1[9],y1[13],y1[23]);
  k2[10]:=so10(y1[6],y1[10],y1[24]);
  k2[11]:=so11(y1[7],y1[11],y1[25]);
  k2[12]:=so12(y1[8],y1[12],y1[26]);
  k2[13]:=so13(y1[9],y1[13],y1[27]);
  k2[14]:=so14(r1,y1[14],y1[15],y1[16]);
  k2[15]:=so15(y1[14],y1[15],y1[17],y1[18],y1[1]);
  k2[16]:=so16(y1[14],y1[16],y1[18],y1[19],y1[2]);
  k2[17]:=so17(y1[15],y1[17],y1[20],y1[21],y1[3]);
  k2[18]:=so18(y1[15],y1[16],y1[18],y1[21],y1[22],y1[4]);
  k2[19]:=so19(y1[16],y1[19],y1[22],y1[23],y1[5]);
  k2[20]:=so20(y1[17],y1[20],y1[24],y1[6]);
  k2[21]:=so21(y1[17],y1[18],y1[21],y1[25],y1[7]);
  k2[22]:=so22(y1[18],y1[19],y1[22],y1[26],y1[8]);
  k2[23]:=so23(y1[19],y1[23],y1[27],y1[9]);
  k2[24]:=so24(y1[20],y1[24],y1[10]);
  k2[25]:=so25(y1[21],y1[25],y1[11]);
  k2[26]:=so26(y1[22],y1[26],y1[12]);
  k2[27]:=so27(y1[23],y1[27],y1[13]);
  for i:=1 to 27 do
  y1[i]:=y[i]+h*(a91*k1[i]+a94*k4[i]+a95*k5[i]+a96*k6[i]+a97*k7[i]+a98*k2[i]);
  r1:=1.0-y1[1]-y1[2]-y1[3]-y1[4]-y1[5]-y1[6]-y1[7]-y1[8]-y1[9]
  -y1[10]-y1[11]-y1[12]-y1[13]-y1[14]-y1[15]-y1[16]-y1[17]-y1[18]
  -y1[19]-y1[20]-y1[21]-y1[22]-y1[23]-y1[24]-y1[25]-y1[26]-y1[27];
  xx:=x+c9*h;
  if CheckBox1.Checked=false then
  if xx>tt0 then ag:=vspl(xx,tt0,tt1,agon) else ag:=0.0
  else ag:=0.0;
  if CheckBox1.Checked=false then
  if xx>tb0 then bg:=vspl(xx,tb0,tb1,bgon) else bg:=0.0
  else bg:=0.0;
  k3[1]:=so1(r1,y1[1],y1[3],y1[4],y1[15]);
  k3[2]:=so2(r1,y1[2],y1[4],y1[5],y1[16]);
  k3[3]:=so3(y1[1],y1[3],y1[6],y1[7],y1[17]);
  k3[4]:=so4(y1[1],y1[2],y1[4],y1[7],y1[8],y1[18]);
  k3[5]:=so5(y1[2],y1[5],y1[8],y1[9],y1[19]);
  k3[6]:=so6(y1[3],y1[6],y1[10],y1[20]);
  k3[7]:=so7(y1[3],y1[4],y1[7],y1[11],y1[21]);
  k3[8]:=so8(y1[4],y1[5],y1[8],y1[12],y1[22]);
  k3[9]:=so9(y1[5],y1[9],y1[13],y1[23]);
  k3[10]:=so10(y1[6],y1[10],y1[24]);
  k3[11]:=so11(y1[7],y1[11],y1[25]);
  k3[12]:=so12(y1[8],y1[12],y1[26]);
  k3[13]:=so13(y1[9],y1[13],y1[27]);
  k3[14]:=so14(r1,y1[14],y1[15],y1[16]);
  k3[15]:=so15(y1[14],y1[15],y1[17],y1[18],y1[1]);
  k3[16]:=so16(y1[14],y1[16],y1[18],y1[19],y1[2]);
  k3[17]:=so17(y1[15],y1[17],y1[20],y1[21],y1[3]);
  k3[18]:=so18(y1[15],y1[16],y1[18],y1[21],y1[22],y1[4]);
  k3[19]:=so19(y1[16],y1[19],y1[22],y1[23],y1[5]);
  k3[20]:=so20(y1[17],y1[20],y1[24],y1[6]);
  k3[21]:=so21(y1[17],y1[18],y1[21],y1[25],y1[7]);
  k3[22]:=so22(y1[18],y1[19],y1[22],y1[26],y1[8]);
  k3[23]:=so23(y1[19],y1[23],y1[27],y1[9]);
  k3[24]:=so24(y1[20],y1[24],y1[10]);
  k3[25]:=so25(y1[21],y1[25],y1[11]);
  k3[26]:=so26(y1[22],y1[26],y1[12]);
  k3[27]:=so27(y1[23],y1[27],y1[13]);
  for i:=1 to 27 do
  y1[i]:=y[i]+h*(a101*k1[i]+a104*k4[i]+a105*k5[i]+a106*k6[i]+a107*k7[i]+a108*k2[i]+a109*k3[i]);
{ Compute Imtermeniate Sums To Save Memory }
  for i:=1 to 27 do
  begin
    y11s:=a111*k1[i]+a114*k4[i]+a115*k5[i]+a116*k6[i]+a117*k7[i]+a118*k2[i]+a119*k3[i];
    y12s:=a121*k1[i]+a124*k4[i]+a125*k5[i]+a126*k6[i]+a127*k7[i]+a128*k2[i]+a129*k3[i];
    k4[i]:=a131*k1[i]+a134*k4[i]+a135*k5[i]+a136*k6[i]+a137*k7[i]+a138*k2[i]+a139*k3[i];
    k5[i]:=b1*k1[i]+b6*k6[i]+b7*k7[i]+b8*k2[i]+b9*k3[i];
    k6[i]:=bh1*k1[i]+bh6*k6[i]+bh7*k7[i]+bh8*k2[i]+bh9*k3[i];
    k2[i]:=y11s;
    k3[i]:=y12s;
  end;
{ The Last 4 stages }
  r1:=1.0-y1[1]-y1[2]-y1[3]-y1[4]-y1[5]-y1[6]-y1[7]-y1[8]-y1[9]
  -y1[10]-y1[11]-y1[12]-y1[13]-y1[14]-y1[15]-y1[16]-y1[17]-y1[18]
  -y1[19]-y1[20]-y1[21]-y1[22]-y1[23]-y1[24]-y1[25]-y1[26]-y1[27];
  xx:=x+c10*h;
  if CheckBox1.Checked=false then
  if xx>tt0 then ag:=vspl(xx,tt0,tt1,agon) else ag:=0.0
  else ag:=0.0;
  if CheckBox1.Checked=false then
  if xx>tb0 then bg:=vspl(xx,tb0,tb1,bgon) else bg:=0.0
  else bg:=0.0;
  k7[1]:=so1(r1,y1[1],y1[3],y1[4],y1[15]);
  k7[2]:=so2(r1,y1[2],y1[4],y1[5],y1[16]);
  k7[3]:=so3(y1[1],y1[3],y1[6],y1[7],y1[17]);
  k7[4]:=so4(y1[1],y1[2],y1[4],y1[7],y1[8],y1[18]);
  k7[5]:=so5(y1[2],y1[5],y1[8],y1[9],y1[19]);
  k7[6]:=so6(y1[3],y1[6],y1[10],y1[20]);
  k7[7]:=so7(y1[3],y1[4],y1[7],y1[11],y1[21]);
  k7[8]:=so8(y1[4],y1[5],y1[8],y1[12],y1[22]);
  k7[9]:=so9(y1[5],y1[9],y1[13],y1[23]);
  k7[10]:=so10(y1[6],y1[10],y1[24]);
  k7[11]:=so11(y1[7],y1[11],y1[25]);
  k7[12]:=so12(y1[8],y1[12],y1[26]);
  k7[13]:=so13(y1[9],y1[13],y1[27]);
  k7[14]:=so14(r1,y1[14],y1[15],y1[16]);
  k7[15]:=so15(y1[14],y1[15],y1[17],y1[18],y1[1]);
  k7[16]:=so16(y1[14],y1[16],y1[18],y1[19],y1[2]);
  k7[17]:=so17(y1[15],y1[17],y1[20],y1[21],y1[3]);
  k7[18]:=so18(y1[15],y1[16],y1[18],y1[21],y1[22],y1[4]);
  k7[19]:=so19(y1[16],y1[19],y1[22],y1[23],y1[5]);
  k7[20]:=so20(y1[17],y1[20],y1[24],y1[6]);
  k7[21]:=so21(y1[17],y1[18],y1[21],y1[25],y1[7]);
  k7[22]:=so22(y1[18],y1[19],y1[22],y1[26],y1[8]);
  k7[23]:=so23(y1[19],y1[23],y1[27],y1[9]);
  k7[24]:=so24(y1[20],y1[24],y1[10]);
  k7[25]:=so25(y1[21],y1[25],y1[11]);
  k7[26]:=so26(y1[22],y1[26],y1[12]);
  k7[27]:=so27(y1[23],y1[27],y1[13]);
  for i:=1 to 27 do
  y1[i]:=y[i]+h*(k2[i]+a1110*k7[i]);
  r1:=1.0-y1[1]-y1[2]-y1[3]-y1[4]-y1[5]-y1[6]-y1[7]-y1[8]-y1[9]
  -y1[10]-y1[11]-y1[12]-y1[13]-y1[14]-y1[15]-y1[16]-y1[17]-y1[18]
  -y1[19]-y1[20]-y1[21]-y1[22]-y1[23]-y1[24]-y1[25]-y1[26]-y1[27];
  xx:=x+c11*h;
  if CheckBox1.Checked=false then
  if xx>tt0 then ag:=vspl(xx,tt0,tt1,agon) else ag:=0.0
  else ag:=0.0;
  if CheckBox1.Checked=false then
  if xx>tb0 then bg:=vspl(xx,tb0,tb1,bgon) else bg:=0.0
  else bg:=0.0;
  k2[1]:=so1(r1,y1[1],y1[3],y1[4],y1[15]);
  k2[2]:=so2(r1,y1[2],y1[4],y1[5],y1[16]);
  k2[3]:=so3(y1[1],y1[3],y1[6],y1[7],y1[17]);
  k2[4]:=so4(y1[1],y1[2],y1[4],y1[7],y1[8],y1[18]);
  k2[5]:=so5(y1[2],y1[5],y1[8],y1[9],y1[19]);
  k2[6]:=so6(y1[3],y1[6],y1[10],y1[20]);
  k2[7]:=so7(y1[3],y1[4],y1[7],y1[11],y1[21]);
  k2[8]:=so8(y1[4],y1[5],y1[8],y1[12],y1[22]);
  k2[9]:=so9(y1[5],y1[9],y1[13],y1[23]);
  k2[10]:=so10(y1[6],y1[10],y1[24]);
  k2[11]:=so11(y1[7],y1[11],y1[25]);
  k2[12]:=so12(y1[8],y1[12],y1[26]);
  k2[13]:=so13(y1[9],y1[13],y1[27]);
  k2[14]:=so14(r1,y1[14],y1[15],y1[16]);
  k2[15]:=so15(y1[14],y1[15],y1[17],y1[18],y1[1]);
  k2[16]:=so16(y1[14],y1[16],y1[18],y1[19],y1[2]);
  k2[17]:=so17(y1[15],y1[17],y1[20],y1[21],y1[3]);
  k2[18]:=so18(y1[15],y1[16],y1[18],y1[21],y1[22],y1[4]);
  k2[19]:=so19(y1[16],y1[19],y1[22],y1[23],y1[5]);
  k2[20]:=so20(y1[17],y1[20],y1[24],y1[6]);
  k2[21]:=so21(y1[17],y1[18],y1[21],y1[25],y1[7]);
  k2[22]:=so22(y1[18],y1[19],y1[22],y1[26],y1[8]);
  k2[23]:=so23(y1[19],y1[23],y1[27],y1[9]);
  k2[24]:=so24(y1[20],y1[24],y1[10]);
  k2[25]:=so25(y1[21],y1[25],y1[11]);
  k2[26]:=so26(y1[22],y1[26],y1[12]);
  k2[27]:=so27(y1[23],y1[27],y1[13]);
  xph:=x+h;
  for i:=1 to 27 do
  y1[i]:=y[i]+h*(k3[i]+a1210*k7[i]+a1211*k2[i]);
  r1:=1.0-y1[1]-y1[2]-y1[3]-y1[4]-y1[5]-y1[6]-y1[7]-y1[8]-y1[9]
  -y1[10]-y1[11]-y1[12]-y1[13]-y1[14]-y1[15]-y1[16]-y1[17]-y1[18]
  -y1[19]-y1[20]-y1[21]-y1[22]-y1[23]-y1[24]-y1[25]-y1[26]-y1[27];
  xx:=xph;
  if CheckBox1.Checked=false then
  if xx>tt0 then ag:=vspl(xx,tt0,tt1,agon) else ag:=0.0
  else ag:=0.0;
  if CheckBox1.Checked=false then
  if xx>tb0 then bg:=vspl(xx,tb0,tb1,bgon) else bg:=0.0
  else bg:=0.0;
  k3[1]:=so1(r1,y1[1],y1[3],y1[4],y1[15]);
  k3[2]:=so2(r1,y1[2],y1[4],y1[5],y1[16]);
  k3[3]:=so3(y1[1],y1[3],y1[6],y1[7],y1[17]);
  k3[4]:=so4(y1[1],y1[2],y1[4],y1[7],y1[8],y1[18]);
  k3[5]:=so5(y1[2],y1[5],y1[8],y1[9],y1[19]);
  k3[6]:=so6(y1[3],y1[6],y1[10],y1[20]);
  k3[7]:=so7(y1[3],y1[4],y1[7],y1[11],y1[21]);
  k3[8]:=so8(y1[4],y1[5],y1[8],y1[12],y1[22]);
  k3[9]:=so9(y1[5],y1[9],y1[13],y1[23]);
  k3[10]:=so10(y1[6],y1[10],y1[24]);
  k3[11]:=so11(y1[7],y1[11],y1[25]);
  k3[12]:=so12(y1[8],y1[12],y1[26]);
  k3[13]:=so13(y1[9],y1[13],y1[27]);
  k3[14]:=so14(r1,y1[14],y1[15],y1[16]);
  k3[15]:=so15(y1[14],y1[15],y1[17],y1[18],y1[1]);
  k3[16]:=so16(y1[14],y1[16],y1[18],y1[19],y1[2]);
  k3[17]:=so17(y1[15],y1[17],y1[20],y1[21],y1[3]);
  k3[18]:=so18(y1[15],y1[16],y1[18],y1[21],y1[22],y1[4]);
  k3[19]:=so19(y1[16],y1[19],y1[22],y1[23],y1[5]);
  k3[20]:=so20(y1[17],y1[20],y1[24],y1[6]);
  k3[21]:=so21(y1[17],y1[18],y1[21],y1[25],y1[7]);
  k3[22]:=so22(y1[18],y1[19],y1[22],y1[26],y1[8]);
  k3[23]:=so23(y1[19],y1[23],y1[27],y1[9]);
  k3[24]:=so24(y1[20],y1[24],y1[10]);
  k3[25]:=so25(y1[21],y1[25],y1[11]);
  k3[26]:=so26(y1[22],y1[26],y1[12]);
  k3[27]:=so27(y1[23],y1[27],y1[13]);
  for i:=1 to 27 do
  y1[i]:=y[i]+h*(k4[i]+a1310*k7[i]+a1311*k2[i]);
  r1:=1.0-y1[1]-y1[2]-y1[3]-y1[4]-y1[5]-y1[6]-y1[7]-y1[8]-y1[9]
  -y1[10]-y1[11]-y1[12]-y1[13]-y1[14]-y1[15]-y1[16]-y1[17]-y1[18]
  -y1[19]-y1[20]-y1[21]-y1[22]-y1[23]-y1[24]-y1[25]-y1[26]-y1[27];
  k4[1]:=so1(r1,y1[1],y1[3],y1[4],y1[15]);
  k4[2]:=so2(r1,y1[2],y1[4],y1[5],y1[16]);
  k4[3]:=so3(y1[1],y1[3],y1[6],y1[7],y1[17]);
  k4[4]:=so4(y1[1],y1[2],y1[4],y1[7],y1[8],y1[18]);
  k4[5]:=so5(y1[2],y1[5],y1[8],y1[9],y1[19]);
  k4[6]:=so6(y1[3],y1[6],y1[10],y1[20]);
  k4[7]:=so7(y1[3],y1[4],y1[7],y1[11],y1[21]);
  k4[8]:=so8(y1[4],y1[5],y1[8],y1[12],y1[22]);
  k4[9]:=so9(y1[5],y1[9],y1[13],y1[23]);
  k4[10]:=so10(y1[6],y1[10],y1[24]);
  k4[11]:=so11(y1[7],y1[11],y1[25]);
  k4[12]:=so12(y1[8],y1[12],y1[26]);
  k4[13]:=so13(y1[9],y1[13],y1[27]);
  k4[14]:=so14(r1,y1[14],y1[15],y1[16]);
  k4[15]:=so15(y1[14],y1[15],y1[17],y1[18],y1[1]);
  k4[16]:=so16(y1[14],y1[16],y1[18],y1[19],y1[2]);
  k4[17]:=so17(y1[15],y1[17],y1[20],y1[21],y1[3]);
  k4[18]:=so18(y1[15],y1[16],y1[18],y1[21],y1[22],y1[4]);
  k4[19]:=so19(y1[16],y1[19],y1[22],y1[23],y1[5]);
  k4[20]:=so20(y1[17],y1[20],y1[24],y1[6]);
  k4[21]:=so21(y1[17],y1[18],y1[21],y1[25],y1[7]);
  k4[22]:=so22(y1[18],y1[19],y1[22],y1[26],y1[8]);
  k4[23]:=so23(y1[19],y1[23],y1[27],y1[9]);
  k4[24]:=so24(y1[20],y1[24],y1[10]);
  k4[25]:=so25(y1[21],y1[25],y1[11]);
  k4[26]:=so26(y1[22],y1[26],y1[12]);
  k4[27]:=so27(y1[23],y1[27],y1[13]);
{  nfcn:=nfcn+13;}
  for i:=1 to 27 do
  begin
    k5[i]:=y[i]+h*(k5[i]+b10*k7[i]+b11*k2[i]+b12*k3[i]+b13*k4[i]);
    k6[i]:=y[i]+h*(k6[i]+bh10*k7[i]+bh11*k2[i]+bh12*k3[i]);
  end;
{ Error estimation }
  Err:=0;
  for i:=1 to 27 do
  begin
    denom:=max(max(1e-6,abs(k5[i])),max(abs(y[i]),2.0*uround/eps));
    err:=err+sqr((k5[i]-k6[i])/denom);
  end;
  err:=sqrt(err/(5.0));
{ Computation HNEW }
  fac:=max((1.0/6.0),min(3.0,Power(err/eps,1.0/8.0)/0.9));
  hnew:=h/fac; nstep:=nstep+1;
  if err>eps then goto 51;
{ Step is accepted }
{  inc(naccpt);}
  for i:=1 to 27 do
  y[i]:=k5[i];
  yh[0]:=1.0-y[1]-y[2]-y[3]-y[4]-y[5]-y[6]-y[7]-y[8]-y[9]
  -y[10]-y[11]-y[12]-y[13]-y[14]-y[15]-y[16]-y[17]-y[18]
  -y[19]-y[20]-y[21]-y[22]-y[23]-y[24]-y[25]-y[26]-y[27];
  for ii:=1 to 27 do yh[ii]:=y[ii];
  x:=xph;
  if CheckBox1.Checked=true then ddp:=500.0*dp
  else
if agon>=bgon then
  if x<tt0 then
  if x<2.0 then ddp:=10.0*dp else ddp:=100.0*dp
  else
  if x<(tt0+tt1) then if x<(tt0+6.0) then ddp:=dp else ddp:=10.0*dp
  else if x<(tt0+6.0) then ddp:=dp else ddp:=100.0*dp
else
  if x<tb0 then
  if x<2.0 then ddp:=10.0*dp else ddp:=100.0*dp
  else
  if x<(tb0+tb1) then if x<(tb0+6.0) then ddp:=dp else ddp:=10.0*dp
  else if x<(tb0+6.0) then ddp:=dp else ddp:=100.0*dp;

  if (x-xold)>ddp then
  begin
{Writeln(fff,'ag = ',yh[8]:10:5);
for ii:=0 to 7 do
Write(fff,yh[ii]:10:5);
Writeln(fff);}
  xold:=xold+ddp;
  jj:=jj+1;
  ProgressBar1.StepBy(1);
  tx[jj]:=x;
  rr[jj]:=yh[0];
  ar[jj]:=y[1];
  br[jj]:=y[2];
  a2r[jj]:=y[3];
  abr[jj]:=y[4];
  b2r[jj]:=y[5];
  a3r[jj]:=y[6];
  a2br[jj]:=y[7];
  ab2r[jj]:=y[8];
  b3r[jj]:=y[9];
  a3ro[jj]:=y[10];
  a2bro[jj]:=y[11];
  ab2ro[jj]:=y[12];
  b3ro[jj]:=y[13];
  d[jj]:=y[14];
  ad[jj]:=y[15];
  bd[jj]:=y[16];
  a2d[jj]:=y[17];
  abd[jj]:=y[18];
  b2d[jj]:=y[19];
  a3d[jj]:=y[20];
  a2bd[jj]:=y[21];
  ab2d[jj]:=y[22];
  b3d[jj]:=y[23];
  a3do[jj]:=y[24];
  a2bdo[jj]:=y[25];
  ab2do[jj]:=y[26];
  b3do[jj]:=y[27];
sumr[jj]:=rr[jj]+ar[jj]+br[jj]+a2r[jj]+abr[jj]+b2r[jj]
+a3r[jj]+a2br[jj]+ab2r[jj]+b3r[jj];
sumo[jj]:=a3ro[jj]+a2bro[jj]+ab2ro[jj]+b3ro[jj];
sumd[jj]:=d[jj]+ad[jj]+bd[jj]+a2d[jj]+abd[jj]+b2d[jj]+a3d[jj]+a2bd[jj]
+ab2d[jj]+b3d[jj]+a3do[jj]+a2bdo[jj]+ab2do[jj]+b3do[jj];
  {writeln(fff,tx[jj]:15:5,a3ro[jj]:10:5,rm[jj]:10:5);
    Assignfile(fff,'otl.dat');
  Append(fff);
  Writeln(fff,tx[jj]:10:3,a2ro[jj]:8:3,y[5]:8:3);
  CloseFile(fff);}
  end;
  if abs(hnew)>hmax then hnew:=posneg*hmax;
  if reject then hnew:=posneg*min(abs(hnew),abs(h));
  reject:=false;
  h:=hnew;
  if (xend-x)<hnew then goto 79;
  goto 1;
51:
  reject:=true;
  h:=hnew;
{  if naccpt>=1 then inc(nreject);}
{  dec(nfcn);}
 if nstep<nmax then goto 2;
79:
{     writeln(fff,'end ',j:5,jj:8);}
        {end of Runge-Kutta calc}
        NJ:=jj; am:=-10.0; ta1:=-10.0; ta2:=-10.0; rt:=-10.0;
        if agon>=bgon then
        begin aka:=1.0; bka:=0.2; end
        else
        begin bka:=1.0; aka:=0.2; end;
         for ii:=1 to jj do
         begin
         Series3.AddXY(tx[ii],sumo[ii],'',clBlue);
         if tx[ii]>tt0 then ag:=vspl(tx[ii],tt0,tt1,aka) else ag:=0.0;
         Series2.AddXY(tx[ii],ag,'',clYellow);
         if tx[ii]>tb0 then bg:=vspl(tx[ii],tb0,tb1,bka) else bg:=0.0;
         Series5.AddXY(tx[ii],bg,'',clAqua);
         Series1.AddXY(tx[ii],sumr[ii],'',clRed);
         Series4.AddXY(tx[ii],sumd[ii],'',clLime);
         {Series4.AddXY(tx[ii],ar[ii],'',clFuchsia);
         Series5.AddXY(tx[ii],a2r[ii],'',clAqua);
         Series6.AddXY(tx[ii],a3r[ii],'',clOlive);}
          end;
if CheckBox1.Checked=false then
begin
         im:=0; imi:=jj-1;
         for ii:=1 to jj do
{         writeln(fff,j:5,ii:8,tx[ii]:15:5,a2ro[ii]:8:5);}
         if sumo[ii]>am then begin am:=sumo[ii]; im:=ii; end;
         if agon>=bgon then
         begin
         for ii:=2 to jj-1 do
         if (tx[ii-1]<(tt0+tt1)) AND (tx[ii+1]>(tt0+tt1)) then imi:=ii
         end
         else
         begin
         for ii:=2 to jj-1 do
         if (tx[ii-1]<(tb0+tb1)) AND (tx[ii+1]>(tb0+tb1)) then imi:=ii
         end;
         ami:=sumo[imi]; amo:=am-ami;
         for ii:=2 to im-1 do
         if (sumo[ii-1]<0.1*am) AND (sumo[ii+1]>0.1*am) then t1:=tx[ii];
         for ii:=2 to im-1 do
         if (sumo[ii-1]<0.9*am) AND (sumo[ii+1]>0.9*am) then t9:=tx[ii];
         rt:=t9-t1;

         for ii:=im+1 to jj-1 do
         if ((sumo[ii-1]-ami)>0.2*amo) AND ((sumo[ii+1]-ami)<0.2*amo) then
         begin t2:=tx[ii]; a2:=sumo[ii]; end;
         if agon>=bgon then
         if (t2<(tt0+tt1)) AND (a2>ami) AND (ami>0.0) then
         ta2:=(t2-tt0-tt1)/ln(ami/a2)
         else ta2:=-1.0
         else
         if (t2<(tb0+tb1)) AND (a2>ami) AND (ami>0.0) then
         ta2:=(t2-tb0-tb1)/ln(ami/a2)
         else ta2:=-1.0;

         for ii:=im+1 to jj-1 do
         if ((sumo[ii-1]-ami)>0.8*amo) AND ((sumo[ii+1]-ami)<0.8*amo) then
         begin t8:=tx[ii]; a8:=sumo[ii]; end;
         for ii:=im+1 to jj-1 do
         if ((sumo[ii-1]-ami)>0.6*amo) AND ((sumo[ii+1]-ami)<0.6*amo) then
         begin t5:=tx[ii]; a5:=sumo[ii]; end;
         if agon>=bgon then
         begin
         a8:=a8-ami*exp((tt0+tt1-t8)/ta2);
         a5:=a5-ami*exp((tt0+tt1-t5)/ta2);
         end
         else
         begin
         a8:=a8-ami*exp((tb0+tb1-t8)/ta2);
         a5:=a5-ami*exp((tb0+tb1-t5)/ta2);
         end;
         if (t8>0.0) AND (t5>t8) AND (a5>0.0) AND (a8>a5) then
         ta1:=(t8-t5)/ln(a5/a8)
         else ta1:=-1.0;
end;
{  Assignfile(fff,'otl.dat');
  Append(fff); }
99:  {Writeln(fff,j:5,am[j]:10:2,rt[j]:10:2,ta[j]:10:2)};
{  CloseFile(fff);}

Edit13.Text:=FloatToStrF(yh[0],ffFixed,7,4);
Edit14.Text:=FloatToStrF(yh[1],ffFixed,7,4);
Edit15.Text:=FloatToStrF(yh[2],ffFixed,7,4);
Edit16.Text:=FloatToStrF(yh[3],ffFixed,7,4);
Edit23.Text:=FloatToStrF(yh[4],ffFixed,7,4);
Edit30.Text:=FloatToStrF(yh[5],ffFixed,7,4);
Edit33.Text:=FloatToStrF(yh[6],ffFixed,7,4);
Edit34.Text:=FloatToStrF(yh[7],ffFixed,7,4);
Edit35.Text:=FloatToStrF(yh[8],ffFixed,7,4);
Edit36.Text:=FloatToStrF(yh[9],ffFixed,7,4);
Edit41.Text:=FloatToStrF(yh[10],ffFixed,7,4);
Edit43.Text:=FloatToStrF(yh[11],ffFixed,7,4);
Edit45.Text:=FloatToStrF(yh[12],ffFixed,7,4);
Edit52.Text:=FloatToStrF(yh[13],ffFixed,7,4);
Edit77.Text:=FloatToStrF(yh[14],ffFixed,7,4);
Edit78.Text:=FloatToStrF(yh[15],ffFixed,7,4);
Edit79.Text:=FloatToStrF(yh[16],ffFixed,7,4);
Edit80.Text:=FloatToStrF(yh[17],ffFixed,7,4);
Edit81.Text:=FloatToStrF(yh[18],ffFixed,7,4);
Edit82.Text:=FloatToStrF(yh[19],ffFixed,7,4);
Edit83.Text:=FloatToStrF(yh[20],ffFixed,7,4);
Edit84.Text:=FloatToStrF(yh[21],ffFixed,7,4);
Edit85.Text:=FloatToStrF(yh[22],ffFixed,7,4);
Edit86.Text:=FloatToStrF(yh[23],ffFixed,7,4);
Edit87.Text:=FloatToStrF(yh[24],ffFixed,7,4);
Edit88.Text:=FloatToStrF(yh[25],ffFixed,7,4);
Edit89.Text:=FloatToStrF(yh[26],ffFixed,7,4);
Edit90.Text:=FloatToStrF(yh[27],ffFixed,7,4);
Edit91.Text:=FloatToStrF(sumr[jj],ffFixed,7,4);
Edit92.Text:=FloatToStrF(sumo[jj],ffFixed,7,4);
Edit93.Text:=FloatToStrF(sumd[jj],ffFixed,7,4);

Edit24.Text:=FloatToStrF(am,ffFixed,7,4);
Edit50.Text:=FloatToStrF(100.0*am/mam,ffFixed,6,1);
Edit25.Text:=FloatToStrF(ta1,ffFixed,7,4);
Edit44.Text:=FloatToStrF(ta2,ffFixed,7,4);
Edit32.Text:=FloatToStrF(rt,ffFixed,7,4);
{CloseFile(fff);}
end;

procedure TForm1.Button5Click(Sender: TObject);
begin
rish:=StrToFloat(Edit13.Text);
arish:=StrToFloat(Edit14.Text);
brish:=StrToFloat(Edit15.Text);
a2rish:=StrToFloat(Edit16.Text);
abrish:=StrToFloat(Edit23.Text);
b2rish:=StrToFloat(Edit30.Text);
a3rish:=StrToFloat(Edit33.Text);
a2brish:=StrToFloat(Edit34.Text);
ab2rish:=StrToFloat(Edit35.Text);
b3rish:=StrToFloat(Edit36.Text);
a3roish:=StrToFloat(Edit41.Text);
a2broish:=StrToFloat(Edit43.Text);
ab2roish:=StrToFloat(Edit45.Text);
b3roish:=StrToFloat(Edit52.Text);
dish:=StrToFloat(Edit77.Text);
adish:=StrToFloat(Edit78.Text);
bdish:=StrToFloat(Edit79.Text);
a2dish:=StrToFloat(Edit80.Text);
abdish:=StrToFloat(Edit81.Text);
b2dish:=StrToFloat(Edit82.Text);
a3dish:=StrToFloat(Edit83.Text);
a2bdish:=StrToFloat(Edit84.Text);
ab2dish:=StrToFloat(Edit85.Text);
b3dish:=StrToFloat(Edit86.Text);
a3doish:=StrToFloat(Edit87.Text);
a2bdoish:=StrToFloat(Edit88.Text);
ab2doish:=StrToFloat(Edit89.Text);
b3doish:=StrToFloat(Edit90.Text);
end;

procedure TForm1.Button6Click(Sender: TObject);
var su:real;
begin
Edit13.Text:=FloatToStr(rish);
Edit14.Text:=FloatToStr(arish);
Edit15.Text:=FloatToStr(brish);
Edit16.Text:=FloatToStr(a2rish);
Edit23.Text:=FloatToStr(abrish);
Edit30.Text:=FloatToStr(b2rish);
Edit33.Text:=FloatToStr(a3rish);
Edit34.Text:=FloatToStr(a2brish);
Edit35.Text:=FloatToStr(ab2rish);
Edit36.Text:=FloatToStr(b3rish);
Edit41.Text:=FloatToStr(a3roish);
Edit43.Text:=FloatToStr(a2broish);
Edit45.Text:=FloatToStr(ab2roish);
Edit52.Text:=FloatToStr(b3roish);
Edit77.Text:=FloatToStr(dish);
Edit78.Text:=FloatToStr(adish);
Edit79.Text:=FloatToStr(bdish);
Edit80.Text:=FloatToStr(a2dish);
Edit81.Text:=FloatToStr(abdish);
Edit82.Text:=FloatToStr(b2dish);
Edit83.Text:=FloatToStr(a3dish);
Edit84.Text:=FloatToStr(a2bdish);
Edit85.Text:=FloatToStr(ab2dish);
Edit86.Text:=FloatToStr(b3dish);
Edit87.Text:=FloatToStr(a3doish);
Edit88.Text:=FloatToStr(a2bdoish);
Edit89.Text:=FloatToStr(ab2doish);
Edit90.Text:=FloatToStr(b3doish);
su:=rish+arish+brish+a2rish+abrish+b2rish+a3rish+a2brish+ab2rish+b3rish;
Edit91.Text:=FloatToStr(su);
su:=a3roish+a2broish+ab2roish+b3roish;
Edit92.Text:=FloatToStr(su);
su:=dish+adish+bdish+a2dish+abdish+b2dish+a3dish+a2bdish+ab2dish+b3dish+
a3doish+a2bdoish+ab2doish+b3doish;
Edit93.Text:=FloatToStr(su);
end;

procedure TForm1.Button8Click(Sender: TObject);
var foo:TextFile;
pao:real;
pmo,pjo,pko:integer;
pps,ppsn:string;
begin
ForceCurrentDirectory:=True;
SaveDialog1.FileName:=Edit48.Text;
SaveDialog1.Filter:='Parameter files(*.tpt)|*.tpt';
if SaveDialog1.Execute then
begin
pps:=SaveDialog1.FileName;
pko:=Length(pps);
pmo:=1;
for pjo:=1 to pko do if pps[pjo]='\' then pmo:=pjo;
ppsn:=pps;
Delete(ppsn,1,pmo);
AssignFile(foo,pps);
Rewrite(foo);
pao:=StrToFloat(Edit1.Text);
writeln(foo,pao);
pao:=StrToFloat(Edit2.Text);
writeln(foo,pao);
pao:=StrToFloat(Edit3.Text);
writeln(foo,pao);
pao:=StrToFloat(Edit4.Text);
writeln(foo,pao);
pao:=StrToFloat(Edit5.Text);
writeln(foo,pao);
pao:=StrToFloat(Edit6.Text);
writeln(foo,pao);
pao:=StrToFloat(Edit7.Text);
writeln(foo,pao);
pao:=StrToFloat(Edit8.Text);
writeln(foo,pao);
pao:=StrToFloat(Edit19.Text);
writeln(foo,pao);
pao:=StrToFloat(Edit20.Text);
writeln(foo,pao);
pao:=StrToFloat(Edit21.Text);
writeln(foo,pao);
pao:=StrToFloat(Edit22.Text);
writeln(foo,pao);
pao:=StrToFloat(Edit26.Text);
writeln(foo,pao);
pao:=StrToFloat(Edit27.Text);
writeln(foo,pao);
pao:=StrToFloat(Edit28.Text);
writeln(foo,pao);
pao:=StrToFloat(Edit29.Text);
writeln(foo,pao);
pao:=StrToFloat(Edit37.Text);
writeln(foo,pao);
pao:=StrToFloat(Edit38.Text);
writeln(foo,pao);
pao:=StrToFloat(Edit39.Text);
writeln(foo,pao);
pao:=StrToFloat(Edit40.Text);
writeln(foo,pao);
pao:=StrToFloat(Edit51.Text);
writeln(foo,pao);
pao:=StrToFloat(Edit12.Text);
writeln(foo,pao);
pao:=StrToFloat(Edit17.Text);
writeln(foo,pao);
pao:=StrToFloat(Edit18.Text);
writeln(foo,pao);
pao:=StrToFloat(Edit53.Text);
writeln(foo,pao);
pao:=StrToFloat(Edit54.Text);
writeln(foo,pao);
pao:=StrToFloat(Edit55.Text);
writeln(foo,pao);
pao:=StrToFloat(Edit56.Text);
writeln(foo,pao);
pao:=StrToFloat(Edit57.Text);
writeln(foo,pao);
pao:=StrToFloat(Edit58.Text);
writeln(foo,pao);
pao:=StrToFloat(Edit59.Text);
writeln(foo,pao);
pao:=StrToFloat(Edit60.Text);
writeln(foo,pao);
pao:=StrToFloat(Edit61.Text);
writeln(foo,pao);
pao:=StrToFloat(Edit62.Text);
writeln(foo,pao);
pao:=StrToFloat(Edit63.Text);
writeln(foo,pao);
pao:=StrToFloat(Edit64.Text);
writeln(foo,pao);
pao:=StrToFloat(Edit65.Text);
writeln(foo,pao);
pao:=StrToFloat(Edit66.Text);
writeln(foo,pao);
pao:=StrToFloat(Edit67.Text);
writeln(foo,pao);
pao:=StrToFloat(Edit68.Text);
writeln(foo,pao);
pao:=StrToFloat(Edit69.Text);
writeln(foo,pao);
pao:=StrToFloat(Edit70.Text);
writeln(foo,pao);
pao:=StrToFloat(Edit71.Text);
writeln(foo,pao);
pao:=StrToFloat(Edit72.Text);
writeln(foo,pao);
pao:=StrToFloat(Edit94.Text);
writeln(foo,pao);
pao:=StrToFloat(Edit95.Text);
writeln(foo,pao);
pao:=StrToFloat(Edit96.Text);
writeln(foo,pao);
pao:=StrToFloat(Edit97.Text);
writeln(foo,pao);
pao:=StrToFloat(Edit98.Text);
writeln(foo,pao);
pao:=StrToFloat(Edit99.Text);
writeln(foo,pao);
pao:=StrToFloat(Edit100.Text);
writeln(foo,pao);
pao:=StrToFloat(Edit101.Text);
writeln(foo,pao);
pao:=StrToFloat(Edit73.Text);
writeln(foo,pao);
pao:=StrToFloat(Edit74.Text);
writeln(foo,pao);
pao:=StrToFloat(Edit75.Text);
writeln(foo,pao);
pao:=StrToFloat(Edit76.Text);
writeln(foo,pao);
pao:=StrToFloat(Edit9.Text);
writeln(foo,pao);
pao:=StrToFloat(Edit10.Text);
writeln(foo,pao);
pao:=StrToFloat(Edit31.Text);
writeln(foo,pao);
pao:=StrToFloat(Edit42.Text);
writeln(foo,pao);
pao:=StrToFloat(Edit46.Text);
writeln(foo,pao);
pao:=StrToFloat(Edit47.Text);
writeln(foo,pao);
Edit48.Text:=ppsn;
CloseFile(foo);
end;

end;

procedure TForm1.Button7Click(Sender: TObject);
var fpp:TextFile;
pkp,pmp,pjp:integer;
ppr:real;
psp,pspn:string;
begin
ForceCurrentDirectory:=True;
OpenDialog1.DefaultExt:='tpt';
OpenDialog1.Filter:='Parameter files(*.tpt)|*.tpt';
if OpenDialog1.Execute then
begin
psp:=OpenDialog1.FileName;
pkp:=Length(psp);
pmp:=1;
for pjp:=1 to pkp do if psp[pjp]='\' then pmp:=pjp;
pspn:=psp;
Delete(pspn,1,pmp);
AssignFile(fpp,psp);
Reset(fpp);
readln(fpp,ppr);
Edit1.Text:=FloatToStrF(ppr,ffFixed,9,1);
readln(fpp,ppr);
Edit2.Text:=FloatToStrF(ppr,ffFixed,9,1);
readln(fpp,ppr);
Edit3.Text:=FloatToStrF(ppr,ffFixed,9,1);
readln(fpp,ppr);
Edit4.Text:=FloatToStrF(ppr,ffFixed,9,1);
readln(fpp,ppr);
Edit5.Text:=FloatToStrF(ppr,ffFixed,9,1);
readln(fpp,ppr);
Edit6.Text:=FloatToStrF(ppr,ffFixed,9,5);
readln(fpp,ppr);
Edit7.Text:=FloatToStrF(ppr,ffFixed,9,1);
readln(fpp,ppr);
Edit8.Text:=FloatToStrF(ppr,ffFixed,9,5);
readln(fpp,ppr);
Edit19.Text:=FloatToStrF(ppr,ffFixed,9,2);
readln(fpp,ppr);
Edit20.Text:=FloatToStrF(ppr,ffFixed,9,4);
readln(fpp,ppr);
Edit21.Text:=FloatToStrF(ppr,ffFixed,9,2);
readln(fpp,ppr);
Edit22.Text:=FloatToStrF(ppr,ffFixed,9,4);
readln(fpp,ppr);
Edit26.Text:=FloatToStrF(ppr,ffFixed,9,2);
readln(fpp,ppr);
Edit27.Text:=FloatToStrF(ppr,ffFixed,9,4);
readln(fpp,ppr);
Edit28.Text:=FloatToStrF(ppr,ffFixed,9,2);
readln(fpp,ppr);
Edit29.Text:=FloatToStrF(ppr,ffFixed,9,4);
readln(fpp,ppr);
Edit37.Text:=FloatToStrF(ppr,ffFixed,9,2);
readln(fpp,ppr);
Edit38.Text:=FloatToStrF(ppr,ffFixed,9,4);
readln(fpp,ppr);
Edit39.Text:=FloatToStrF(ppr,ffFixed,9,2);
readln(fpp,ppr);
Edit40.Text:=FloatToStrF(ppr,ffFixed,9,4);
readln(fpp,ppr);
Edit51.Text:=FloatToStrF(ppr,ffFixed,9,2);
readln(fpp,ppr);
Edit12.Text:=FloatToStrF(ppr,ffFixed,9,4);
readln(fpp,ppr);
Edit17.Text:=FloatToStrF(ppr,ffFixed,9,2);
readln(fpp,ppr);
Edit18.Text:=FloatToStrF(ppr,ffFixed,9,4);
readln(fpp,ppr);
Edit53.Text:=FloatToStrF(ppr,ffFixed,9,5);
readln(fpp,ppr);
Edit54.Text:=FloatToStrF(ppr,ffFixed,9,5);
readln(fpp,ppr);
Edit55.Text:=FloatToStrF(ppr,ffFixed,9,5);
readln(fpp,ppr);
Edit56.Text:=FloatToStrF(ppr,ffFixed,9,5);
readln(fpp,ppr);
Edit57.Text:=FloatToStrF(ppr,ffFixed,9,5);
readln(fpp,ppr);
Edit58.Text:=FloatToStrF(ppr,ffFixed,9,5);
readln(fpp,ppr);
Edit59.Text:=FloatToStrF(ppr,ffFixed,9,5);
readln(fpp,ppr);
Edit60.Text:=FloatToStrF(ppr,ffFixed,9,5);
readln(fpp,ppr);
Edit61.Text:=FloatToStrF(ppr,ffFixed,9,5);
readln(fpp,ppr);
Edit62.Text:=FloatToStrF(ppr,ffFixed,9,5);
readln(fpp,ppr);
Edit63.Text:=FloatToStrF(ppr,ffFixed,9,5);
readln(fpp,ppr);
Edit64.Text:=FloatToStrF(ppr,ffFixed,9,5);
readln(fpp,ppr);
Edit65.Text:=FloatToStrF(ppr,ffFixed,9,5);
readln(fpp,ppr);
Edit66.Text:=FloatToStrF(ppr,ffFixed,9,5);
readln(fpp,ppr);
Edit67.Text:=FloatToStrF(ppr,ffFixed,9,5);
readln(fpp,ppr);
Edit68.Text:=FloatToStrF(ppr,ffFixed,9,5);
readln(fpp,ppr);
Edit69.Text:=FloatToStrF(ppr,ffFixed,9,5);
readln(fpp,ppr);
Edit70.Text:=FloatToStrF(ppr,ffFixed,9,5);
readln(fpp,ppr);
Edit71.Text:=FloatToStrF(ppr,ffFixed,9,5);
readln(fpp,ppr);
Edit72.Text:=FloatToStrF(ppr,ffFixed,9,5);
readln(fpp,ppr);
Edit94.Text:=FloatToStrF(ppr,ffFixed,9,5);
readln(fpp,ppr);
Edit95.Text:=FloatToStrF(ppr,ffFixed,9,5);
readln(fpp,ppr);
Edit96.Text:=FloatToStrF(ppr,ffFixed,9,5);
readln(fpp,ppr);
Edit97.Text:=FloatToStrF(ppr,ffFixed,9,5);
readln(fpp,ppr);
Edit98.Text:=FloatToStrF(ppr,ffFixed,9,5);
readln(fpp,ppr);
Edit99.Text:=FloatToStrF(ppr,ffFixed,9,5);
readln(fpp,ppr);
Edit100.Text:=FloatToStrF(ppr,ffFixed,9,5);
readln(fpp,ppr);
Edit101.Text:=FloatToStrF(ppr,ffFixed,9,5);
readln(fpp,ppr);
Edit73.Text:=FloatToStrF(ppr,ffFixed,9,7);
readln(fpp,ppr);
Edit74.Text:=FloatToStrF(ppr,ffFixed,7,2);
readln(fpp,ppr);
Edit75.Text:=FloatToStrF(ppr,ffFixed,7,2);
readln(fpp,ppr);
Edit76.Text:=FloatToStrF(ppr,ffFixed,7,2);
readln(fpp,ppr);
Edit9.Text:=FloatToStrF(ppr,ffFixed,9,7);
readln(fpp,ppr);
Edit10.Text:=FloatToStrF(ppr,ffFixed,7,2);
readln(fpp,ppr);
Edit31.Text:=FloatToStrF(ppr,ffFixed,7,2);
readln(fpp,ppr);
Edit42.Text:=FloatToStrF(ppr,ffFixed,7,2);
readln(fpp,ppr);
Edit46.Text:=FloatToStrF(ppr,ffFixed,8,4);
readln(fpp,ppr);
Edit47.Text:=FloatToStrF(ppr,ffFixed,8,4);
Edit48.Text:=pspn;
CloseFile(fpp);
end;

end;

end.

Loading data, please wait...