1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
|
program arrays_60
implicit none
! Declare variables
integer, dimension(5) :: array
real :: scalar_value
integer :: i
real, dimension(5) :: x, y, z
array = (/ 1, 0, 3, 0, 5 /)
x = (/ 1.0, 2.0, 3.0, 4.0, 5.0 /)
y = (/ 2.0, 4.0, 6.0, 8.0, 10.0 /)
scalar_value = 2.5
print *, "Result:", square(scalar_value)
z = array_function(x, y)
print *, "Result (z):", z
if( any(abs(z - 2.41421366) > 1e-6) ) error stop
where (array == 0)
array = 10
end where
do i = 1, size(array)
print *, "array(", i, ") squared =", square(real(array(i)))
end do
contains
real function square(a)
real, intent(in) :: a
square = a * a
end function
real function array_function(a, b)
real, intent(in) :: a(:), b(:)
real :: temp(size(a))
array_function = 0.0
temp = a + sqrt(b)
array_function = temp(1)
end function
end program arrays_60
|
