[incarnus]

log(p)=-A/T+B+C*log(T)+10^-3*D*T
надо отсюда найти Т(температуру) при p=1,33-все коэффициенты известны и введены в массивы для 70 элементов

[ART]

Каким именно методом? И где эти коэффициенты известные???

[incarnus]

я не знаю каким методом..просто это уравнение решается разве чис методами и никак иначе....
а коэффициенты могу скинуть!

[incarnus]

Код:

//ag
a[0]:=14710;
b[0]:=11.66;
c[0]:=-0.755;
d[0]:=0;
tmin[0]:=298;
tk[0]:=1234;
tmax[0]:=2400;
a1[0]:=14260;
b1[0]:=12.23;
c1[0]:=-1.055;
d1[0]:=0;
tkip[0]:=954.493;
//al
a[1]:=16450;
b[1]:=12.36;
c[1]:=-1.023;
d[1]:=0;
tmin[1]:=1200;
tmax[1]:=2800;
tk[1]:=2800;
tkip[1]:=1078.077;
//am
a[2]:=13700;
b[2]:=13.97;
c[2]:=-1.0;
d[2]:=0;
tmin[2]:=1103;
tmax[2]:=1453;
tk[2]:=1453;
tkip[2]:=802.743;
//as4
a[3]:=6160;
b[3]:=9.82;
c[3]:=0;
d[3]:=0;
tmin[3]:=600;
tmax[3]:=900;
tk[3]:=900;
tkip[3]:=389.355;
//au
a[4]:=19820;
b[4]:=10.81;
c[4]:=-0.306;
d[4]:=-0.16;
tmin[4]:=298;
tmax[4]:=3240;
tk[4]:=1336;
a1[4]:=19280;
b1[4]:=12.38;
c1[4]:=-1.01;
d1[4]:=0;
tkip[4]:=1265.708;
//b
a[5]:=29900;
b[5]:=13.88;
c[5]:=-1.0;
d[5]:=0;
tmin[5]:=1000;
tmax[5]:=2573;
tk[5]:=2573;
tkip[5]:=1798.370;
//ba
a[6]:=9730;
b[6]:=7.83;
c[6]:=0;
d[6]:=0;
tmin[6]:=750;
tmax[6]:=1200;
tk[6]:=983;
a1[6]:=9340;
b1[6]:=7.42;
c1[6]:=0;
d1[6]:=0;
tkip[6]:=703.490;
//be
a[7]:=10734;
b[7]:=9.067;
c[7]:=0;
d[7]:=-0.145;
tmin[7]:=900;
tmax[7]:=1557;
tk[7]:=1557;
tkip[7]:=717.320;
//bi
a[8]:=10400;
b[8]:=12.35;
c[8]:=-1.26;
d[8]:=0;
tmin[8]:=545;
tmax[8]:=1833;
tk[8]:=1833;
tkip[8]:=704.476;
//bi2
a[9]:=10730;
b[9]:=18.1;
c[9]:=-3.02;
d[9]:=0;
tmin[9]:=545;
tmax[9]:=1833;
tk[9]:=1833;
tkip[9]:=691.116 ;
//ca
a[10]:=10300;
b[10]:=14.97;
c[10]:=-1.67;
d[10]:=0;
tmin[10]:=713;
tmax[10]:=1757;
tk[10]:=1112;
a1[10]:=9600;
b1[10]:=12.55;
c1[10]:=-1.21;
d1[10]:=0;
tkip[10]:=642.581;
//ce
a[11]:=20305;
b[11]:=8.305;
c[11]:=0;
d[11]:=0;
tmin[11]:=1611;
tmax[11]:=2038;
tk[11]:=2038;
tkip[11]:=1419.330;
//cd
a[12]:=5908;
b[12]:=9.717;
c[12]:=-0.232;
d[12]:=-0.284;
tmin[12]:=450;
tmax[12]:=1038;
tk[12]:=594;
a1[12]:=5819;
b1[12]:=12.287;
c1[12]:=-1.257;
d1[12]:=0;
tkip[12]:=393.758;
//co
a[13]:=22210;
b[13]:=10.817;
c[13]:=0;
d[13]:=-0.223;
tmin[13]:=1000;
tmax[13]:=1772;
tk[13]:=1772;
tkip[13]:=1344.577;
//cr
a[14]:=20680;
b[14]:=14.56;
c[14]:=-1.31;
d[14]:=0;
tmin[14]:=298;
tmax[14]:=2130;
tk[14]:=2130;
tkip[14]:=1253.130;
//cs
a[15]:=4075;
b[15]:=11.38;
c[15]:=-1.45;
d[15]:=0;
tmin[15]:=280;
tmax[15]:=1000;
tk[15]:=1000;
tkip[15]:=295.308;
//cu
a[16]:=17870;
b[16]:=10.63;
c[16]:=-0.236;
d[16]:=-0.16;
tmin[16]:=298;
tmax[16]:=2870;
tk[16]:=1356;
a1[16]:=17650;
b1[16]:=13.39;
c1[16]:=-1.273;
d1[16]:=0;
tkip[16]:=1136.177;
//eu
a[17]:=8980;
b[17]:=8.16;
c[17]:=0;
d[17]:=0;
tmin[17]:=696;
tmax[17]:=900;
tk[17]:=900;
tkip[17]:=634.134;
//fe
a[18]:=21080;
b[18]:=16.89;
c[18]:=-2.14;
d[18]:=0;
tmin[18]:=900;
tmax[18]:=3000;
tk[18]:=1812;
a1[18]:=19710;
b1[18]:=13.27;
c1[18]:=-1.27;
d1[18]:=0;
tkip[18]:=1298.995;
//ga
a[19]:=14700;
b[19]:=10.07;
c[19]:=-0.5;
d[19]:=0;
tmin[19]:=303;
tmax[19]:=2676;
tk[19]:=2676;
tkip[19]:=1008.984 ;
//ge
a[20]:=20150;
b[20]:=13.28;
c[20]:=-0.91;
d[20]:=0;
tmin[20]:=298;
tmax[20]:=3103;
tk[20]:=1211;
a1[20]:=18700;
b1[20]:=12.87;
c1[20]:=-1.16;
d1[20]:=0;
tkip[20]:=1223.052;
//hf(a)
a[21]:=32000;
b[21]:=11.81;
c[21]:=-0.5;
d[21]:=0;
tmin[21]:=298;
tmax[21]:=2023;
tk[21]:=2023;
tkip[21]:=191.968;
//hf(b)
a[22]:=31630;
b[22]:=11.63;
c[22]:=-0.5;
d[22]:=0;
tmin[22]:=2023;
tmax[22]:=2503;
tk[22]:=2503;
tkip[22]:=1978.999;
//hf
a[23]:=29830;
b[23]:=9.20;
c[23]:=0;
d[23]:=0;
tmin[23]:=2503;
tmax[23]:=5470;
tk[23]:=5470;
tkip[23]:=1962.364;
//hg
a[24]:=3308;
b[24]:=10.373;
c[24]:=-0.8;
d[24]:=0;
tmin[24]:=298;
tmax[24]:=630;
tk[24]:=630;
tkip[24]:=228.340;
//i2
a[25]:=3578;
b[25]:=17.72;
c[25]:=-2.51;
d[25]:=0;
tmin[25]:=298;
tmax[25]:=458;
tk[25]:=387;
a1[25]:=30205;
b1[25]:=23.65;
c1[25]:=-5.18;
d1[25]:=0;
tkip[25]:=199.342;
//in
a[26]:=12580;
b[26]:=9.79;
c[26]:=-0.45;
d[26]:=0;
tmin[26]:=429;
tmax[26]:=2353;
tk[26]:=2353;
tkip[26]:=869.482;
//Ir
a[27]:=35070;
b[27]:=13.18;
c[27]:=-0.7;
d[27]:=0;
tmin[27]:=298;
tmax[27]:=2683;
tk[27]:=2683;
tkip[27]:=2080.315 ;
//k
a[28]:=4770;
b[28]:=11.58;
c[28]:=-1.370;
d[28]:=0;
tmin[28]:=350;
tmax[28]:=1050;
tk[28]:=1050;
tkip[28]:=337.908 ;
//la
a[29]:=22120;
b[29]:=10.39;
c[29]:=-0.33;
d[29]:=0;
tmin[29]:=298;
tmax[29]:=3730;
tk[29]:=1194;
a1[29]:=21530;
b1[29]:=9.89;
c1[29]:=-0.33;
d1[29]:=0;
tkip[29]:=1450.047;
//li
a[30]:=8415;
b[30]:=11.34;
c[30]:=-1;
d[30]:=0;
tmin[30]:=454;
tmax[30]:=1620;
tk[30]:=1620;
tkip[30]:=577.184 ;
//mg
a[31]:=7780;
b[31]:=11.41;
c[31]:=-0.855;
d[31]:=0;
tmin[31]:=298;
tmax[31]:=1363;
tk[31]:=922;
a1[31]:=7550;
b1[31]:=12.79;
c1[31]:=-1.41;
d1[31]:=0;
tkip[31]:=515.513;
//mn
a[32]:=14850;
b[32]:=17.88;
c[32]:=-2.52;
d[32]:=0;
tmin[32]:=993;
tmax[32]:=2235;
tk[32]:=1517;
a1[32]:=13900;
b1[32]:=17.27;
c1[32]:=-2.52;
d1[32]:=0;
tkip[32]:=903.753 ;
//mo
a[33]:=34700;
b[33]:=11.66;
c[33]:=-0.236;
d[33]:=-0.145;
tmin[33]:=298;
tmax[33]:=2890;
tk[33]:=2890;
tkip[33]:=2093.717 ;
//na
a[34]:=5700;
b[34]:=11.33;
c[34]:=-1.718;
d[34]:=0;
tmin[34]:=400;
tmax[34]:=1200;
tk[34]:=1200;
tkip[34]:=446.041 ;
//n2
a[35]:=6540;
b[35]:=10.7;
c[35]:=0;
d[35]:=0;
tmin[35]:=63;
tmax[35]:=77;
tk[35]:=77;
tkip[35]:=391.600 ;
//nb
a[36]:=37650;
b[36]:=8.94;
c[36]:=0.715;
d[36]:=-0.166;
tmin[36]:=298;
tmax[36]:=2741;
tk[36]:=2741;
tkip[36]:=2219.293 ;
//ni
a[37]:=22500;
b[37]:=13.6;
c[37]:=-0.96;
d[37]:=0;
tmin[37]:=298;
tmax[37]:=3005;
tk[37]:=1726;
a1[37]:=22400;
b1[37]:=16.95;
c1[37]:=-2.01;
d1[37]:=0;
tkip[37]:=1355.931 ;
//p4(yellow)
a[38]:=3530;
b[38]:=19.09;
c[38]:=-3.5;
d[38]:=0;
tmin[38]:=298;
tmax[38]:=317;
tk[38]:=317;
tkip[38]:=207.915 ;
//p4
a[39]:=2740;
b[39]:=7.84;
c[39]:=0;
d[39]:=0;
tmin[39]:=317;
tmax[39]:=553;
tk[39]:=553;
tkip[39]:=197.967 ;
//pb
a[40]:=10130;
b[40]:=11.16;
c[40]:=-0.985;
d[40]:=0;
tmin[40]:=600;
tmax[40]:=2030;
tk[40]:=2030;
tkip[40]:=705.686 ;
//pd
a[41]:=19800;
b[41]:=11.82;
c[41]:=-0.755;
d[41]:=0;
tmin[41]:=298;
tmax[41]:=3413;
tk[41]:=1825;
a1[41]:=17500;
b1[41]:=4.81;
c1[41]:=1;
d1[41]:=0;
tkip[41]:=1279.491 ;
//pr
a[42]:=17190;
b[42]:=8.1;
c[42]:=0;
d[42]:=0;
tmin[42]:=1425;
tmax[42]:=1692;
tk[42]:=1692;
tkip[42]:=1219.087 ;
//pt
a[43]:=29200;
b[43]:=13.24;
c[43]:=-0.855;
d[43]:=0;
tmin[43]:=298;
tmax[43]:=4100;
tk[43]:=2045;
a1[43]:=28500;
b1[43]:=14.3;
c1[43]:=-1.26;
d1[43]:=0;
tkip[43]:=1773.682 ;
//pu
a[44]:=17590;
b[44]:=7.9;
c[44]:=0;
d[44]:=0;
tmin[44]:=1392;
tmax[44]:=1793;
tk[44]:=1793;
tkip[44]:=1265.402 ;
//rb
a[45]:=4560;
b[45]:=12;
c[45]:=-1.45;
d[45]:=0;
tmin[45]:=312;
tmax[45]:=952;
tk[45]:=952;
tkip[45]:=317.246 ;
//re
a[46]:=40800;
b[46]:=14.2;
c[46]:=-1.16;
d[46]:=0;
tmin[46]:=298;
tmax[46]:=3000;
tk[46]:=3000;
tkip[46]:=2509.663 ;
//rh
a[47]:=29360;
b[47]:=13.5;
c[47]:=-0.88;
d[47]:=0;
tmin[47]:=298;
tmax[47]:=2239;
tk[47]:=2239;
tkip[47]:=1764.021 ;
//ru
a[48]:=33550;
b[48]:=10.76;
c[48]:=0;
d[48]:=0;
tmin[48]:=2000;
tmax[48]:=2500;
tk[48]:=2500;
tkip[48]:=2001.704 ;
//s2
a[49]:=6975;
b[49]:=16.22;
c[49]:=-1.53;
d[49]:=-1;
tmin[49]:=386;
tmax[49]:=718;
tk[49]:=718;
tkip[49]:=390.423 ;
//s(x)
a[50]:=4830;
b[50]:=23.88;
c[50]:=-5;
d[50]:=0;
tmin[50]:=386;
tmax[50]:=718;
tk[50]:=718;
tkip[50]:=272.866 ;
//sb(x)
a[51]:=11560;
b[51]:=22.4;
c[51]:=-3.52;
d[51]:=0;
tmin[51]:=298;
tmax[51]:=904;
tk[51]:=904;
tkip[51]:=622.680 ;
//sb(2)
a[52]:=11170;
b[52]:=18.54;
c[52]:=-3.02;
d[52]:=0;
tmin[52]:=904;
tmax[52]:=1908;
tk[52]:=1908;
tkip[52]:=700.412 ;
//sc(b)
a[53]:=19700;
b[53]:=13.07;
c[53]:=-1;
d[53]:=0;
tmin[53]:=1607;
tmax[53]:=1814;
tk[53]:=1814;
tkip[53]:=1232.805 ;
//se(x)
a[54]:=4990;
b[54]:=8.09;
c[54]:=0;
d[54]:=0;
tmin[54]:=493;
tmax[54]:=958;
tk[54]:=958;
tkip[54]:=354.134 ;
//si
a[55]:=20900;
b[55]:=10.84;
c[55]:=-0.565;
d[55]:=0;
tmin[55]:=1683;
tmax[55]:=2628;
tk[55]:=2628;
tkip[55]:=1387.285 ;
//sm
a[56]:=11170;
b[56]:=13.76;
c[56]:=-1.56;
d[56]:=0;
tmin[56]:=298;
tmax[56]:=1350;
tk[56]:=1350;
tkip[56]:=730.371 ;
//sn
a[57]:=15500;
b[57]:=8.23;
c[57]:=0;
d[57]:=0;
tmin[57]:=505;
tmax[57]:=2543;
tk[57]:=2543;
tkip[57]:=1089.193 ;
//sr
a[58]:=9450;
b[58]:=13.08;
c[58]:=-1.31;
d[58]:=0;
tmin[58]:=813;
tmax[58]:=1657;
tk[58]:=1042;
a1[58]:=9000;
b1[58]:=12.63;
c1[58]:=-1.31;
d1[58]:=0;
tkip[58]:=612.457 ;
//ta
a[59]:=40800;
b[59]:=10.29;
c[59]:=0;
d[59]:=0;
tmin[59]:=298;
tmax[59]:=3269;
tk[59]:=3269;
tkip[59]:=2504.494 ;
//te2
a[60]:=9175;
b[60]:=19.68;
c[60]:=-2.71;
d[60]:=0;
tmin[60]:=298;
tmax[60]:=5698;
tk[60]:=3269;
a1[60]:=7830;
b1[60]:=22.69;
c1[60]:=-4.27;
d1[60]:=0;
tkip[60]:=499.520;
//th
a[61]:=30200;
b[61]:=12.95;
c[61]:=-1;
d[61]:=0;
tmin[61]:=298;
tmax[61]:=2023;
tk[61]:=2023;
tkip[61]:=1927.784 ;

[incarnus]

Код:

//ti(b)
a[62]:=24400;
b[62]:=13.17;
c[62]:=-0.91;
d[62]:=0;
tmin[62]:=1155;
tmax[62]:=1933;
tk[62]:=1933;
tkip[62]:=1498.696;
//ti
a[63]:=23200;
b[63]:=11.74;
c[63]:=-0.66;
d[63]:=0;
tmin[63]:=1933;
tmax[63]:=3560;
tk[63]:=3560;
tkip[63]:=1482.633;
//tl
a[64]:=9300;
b[64]:=11.1;
c[64]:=-0.892;
d[64]:=0;
tmin[64]:=700;
tmax[64]:=1800;
tk[64]:=1800;
tkip[64]:=637.012 ;
//tm
a[65]:=12550;
b[65]:=9.18;
c[65]:=0;
d[65]:=0;
tmin[65]:=807;
tmax[65]:=1219;
tk[65]:=1219;
tkip[65]:=826.707;
//u
a[66]:=25580;
b[66]:=18.58;
c[66]:=-2.62;
d[66]:=0;
tmin[66]:=298;
tmax[66]:=4200;
tk[66]:=1405;
a1[66]:=24090;
b1[66]:=13.2;
c1[66]:=-1.26;
d1[66]:=0;
tkip[66]:=1588.222;
//v
a[67]:=26900;
b[67]:=10.12;
c[67]:=0.33;
d[67]:=-0.265;
tmin[67]:=298;
tmax[67]:=2160;
tk[67]:=2160;
tkip[67]:=1605.675;
//w
a[68]:=44000;
b[68]:=8.76;
c[68]:=0.5;
d[68]:=0;
tmin[68]:=298;
tmax[68]:=3680;
tk[68]:=3680;
tkip[68]:=2670.868;
//y
a[69]:=22230;
b[69]:=11.835;
c[69]:=-0.66;
d[69]:=0;
tmin[69]:=298;
tmax[69]:=3611;
tk[69]:=1795;
a1[69]:=22280;
b1[69]:=16.13;
c1[69]:=-1.97;
d1[69]:=0;
tkip[69]:=1410.795;
//zn
a[70]:=6883;
b[70]:=9.418;
c[70]:=-0.0503;
d[70]:=-0.33;
tmin[70]:=473;
tmax[70]:=1000;
tk[70]:=692;
a1[70]:=6670;
b1[70]:=12;
c1[70]:=-1.126;
d1[70]:=0;
tkip[70]:=454.775;
//zr
a[71]:=31820;
b[71]:=11.78;
c[71]:=-0.5;
d[71]:=0;
tmin[71]:=1125;
tmax[71]:=4650;
tk[71]:=2125;
a1[71]:=30300;
b1[71]:=9.38;
c1[71]:=0;
d1[71]:=0;
tkip[71]:=1972.327;

Админ: Пользуемся тегами для оформления кода!

[AND_REY]

Делай цикл Repeat Until в нём изменяй температуру от мин. до макс. с какимто шагом 0.1 к примеру. Для каждого изменения температуры считай формулу и смотри чтоб разница результата формулы и log(p) была минимально допустима. Результатом будет температура при которой разница меньше какого нибуть ебсилон (ошибки).
Ещё надо сделать внешний цикл For 1..70 чтоб вышесказанное посчитать 70 раз для разных коэфициентов, наверное взятых из файла.