How to get double precision pi constant in Fortran

<code>

program testpi
implicit none
double precision::mypi
mypi=acos(dble(-1.0));
write(*,*) mypi
write(*,*) sin(mypi)
end program testpi

</code>

This outputs:

<code>

3.14159265358979
1.224646799147353E-016

</code>

Note that the sin(pi) is close to the value given in Mathworks:

http://www.mathworks.com/access/helpdesk/help/techdoc/ref/pi.html

Advertisements

Leave a Reply

Fill in your details below or click an icon to log in:

WordPress.com Logo

You are commenting using your WordPress.com account. Log Out / Change )

Twitter picture

You are commenting using your Twitter account. Log Out / Change )

Facebook photo

You are commenting using your Facebook account. Log Out / Change )

Google+ photo

You are commenting using your Google+ account. Log Out / Change )

Connecting to %s

%d bloggers like this: