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

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

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

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

Приложение I. Программа для выполнения расчета нормального потенциала цепи без переноса Е*

Программа для выполнения расчета нормального потенциала цепи без переноса Е^*

10 rem “Линейная регрессия”

100 cls

1000 dim x(100), y(100)

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

1101 if N<2 goto 7040

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

2200 for i=1 to N

2300 s2=s2+x(i)

2400 s3=s3+y(i)

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

2600 s5=s5+y(i)×x(i)

2610 s6=s6+y(i)×x(i)

2620 s7=s2/N

2630 s8=s3/N

2700 next i

3000 d1=s1×s4-s2×s2

3100 d2=s3×s4-s5×s2

3200 d3=s1×s5-s2×s3

4000 a=d2/d1

4100 b=d3/d1

5000 print “Уравнение регрессии”

5100 print “ y = “;: print using “####.####”;a;: print”+ “;: print using “####.####”; b;:print “ x ”

6000 s=0

6100 for i=1 to N

6200 s=s+(y(i)-(a×b×x(i)))^2

6300 next i

6500 print “Сумма квадратов отклонений =“;: print using “####.####”;s

6550 z1=s5/s1-s7×s8

6560 z2=sqr(s4/s1-s7×s7)

6570 z3=sqr(s6/s1-s8^2)

6580 r=z1/(z2×z3)

6581 if r>1 then r=2-r

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

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

6660 a$=input$(1)

6670 if a$=“Y” or a$=“y” goto 6700

6680 if a$=“n” or a$=“N” then 9999 else 6650

6700 print

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

6800 for i=1 to N

6850 print using “####.####”; x(i),y(i),a+b×x(i)

6900 next i

7040 goto 9999

9999 end