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 |
|