fortran cosθをベクトリから求めるプログラム

質問

ベクトルからcosθを求めるペログラムを作ってみたのですが、実際計算した値と実行結果の値が一致しないのですが。プログラム上に問題があるのでしょうか？教えて下さい。よろしくお願いします。
module subprogs
implicit none
contains
function vec_cos(a,b) result(vcos)
real(8), intent(in) :: a(:), b(:)
real(8) ab, vcos
if (size(a) /= size(b)) stop ' er : size(a) /= size(b) '
ab = dot_product(a,a) * dot_product(b,b)
if (ab == 0.0d0) then
vcos = 0.0d0
else
vcos = dot_product(a,b)/sqrt(ab)
endif
end function vec_cos
end module subprogs

program main
use subprogs
implicit none
real(8) :: x(1:2) = (/1.2d0, 3.4d0/), y(1:2) = (/5.6d0, 7.8d0/)
write(*,*) 'cos = ', vec_cos(x,y)
end program main

実行結果　cos = 0.9601163787292428

回答No.1

Tacosan · ベストアンサー

その「実際計算した値」というのはどのように計算したもので値はいくつになったのでしょうか?

お礼

cos = ( 1.2*5.6+3.4*7.8)/sqrt((1.2*1.2+3.4*3.4)*(5.6*5.6+7.8*7.8))=0.9531....
再度計算し直したらちゃんと実行結果の値と同じになりました。こちらのケアレスミスでした。ありがとうございました。