nSquareRoots.f95

program nSquareRoots
    parameter (n = 100)
    real, dimension (n) :: squares, squareRoots
    squares = (/ (i, i = 1, n) /)
    squareRoots = (/ (customSqrt(squares(i), 0.001), i = 1, n) /)
    Print *, "Squares:"
    Print *, squares
    Print *, "Square roots:"
    Print *, squareRoots
end program nSquareRoots

real function customSqrt(a, e)
    x = 0.5 * a
    y = (x*x + a)/(2*x)
    do while (abs(x - y) > e)
        x = y
        y = (x*x + a)/(2*x)
    end do
    customSqrt = y
    return
end function customSqrt