program array_func
real, dimension(0:50) :: x
real, dimension(0:50) :: y
integer :: i=0, n=50, max_i, min_i
real :: h, a=0, b=1, m, d, plus, minus, e, k
h=(b-a)/n
do while (i<=n)    !задание массива x
    x(i) = a+i*h
    i = i+1
end do             !задание массива x
i = 0              !задание массива y
do while (i<=n)
    y(i) = FUNC(x(i))
    i = i+1
end do             !задание массива y
i = 1              !поиск максимального элемента
m = y(0)
max_i = 0
do while (i<=n)
    if (m < y(i)) then
        m = y(i)
        max_i = i
    end if
    i = i+1
end do
print*, "maximalniy element massiva", m
print*, "index maximalnogo elementa", max_i      !поиск максимального элемента
i = 1                                            !поиск минимального элемента
m = y(0)
min_i = 0
do while (i<=n)
    if (m > y(i)) then
        m = y(i)
        min_i = i
    end if
    i = i+1
end do
print*, "minimalniy element massiva", m
print*, "index minimalnogo elementa", min_i      !поиск минимального элемента
i = 0                                            !поиск значений
m = 0
do while (i<=n)
    m = m + y(i)
    i = i+1
end do
e = m/(n+1)
print*, "srednee znachenie", e
m = 0
i = 0
do while (i<=n)
    m = m + (y(i))**2
    i = i+1
end do
d = m/(n+1)
print*, "sredniy quadrat", d
d = sqrt(d)
print*, "srednequadratichnoe znachenie", d       !поиск значений
i = 0                                             !относитльное количество+и-
d = 0
m = 0
do while (i<=n)
    if (y(i) < 0) then
        d = d+1
    end if
    if (y(i) > 0) then
        m = m+1
    end if
    i = i+1
end do
plus = m/(n+1)
minus = d/(n+1)
print*, "otnositelnoe chislo pologitelnih", plus
print*, "otnositelnoe chislo otrizatelnih", minus !относитльное количество+и-
m = 0                                             !поиск среднеквадратичного отклонения
i = 0
do while (i<=n)
    m = m + (y(i)-e)**2
    i = i+1
end do
k = sqrt(m/(n+1))
print*, "srednequadratichnoe otklonenie", k       !поиск среднеквадратичного отклонения
end program

real function FUNC(x)               !задание функции
real :: x
FUNC = (1-x)**2-tan(pi*x**2/4)**2
return
end