[수치해석]수치해석

소개글

오일러 공식, 개선된 오일러 공식, runge-kutta애 대한 수치해석
c code와 결과화면

목차

1. Euler Method.
(1)소스코드
(2)화면캡처
2. Improved Euler Method
(1)소스코드
(2)화면캡처
3. Runge-Kutta Method
(1)소스코드
(2)화면캡처
4.exact solution
(1)소스코드
(2)화면캡처

본문내용

1. Euler Method.
(1) source code
c234567

program euler
implicit double precision(a-h,o-z)
external userf
parameter (max=200)
dimension y(0:max)
open(15,file=`euler`,status=`unknown`)


write(*,*)`구간수를 지정해주세요`
read(*,*) n
write(*,*) `구간을 지정해주세요`
read(*,*) a,b
h=(b-a)/n
x0=2
y(0)=1
x=x0

call s_eulr(y,x,h,n,userf)
write(*,*)

do i=0,n
write(*,*) h*i+x0,y(i)
write(15,`(1x,f3.1,4x,f10.9)`) h*i+x0,y(i)
enddo
close(15)

write(*,*)`계산값이 euler.dat 파일에 저장.`
stop
end

subroutine s_eulr(y,x,h,n,func)
implicit double precision (a-h,o-z)
dimension y(0:n)
do i=0,n-1
y(i+1)=y(i)+h*func(x,y(i))
x=x+h
enddo
return

참고 자료

없음