Как я могу исправить эту проблему об этом массиве, которые меняются? - PullRequest
Новая
       15

Как я могу исправить эту проблему об этом массиве, которые меняются?

0 голосов
/ 19 мая 2019

Я хотел проверить метод Тейлора, и у меня есть некоторые проблемы.Я не знаю почему, но когда я пытаюсь напечатать первую строку x, программа выдает: 3.30448736E + 17.Но я не понимаю, как это случилось.

Спасибо за помощь 


    implicit none
    real,parameter :: x1=1,y1=2,h=0.5
    real,dimension(1:20) :: x,y
    integer :: n

    n=1
    x(1)= x1
    y(1)= y1

    do while(n<21)
            y(n+1)= y(n)+h*x(n)*y(n)+0.5*y(n)*(1+x(n)**2)*h**2                   
            x(n+1)=x(n)+h
        n=n+1
    end do

    print*,x(1)
...