simple_plot.f95

program simple_plot
    !implicit none
    real, parameter :: pi = 4*atan(1.0)
    integer, parameter :: n = 100
    real, dimension(1:n) :: x, y
    real :: a = 0., b = 2*pi
    real :: increment
    integer :: i

    increment = (b - a) / (real(n) - 1)

    x(1) = 0.0
    do i = 2, n
        x(i) = x(i-1) + increment
    end do

    y = sin(x)

    print *, x(1:5)

    open(unit=1, file='data.txt')
    do i = 1, n
        write(1,*) x(i), y(i)
    end do

   close(1)

end program simple_plot