Tính diện tích tam giác trên ngôn ngữ lập trình FORTRAN 95

Code tính diện tích tam giác bằng công thức Heron (Heron's formula):

PROGRAM AREA_OF_TRIANGLE

REAL::a,b,c,p,S

PRINT*,'Nhap chieu dai canh thu nhat:'
READ*,a

PRINT*,'Nhap chieu dai canh thu hai:'
READ*,b

PRINT*,'Nhap chieu dai canh thu ba:'
READ*,c

IF (((a+b)<=c).or.((a+c)<=b).or.((b+c)<=a)) THEN
  PRINT*,'Day khong phai la tam giac'
ELSE
  p=(a+b+c)/2
  S=sqrt(p*(p-a)*(p-b)*(p-c))
PRINT*,'Dien tich tam giac la: ',S
END IF

END PROGRAM