Фортран.Функция с промежутком и шагом расчёта.

[Семён Т.]

Дана функция
11-7* √lx-1l ,если x<-3
y(x)= x,если -3≤x<3
7√x+1+12,если 3≤x,и дан промежуток и шаг расчёта: x=-5(0,5)5


program lab2
real x,y
write (*,'('' введите x'')')
read (*,*)x
if (x.LT.-3) then
y=11-7x SQRT (ABS (x-1))
else
if (-3.LE.x. and. x.LT.3.) then
y=x
else
y=7* SQRT (x+1+12)
end if
end if
write (*,'(''x='', f7.2, ''y='',f7.2)') x,y
end


Без шага расчёта и промежутка я написал,но в программу нужно вввести промежуток и шаг расчёта,не могу понять как это сделать.

[TinMan]

Сообщение от Семён Т.

Без шага расчёта и промежутка я написал,но в программу нужно вввести промежуток и шаг расчёта,не могу понять как это сделать.

Цикл.
Вычисляй значения x в цикле, затем вычисляй функцию и выводи на печать.

Если правильно помню, то что-то типа такого (в f77)

      x1 = -5.0
      h = 0.5
      x = x1
      DO

C тут вычисляешь значение функции и печатаешь его

      x = x + 0.5;
      UNTIL (x .gt. x1 + h*0.1)

[Семён Т.]

Вот сделал:

program lab3
real x,y,xn,xk,h.
write (*,'(''Введите xn,xk,h'')')
read (*,*) xn,xk,h.
x=xn
10 if (x.LT.-3.)then
y=11-7*(SQRT(ABS(x-1)))
else
if (-3.LE.x.and.x.LT.3.) then
else
y=x
y=7*SQRT(x+1+12)
end if
end if
write (*,''(x='',f7.3, ''y='',f7.3)')x,y
if (x.LT.xk) then
x=x+h
go to 10
end if
end

Теперь сделать тоже самое через массивы x-ов и y-ов и через цикл DO.

[TinMan]

Сообщение от Семён Т.

Теперь сделать тоже самое через массивы x-ов и y-ов и через цикл DO.

Читать не умеем, только пишем?..

[Семён Т.]

Я не совсем понял твоё объяснение,можешь целиком написать?У меня пока не получается.

[Rin]

Читать не умеем, только пишем?

Думаю TinMan об этом(т.к. все реализовано через цикл DO):

x1 = -5.0
h = 0.5
x = x1
DO

C тут вычисляешь значение функции и печатаешь его

x = x + 0.5;
UNTIL (x .gt. x1 + h*0.1)

Фортран не знаю, прошёлся по мануалам и думаю, что так:

program lab3
real x(21),y(21),xn,xk,h,i.
write (*,'(''Введите xn,xk,h'')')
read (*,*) xn,xk,h.
i=1
DATA x(i)=xn
DO
 if (x(i).LT.-3.)then
DATA y(i)=11-7*(SQRT(ABS(x(i)-1)))
else
if (-3.LE.x(i).and.x(i).LT.3.) then
else
DATA y(i)=x(i)
DATA y(i)=7*SQRT(x(i)+1+12)
end if
end if
write (*,''(x='',f7.3, ''y='',f7.3)')x(i),y(i)
if (i.LT.21) then
DATA x(i+1)=x(i)+h
end if
i=i+1
 UNTIL (i.GE.22)
end

[TinMan]

Сообщение от Rin

Думаю TinMan об этом(т.к. все реализовано через цикл DO):

Фортран не знаю, прошёлся по мануалам и думаю, что так:

Rin, спасибо, +1.
Я писал на Фортране, но жутко давно. Помню, что там нельзя начинать с самого начала строки. Если самый первый символ 'c' - то вся строка комментарий (у меня есть один в коде выше). Код топикстартера не взят в теги, поэтому сдвинут влево.

...
добавлено: фу ты блин, опять не могу поблагодарить человека... нет, это точно ненормально..

[Vago]

repeat ... until в Фортране?.. Ню - ню...
...

Сообщение от Семён Т.

Дана функция
...
7√x+1+12,если 3≤x
...

1+12 будет 13. Что под корнем-то, а что за?

Какой FORTRAN хоть? 66, 77, 90, 2003?

[Семён Т.]

Фортран 90. Под корнем x+1+12,через таблицу символов поставил криво,спешил.

[Vago]

Сообщение от Семён Т.

Фортран 90.

И что, можно использовать конструкции Fortran 90, не ограничивая себя 77-м? Типа free-form набивки, динамического выделения памяти, DO - EXIT - ENDDO и т.п?.. Или вы пользуетесь транслятором с 90-го, а программируете в манере 77-го?

Сообщение от Семён Т.

Под корнем x+1+12,через таблицу символов поставил криво,спешил.

Ничего не понял. Почему тогда не x + 13 ?!

77-й:

      program lab3

         real     x(200), y(200), xn, xk, h
         integer  i

         write ( *, '(" Введите xn, xk, h")' )
         read *, xn, xk, h

         i = 1
         x(i) = xn
c        repeat
   10       if ( x(i) .LT. -3. ) then
               y(i) = 11. - 7. * (SQRT( ABS( x(i)-1. )) )
            else if ( x(i) .LT. 3.) then
               y(i) = x(i)
            else
               y(i) = 7. * SQRT( x(i) + 1. + 12. )
            end if
            write (*, 900 ) i, x(i), i, y(i)
            x(i+1) = x(i)+h      ! Бесстрашно накапливаем погрешность
            i = i+1
            if ( x(i) .LE. xk ) goto 10
c        until ( x .GT. xk )

  900    format ( ' X(', i2, ') = ', f7.1, '   Y(', i2, ') = ', f9.3 )

      end