vecfor_R8P-doctest-13.f90 Source File

program volatile_doctest
use vecfor_R8P
 type(vector_R8P) :: pt(0:3)

 pt(0) = 5.3 * ez_R8P
 pt(1) = ex_R8P
 pt(2) = ey_R8P
 pt(3) = ex_R8P - ey_R8P
 pt(0) = distance_vectorial_to_plane_R8P(pt(0), pt1=pt(1), pt2=pt(2), pt3=pt(3))
 print "(3(F3.1,1X))", pt(0)%x, pt(0)%y, pt(0)%z
endprogram volatile_doctest