Перевод с C++ на Python

#include <iostream>
#include <math.h>
using namespace std;
double F(double x)
{
  return log10(x);
}
int main()
{
  setlocale(LC_ALL, "RUSSIAN");
  int n;
  long double  a, b, h, t1=0, t2=0;
//ввод точек интеграла от а до b
nachalo:  
  cout << "a="; cin >> a;
  cout << "b="; cin >> b;
  cout << "n="; cin >> n;
  if (a < b&&a>0&&b>0) //тут надо ввести свои условия(например: тут интеграл от а до b log10(x). И условия тут будут a<b&&a>0&&b>0. Так как x>0
  {
    h = (b - a) / n;
    double *x = new double[n+1];
    double *y = new double[n+1];
    x[0] = a; y[0] = F(x[0]);
    for (int i = 1; i <= n; i++)
    {  x[i] = x[i - 1] + h;
      y[i] = F(x[i]);      
    }
    cout << "i |   x[i]   |   y[i]   \n";
    for (int i = 0; i <= n; i++)
    {
      cout << i << " |  " << x[i] << "  |  " << y[i] << endl;
    }
    for (int i = 1; i<n;){t1 += y[i]; i += 2;}
    for (int i = 2; i<n;){t2 += y[i]; i += 2;}
    cout << "S от " << a << " до " << b << " log10(x) = " << h / 3 * (y[0] + y[n] + 4 * t1 + 2 * t2) << endl;
    delete[]x;
    delete[]y;  
  }
  else
  {
    cout << "Нет решений. Введите заново\n";
    goto nachalo;
  }
  system("pause");
}

#include <iostream>
#include <math.h>
using namespace std;
double F(double x)
{
  return *****
}
int Fac(int x)
{
  if (x == 0 || x == 1){return 1;}
  else { return x*Fac(x - 1); }
}
double dF(double q, double dy, double d2y, double d3y, double d4y, double d5y)
{return dy + (2 * q - 1)*d2y / Fac(2) + (3 * q*q - 6 * q + 2)*d3y / Fac(3) + (2 * pow(q, 3) - 9 * q*q + 11 * q - 3)*d4y / Fac(4) + (5 * pow(q, 4) - 40 * pow(q, 3) + 105 * q*q - 100 * q + 24)*d5y / Fac(5);}
double d2F(double q, double d2y, double d3y, double d4y, double d5y)
{
  return d2y + (q - 1)*d3y + (6 * q*q - 18 * q + 11)*d4y / Fac(4) + (20 * pow(q, 3) - 120 * q*q + 210 * q - 100)*d5y / Fac(5);
}
int main()
{
  setlocale(LC_ALL, "rus");
  double a, b, h, q, xt; int n;
  cout << "y(x) = x^x\n";
n:  cout << "a = "; cin >> a;
  cout << "b = "; cin >> b;
  if (a > 0 && b > 0)
  {
    cout << "Число отрезков n = "; cin >> n;
    double *x = new double[n + 1]; double *y = new double[n + 1];
    double *dy = new double[n + 1];  double *d2y = new double[n + 1];
    double *d3y = new double[n + 1]; double *d4y = new double[n + 1];
    double *d5y = new double[n + 1];
    h = (b - a) / n;
    x[0] = a; y[0] = F(x[0]);
    for (int i = 1; i <= n; i++)
    {
      x[i] = x[i - 1] + h;
      y[i] = F(x[i]);
    }
 
    for (int i = 0; i < n; i++)  { dy[i] = y[i + 1] - y[i]; }
    for (int i = 0; i < n - 1; i++){ d2y[i] = dy[i + 1] - dy[i]; }
    for (int i = 0; i < n - 2; i++){ d3y[i] = d2y[i + 1] - d2y[i]; }
    for (int i = 0; i < n - 3; i++){ d4y[i] = d3y[i + 1] - d3y[i]; }
    for (int i = 0; i < n - 4; i++){ d5y[i] = d4y[i + 1] - d4y[i]; }
 
    cout << "x = "; cin >> xt;
    q = (xt - x[0]) / h;
    cout << "Ответ (y(" << xt << "))' = " << 1 / h*dF(q, dy[0], d2y[0], d3y[0], d4y[0], d5y[0]) << endl;
    cout << "Ответ (y(" << xt << "))'' = " << 1 / h*d2F(q, d2y[0], d3y[0], d4y[0], d5y[0]) << endl;
 
 
    delete[]x;  delete[]y;  delete[]dy;  delete[]d2y;  delete[]d3y;  delete[]d4y;  delete[]d5y;
  }
  else { goto n; }
  system("pause");
}