Перевод кода Pascal в Python

@mlem · Регистрация: 16.02.2023

Студворк — интернет-сервис помощи студентам

Помогите, пожалуйста, перевести код из паскаля в пайтон

Pascal
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
Program kinel;
 type mas=array[1..2] of real;
 arr=array[1..4] of real;
 const Ca0=0.8;Cb0=0.0;Cc0=0.0;Cd0=0;
 k1:mas=(0.38,0.27);T1=580.0;
 E:mas=(75000.0,56000.0);
 R=8.314; Tn=550.0;Tk=600.0;
 hT=10.0;ttk=10.0;h=0.05;
 var
 c,f:arr;
 k,k0:mas;T,tt:real;i,j,n,jn:integer;
 F1:text;
Procedure pr(c:arr;k:mas;var f:arr);
 begin
 f[1]:=-k[1]*c[1];
 f[2]:=k[1]*c[1]-k[2]*c[2]*c[3]*c[3];
 f[3]:=2*(k[1]*c[1]-k[2]*c[2]*c[3]*c[3]);
 f[4]:=k[2]*c[2]*c[3]*c[3];
 end;
 begin assign(F1,'rel.pas');
 rewrite(F1);n:=Trunc(Round(ttk/h/10));
 for i:=1 to 2 do
 k0[i]:=k1[i]*exp(E[i]/R/T1);
 T:=Tn;jn:=1;
 repeat tt:=0.0;j:=0;
 c[1]:=Ca0;c[2]:=Cb0;c[3]:=Cc0;c[4]:=Cd0;
 for i:=1 to 2 do
 k[i]:=k0[i]*exp(-E[i]/R/T);
 writeln(F1,'Таблица ',jn,' –Расчет кинетики гомогенной химической');
 writeln(F1,' реакции при температуре T=',T:4:0,' C');
 writeln(F1,'Константы k1=',k[1]:4:2);
 writeln(F1,' k2=',k[2]:4:2);
 writeln(F1,'Время,c Ca Cb Cc Cd');
 writeln(F1,tt:5:0,c[1]:6:2,c[2]:6:2,c[3]:6:2,c[4]:6:2);
 while tt<=ttk do
 begin
 tt:=tt+h;j:=j+1;
 pr(c,k,f);
 for i:=1 to 4 do
 c[i]:=c[i]+h*f[i];
 if j=n then
 begin
 writeln(F1,tt:5:0,c[1]:6:2,c[2]:6:2,c[3]:6:2,c[4]:6:2);
 j:=0;
 end;
 end;
 writeln(F1,'Степень превращения Ха=',((Ca0-c[1])/Ca0):6:2);
 writeln(F1);
 T:=T+hT;jn:=jn+1;
 until T>Tk;
 close(f1);
 end.

Рыжий Лис · 17.02.2023, 17:28

Python
1
2
3
4
5
6
7
8
9
# Program kinel;
 
ttk = 10.0
h = 0.05
 
n = int(ttk / h / 10)  # todo округлить
 
with open('rel.pas', 'w') as f:  # rewrite?
    print('Таблица ', jn, ' –Расчет кинетики гомогенной химической', file=f)

@u235 · 18.02.2023, 19:35

Рыжий Лис, похоже, что даже на Паскале программа неправильно работает...

Не по теме:

В свое время, в ВУЗе, нам курс по химической кинетике читали, реакции первого, второго порядка, энергия активации... Было интересно, помню, что в основе лежат уравнения баланса и дифф. уравнения.

@rim41 · 20.02.2023, 01:51

Просьба не кидаться помидорами и отрицательными отзывами за читы, использовал только ради помощи ТС

Python
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
import math
 
 
def pr(c, k, f):
    f[0] = -k[0] * c[0]
    f[1] = k[0] * c[0] - k[1] * c[1] * c[2] ** 2
    f[2] = 2 * (k[0] * c[0] - k[1] * c[1] * c[2] ** 2)
    f[3] = k[1] * c[1] * c[2] ** 2
 
n = int(round(ttk / h / 10))
k0 = [k1[i] * math.exp(E[i] / R / T1) for i in range(2)]
T = Tn
jn = 1
with open("rel.pas", "w") as f1:
    while T < Tk:
        tt = 0.0
        j = 0
        c = [Ca0, Cb0, Cc0, Cd0]
        k = [k0[i] * math.exp(-E[i] / R / T) for i in range(2)]
        f1.write(f"Таблица {jn} - Расчет кинетики гомогенной химической\n")
        f1.write(f"реакции при температуре T = {T:4.0f} C\n")
        f1.write(f"Константы k1 = {k[0]:4.2f}\n")
        f1.write(f"k2 = {k[1]:4.2f}\n")
        f1.write("Время,c Ca Cb Cc Cd\n")
        f1.write(f"{tt:5.0f} {c[0]:6.2f} {c[1]:6.2f} {c[2]:6.2f} {c[3]:6.2f}\n")
        while tt <= ttk:
            tt += h
            j += 1
            pr(c, k, f)
            c = [c[i] + h * f[i] for i in range(4)]
            if j == n:
                f1.write(f"{tt:5.0f} {c[0]:6.2f} {c[1]:6.2f} {c[2]:6.2f} {c[3]:6.2f}\n")
                j = 0
        f1.write(f"Степень превращения Ха = {(Ca0 - c[0]) / Ca0:6.2f}\n\n")
        T += hT
        jn += 1

Новые блоги и статьи Все статьи Все блоги /
сукцессия микоризы: основная теория в виде двух уравнений. anaschu 11.01.2026 https:/ / rutube. ru/ video/ 7a537f578d808e67a3c6fd818a44a5c4/	WordPad для Windows 11 Jel 10.01.2026 WordPad для Windows 11 — это приложение, которое восстанавливает классический текстовый редактор WordPad в операционной системе Windows 11. После того как Microsoft исключила WordPad из. . .	Classic Notepad for Windows 11 Jel 10.01.2026 Old Classic Notepad for Windows 11 Приложение для Windows 11, позволяющее пользователям вернуть классическую версию текстового редактора «Блокнот» из Windows 10. Программа предоставляет более. . .	Почему дизайн решает? Neotwalker 09.01.2026 В современном мире, где конкуренция за внимание потребителя достигла пика, дизайн становится мощным инструментом для успеха бренда. Это не просто красивый внешний вид продукта или сайта — это. . .
Модель микоризы: классовый агентный подход 3 anaschu 06.01.2026 aa0a7f55b50dd51c5ec569d2d10c54f6/ O1rJuneU_ls https:/ / vkvideo. ru/ video-115721503_456239114	Owen Logic: О недопустимости использования связки «аналоговый ПИД» + RegKZR ФедосеевПавел 06.01.2026 Owen Logic: О недопустимости использования связки «аналоговый ПИД» + RegKZR ВВЕДЕНИЕ Введу сокращения: аналоговый ПИД — ПИД регулятор с управляющим выходом в виде числа в диапазоне от 0% до. . .	Модель микоризы: классовый агентный подход 2 anaschu 06.01.2026 репозиторий https:/ / github. com/ shumilovas/ fungi ветка по-частям. коммит Create переделка под биомассу. txt вход sc, но sm считается внутри мицелия. кстати, обьем тоже должен там считаться. . . .	Расчёт токов в цепи постоянного тока igorrr37 05.01.2026 / * Дана цепь постоянного тока с сопротивлениями и источниками (напряжения, ЭДС и тока). Найти токи и напряжения во всех элементах. Программа составляет систему уравнений по 1 и 2 законам Кирхгофа и. . .

@u235 5515 / 2868 / 571 Регистрация: 07.11.2019 Сообщений: 4,757
	18.02.2023, 19:35
	Рыжий Лис, похоже, что даже на Паскале программа неправильно работает... Не по теме: В свое время, в ВУЗе, нам курс по химической кинетике читали, реакции первого, второго порядка, энергия активации... Было интересно, помню, что в основе лежат уравнения баланса и дифф. уравнения. 0