수치해석 6장 레포트

목차

없음

본문내용

program HW_1

implicit none

real :: x2, x3, fx2, fx3, a, b, sp, tf
x2 = 0.3 !! 값 지정
x3 = 0.5
fx2 = 0.404957
fx3 = 0.824360
!!먼저 1차로 spline 한다.
!!이때 ,0.4가 0.3과 0.5 사이에 있기 때문에 0.3과0.5 사이의 구간만 spline 하면 된다.
a=(fx3-fx2)/(x3-x2)
b=fx2
!!(x2,f(x2))와 (x3,f(x3))을 지나는 곡선의 기울기와 y절편을 구한다. 결론적으로 직선을 구하기 위함
sp=a*(0.4-x2)+b
tf=0.4*exp(0.4)
!!산출된 1차 spline에 원하는 x=0.4를 대입한다.
!!추정치를 도출한다. 원래 식 = x*exp(x)도 계산한다.

<중 략>

do n=1,3
a(n)= (s(n+1)-s(n))/(6*d(n))
b(n)= s(n)/2
c(n)= (fx(n+1)-fx(n))/d(n) - d(n)*(2*S(n)+S(n+1))/6
end do
!!근사식과 근사값 쓰기 위한 계수 행렬 값 입력
i=2.05
!!입력값 2.05로 입력
do n=1,3
if(x(n) <= i)then
if(i <= x(n+1))then
r= a(n)*(i-x(n))**3 + b(n)*(i-x(n))**2 + c(n)*(i-x(n)) + fx(n)
exit
end if
end if
end do

print `( A1, i1, A1, f5.2 , A2, f10.5)` ,"S",n,"(", i,")=", r
!!결과 값 산출 및 출력
i=2.15
!!입력값 2.15로 입력
do n=1,3
if(x(n) <= i)then
if(i <= x(n+1))then
r= a(n)*(i-x(n))**3 + b(n)*(i-x(n))**2 + c(n)*(i-x(n)) + fx(n)

참고 자료

없음