Перевести код с Паскаля на Фортран

uses crt;
const
  a = 0;{нижний предел}
  b = 1;{верхний предел}
  e=0.0001;{точность}
function F(x:real):real;
begin
F:=exp(2*x)-exp(x)-ln(abs(x)+1)*cos(x);
end;
 
var n,i:integer;
    h,k,s1,s2: real;
begin
clrscr;
write('Количество разбиений n=');
readln(n);//начальное число
s1:=0;//итоговый интеграл
h:=(b-a)/n;//шаг
s2:=h*F(a);//текущий интеграл
repeat
s1:=s2;//запомним текущий
h:=(b-a)/n;
k:=a;//левый край
for i:=1 to n-1 do
 begin
  k:=k+h;//наращиваем
  if i mod 2 = 0 then s2:=s2+2*f(k)//считаем новый интеграл
  else s2:=s2+4*f(k);
 end;
s2:=(s2+f(a)+f(b))*h/3;//конечное значение
n:=n*2;//если точность не достигнута удваиваем количество разбиений
until abs(s1-s2)<=e;
write('Интеграл=',s2:0:6);
readln
end.

if i mod 2 = 0 then

if i mod 2 = 0 then

program semistrovaia4;
 
const
Nmass = 100;
 
var
X: array [1..Nmass] of integer;
sum, N, i, c1, c2: integer;
max: real;
 
begin
Writeln('Введите число N');
readln(N);
for i := 1 to N do
begin
write('X[', i, ']=');
readln(X[i]);
end;
max := x[1];
c1 := 0;
c2 := 1;
for i := 1 to N do
begin
if (X[i] > max) then
begin
max := x[i];
c2 := i;
end;
if((X [i] <> max) and (x[i] mod 3 = 0) and (c1 < c2)) then
begin
c1 := i;
end;
if(X [i] < 0) then
begin
sum := sum + x[i];
end;
end;
if(c1 > c2) then
x[c1] := abs(sum)
else
writeln('Нет такого элемента , массив не изменный');
writeln(' Замена числа кратного трем,на сумму модуля отрицательных чисел стояшего после мах работает:');
for i := 1 to N do
begin
write('X[',i,']=',X[i],' ');
end;
end.

uses crt;
var a,b,x,y,z,t:real;
begin
writeln('Введите значения a,b,x');
read(a,b,x);
if a>b then y:=a*x*x*x+b*b+x*x*x*x*x
else 
 begin
  if x<>0 then 
   begin
    if cos(a/2)=0 then y:=-1
    else y:=(sin(a/2)/cos(a/2))*b*x 
   end
  else y:=-1;
 end;
if y>0 then write('y>0')
else if y=0 then write('a=',a:0:2,' b=',b:0:2,' x=',x:0:2)
else
 begin
  t:=exp(2*x)+3*a*a*a*b; 
  if t=0 then z:=0
  else if t>0 then  z:=exp(ln(t)/5)
  else z:=-exp(ln(-t)/5);
  write('z=',z:0:2);
 end;
end.

uses crt;
const a=2.8;
      b=-0.3;
      c=4;
var xn,xk,x,h,y:real;
begin
xn:=0.8;
xk:=2.8;
h:=0.2;
x:=xn;
while x<xk+h/2 do
 begin
  if x<1.4 then y:=(a+b*x)/sqrt(x*x+4)
  else if x=1.4 then y:=a*x*x+b*x+c
  else y:=a/x+exp(a*b);
  writeln('x=',x:4:1,'  y=',y:5:2);
  x:=x+h;
 end;
end.

uses crt;
var n,i,k:integer;
begin
repeat
write('Введите натуральное число >=3 n=');
read(n);
until n>2;
i:=3;
k:=0;
while i<=n do
 begin
  k:=k+1;
  write(i,' ');
  if k mod 10=0 then writeln;
  i:=i+3;
 end;
end.

uses crt;
const nmax=15;
var a:array[1..nmax] of real;
    n,i,imn,k:integer;
    mn,mx,s:real;
begin
repeat
write('Размер массива до ',nmax,' n=');
read(n);
until n in [1..nmax];
writeln('Введите элементы массива:');
for i:=1 to n do
 begin
  write('a[',i,']=');
  read(a[i]);
 end;
clrscr;
writeln('Исходный массив:');
mn:=a[1];
mx:=a[1];
imn:=1;
for i:=1 to n do
 begin
  write(a[i]:0:2,' ');
  if a[i]<mn then
   begin
    mn:=a[i];
    imn:=i;
   end;
  if a[i]>mx then mx:=a[i];
 end;
writeln;
writeln('Минимальный элемент=',mn:0:2,' его индекс=',imn);
writeln('Максимальный элемент=',mx:0:2);
repeat
write('Сколько последних элементов включить в сумму от 1 до ',n,' k=');
read(k);
until k in [1..n];
s:=0;
for i:=n downto n-k+1 do
s:=s+a[i];
writeln('сумма последних ',k,' элементов=',s:0:2);
a[imn]:=s;
writeln('замена минимального суммой:');
for i:=1 to n do
write(a[i]:0:2, ' ');
writeln;
if mn>=0 then
 begin
  writeln('Минимальный элемент >=0');
  for i:=1 to n do
  a[i]:=a[i]*mn*mn;
 end
else
 begin
  writeln('Минимальный элемент <0');
  for i:=1 to n do
  a[i]:=a[i]*mx*mx;
 end;
for i:=1 to n do
write(a[i]:0:2,' ');
end.