Перевести программу из C++ в Python

#define _USE_MATH_DEFINES //
#include <cmath> //
#include <iostream> //
#include <fstream> //
#include <iomanip>
using namespace std; //
int main() {
const int m = 350, n = 250;
float x[m], A[m], B[m], C[m], V[m], We[m], t[n], P[n][m]; //
double h = 0.01, tau = 0.045;
float DI = 0, e = 0;
x[0] = 0.0001; //
double x0 = 0.25, sigma = 0.05;
for (int i = 0; i < m; i++) { //
x[i] = x[0] + i * h;
P[0][i] = exp(-pow((x[i] - x0), 2) / (2 * pow(sigma, 2))) / (sqrt(2 * M_PI) * sigma);
}
 
for (int j = 1; j < n; j++) { //
P[j][0] = 0;
P[j][m - 1] = 0;
}
 
for (int j = 1; j < n; j++) {
for (int i = 1; i < m - 1; i++) {
P[j][i] = 0.0001;
}
}
 
double k = 0.5;
double ksi = 0.005; //
int i = 0;
t[i] = 0;
A[i] = tau * (x[i] * ksi / h + (-ksi + 0.5 * x[i] * (1 - x[i] * (t[i] - 2 * k) / 10)) / 2) / (2 * h);
B[i] = 1 + tau * (x[i] * ksi / pow(h, 2) + (0.5 - 0.05 * x[i] * (t[i] - 2 * k) + 0.5 * x[i] * (-1 / 10 * t[i] + 1 / 5 * k)) / 2);
C[i] = tau * (x[i] * ksi / h + (ksi - 0.5 * x[i] * (1 - x[i] * (t[i] - 2 * k) / 10)) / 2) / (2 * h);
 
for (int j = 1; j < n; j++) {
V[0] = C[0] / B[0];
DI = A[0] * P[j - 1][0] + (2 - B[0]) * P[j - 1][1] + C[0] * P[j - 1][2];
We[0] = DI / B[0];
for (int i = 1; i < m - 1; i++)
{
t[j] = t[j - 1] + tau;
A[i] = tau * (x[i] * ksi / h + (-ksi + 0.5 * x[i] * (1 - x[i] * (t[j] - 2 * k) / 10)) / 2) / (2 * h);
B[i] = 1 + tau * (x[i] * ksi / pow(h, 2) + (0.5 - 0.05 * x[i] * (t[j] - 2 * k) + 0.5 * x[i] * (-1 / 10 * t[j] + 1 / 5 * k)) / 2);
C[i] = tau * (x[i] * ksi / h + (ksi - 0.5 * x[i] * (1 - x[i] * (t[j] - 2 * k) / 10)) / 2) / (2 * h);
}
for (int i = 1; i < m - 1; i++) {
e = -A[i] * V[i - 1] + B[i];
V[i] = C[i] / e;
DI = A[i] * P[j - 1][i - 1] + (2 - B[i]) * P[j - 1][i] + C[i] * P[j - 1][i + 1];
We[i] = (DI + A[i] * We[i - 1]) / e;
}
DI = A[m - 1] * P[j - 1][m - 3] + (2 - B[m - 1]) * P[j - 1][m - 2] + C[m - 1] * P[j - 1][m - 1];
P[j][m - 2] = (DI + A[m - 1] * We[m - 2]) / (B[m - 1] - A[m - 1] * V[m - 2]);
for (int i = m - 2; i > 0; i--) {
P[j][i] = V[i] * P[j][i + 1] + We[i];
}
 
}
fstream file;
file.open("out_file.txt", ios::out);
file.setf(ios::right);
file.setf(ios::fixed);
for (i = 0;i < n;i++)
{
for (int j = 0;j < m;j++)
{
cout « setprecision(6) « P[i][j] « " ";
file « setw(10) « P[i][j];
}
cout « endl;
file « endl;
}
return 0;
}