Научная электронная библиотека

Монографии, изданные в издательстве Российской Академии Естествознания

Взаимодействия в растворах электролитов: моделирование сольватационных процессов, равновесий в растворах полиэлектролитов и математическое прогнозирование свойств химических систем

Танганов Б. Б.,

Приложение II. Программа расчета параболической функции y = a₀xⁿ + a₁x^n-1 + ... + a_n для определения стандартных потенциалов Е⁰ (стандартных ЭДС) электрохимических цепей

Программа расчета параболической функции y = a₀×xⁿ + a₁×x^n-1 + ... + a_n для определения стандартных  потенциалов Е⁰ (стандартных ЭДС)  электрохимических цепей

1 rem Расчет параболы

100 cls

1100 input “Число пар значений “; N

1110 if N<4 goto 7500

1120 dim x(N)

1130 dim y(N)

1200 for i=1 to N

1300 print “x(“;i;”),y(“;i;”) “;

1400 input x(i), y(i)

1700 next i

2000 s1=N

2100 s2=0;s3=0;s4=0;s5=0;s6=0;s7=0;s8=0;s9=0

2200 for i=1 to N

2250 s2=s2+y(i)

2300 s3=s3+x(i)

2400 s4=s4+x(i)×x(i)

2500 s5=s5+(x(i)×x(i)×x(i))

2600 s6=s6+(x(i)×x(i)×x(i)×x(i))

2700 s7=s7+x(i)×y(i)

2800 s8=s8+(x(i)×x(i))×y(i)

2820 o(i)=(y(i)-s2/s1)^2

2830 s9=s9+o(i)

2900 next i

3100 L1=(s4×s4-s3×s5)

3200 L2=(s3×s4-s1×s5)

3300 L3=(s4×s7-s2×s5)

3500 L4=s4×s5-s3×s6

3600 L5=s4^2-s1×s6

3700 L6=s4×s8-s2×s6

3800 c=(L3×L4-L1×L6)/(L2×L4-L1×L5)

3810 b=L3/L1-(L2/L1)×c

3811 a=s2/s4-(s3/s4)×b-(s1/s4)×c

3812 m1=o:m3=0

3820 for i=1 to N

3850 s10(i)=(a×x(i)^2+b×x(i)+c)

3860 m1=m1+s10(i)

3870 m2= m1/s1

3875 o1(i)=(s10(i)-m2)^2

3880 m3=m3+o1(i)

3881 next i

3890 R1=m3/s9

3891 R=sqr(R1)

3892 if R<1 goto 3900

3893 R=2-R

3900 cls:print “Уравнение параболы ”

4000 print “Y=“ using “####.#####”;a;:print”×x^2+” using “####.#####”; b;:print”x+”using “####.#####”;c

4100 print “Коэффициент корреляции R =“ using “#.####”;R

4150 print “Хотите сравнить исходные значения y с рассчитанными значениями Y (Y/N)”;

4200 a$=input$(1)

4250 if a$=“N” goto 7500

4251 if a$=“n” goto 7500

4253 if a$=“Y” or a$=“y” then 4400 else 4150

4400 print “ x  y(измер)  Y(теор)  “

4450 for i=1 to N

4500 print using “####.#####”;x(i),y(i),(a×x(i)^2 + b×x(i)+c)

4550 next i

7500 end