terminat0r
Silver Member | Редактировать | Профиль | Сообщение | Цитировать | Сообщить модератору
akaGM
Цитата:| write(*,*) ' интеграл равен = -1.0' |
хм, все таки +1
wlbmJames
функцию для интегрирования всегда описывайте как фортрановскую функцию , потому что если скажут а тепер посчитайте Hypergeometric0F1[1, Cos[Exp[BesselJ[0, x]]]]
то надо будет просто добавить модуль со спецфункциями и заменить одну строчку в программе
Код: function f(x)
implicit none
real(8), intent(in) :: x
real(8) :: f
f=cos(x)
end function f
subroutine resh(ff,a,b,h)
implicit none
real(8), intent(in) :: a,b,h
real(8), intent(out) :: ff
integer :: n,i
real(8) :: f
external f
n=int(abs(b-a)/h)
ff=h*(f(a)+f(b))*0.5d0
do i=1,n-1
ff=ff+h*f(h*i)
end do
end subroutine |
|
Всего записей: 2084 | Зарегистр. 31-03-2002 | Отправлено: 15:11 28-03-2011 | Исправлено: terminat0r, 16:14 28-03-2011 |
|
Ибо задание мне было дано в явном виде, именно в таком, как я выше написал. И собсна цель этого задания, не дстичь максимальной точности в вычислении интегралла, а показать, что я умею пользоваться подпрограммами и чтением/записью в файл. Фортран первый раз увидел 5 дней назад, когда и получил это задание. Первые 2 проги прошли нормально, а с этой запарка вышла. Директиву implicit none первый раз увидел, когда мне тут написали 
