怎么用FORTRAN写：已知三角形三点坐标,求三角形面积?

write(*,*) ' 请输入三角形三个顶点的坐标 x1,y1,x2,y2,x3,y3：'
read(*,*) x1,y1,x2,y2,x3,y3
a=sqrt((x2-x1)^2+(y2-y1)^2)
b=sqrt((x3-x1)^2+(y3-y1)^2)
c=sqrt((x3-x2)^2+(y3-y2)^2)
L=(a+b+c)/2
s=sqrt(L*(L-a)*(L-b)*(L-c))
write(*,*) ' 该三角形的面积为：',s
end