!f(x)=(e^(-x*x)-e^(-x*x/4)+(3x*x/4)(1-5x*x/8))(128/(21x^6))
 !f5 = (exp(-x**2)-exp(-x**2/4)+3/4*x**2*(1-5*x**2/8))(128.0/(21*x**6))
program V5
real*8 :: x, y, y_m, a, b
integer i, n
open(10,file = 'data.dat')   ! открываем
open(30,file="result.dat")
read(10,*) a, b, n
dx = (a+b)/(n + 0.d0)
x = a
do i = 1, n
    y = f5(x)
    y_m = macloren_f5(x)
    x = x + dx
    write(30,*) x, y, y_m, y_m-y
    write(*,*) i
enddo
contains
    function f5(x)
        real*8 :: f5, x
        f5 =(exp(-(x**2)) - exp(-0.d25*(x**2)) + 0.d75*x**2*(1-5.d0/8.d0*x**2))*(128.d0/(21.d0*x**6))
    end function f5

    function macloren_f5(x)
        real*8 :: macloren_f5, x
        macloren_f5 = -1.d0+85.d0/336.d0*x**2 - 341.d0/6720.d0*x**4 + 5461.d0/4515840.0*x**8
    end function macloren_f5
end program V5