인하대학교/기계공학/수치해석/주현철교수님/Homework2

소개

Q ) 1 -1


// Use The bisection method

#include
#include

// y = ln(x^4) - 0.7

long double func ( long double x ) // 주어진 함수를 이용하기 위해 따로 정의하여 사용
{
return( log(pow(x,4) ) - 0.7 ) ; // 주어진 함수 f = ln(x^4) - 0.7

}

int main()
{
FILE *fp; // 출력을 파일로 저장하기 위하여 포인터 사용
long double xl = 0.5 ; // 초기값 중에 작은 것
long double xu = 2.0 ; // 초기값 중에 큰 것
int iter =0 ;
long double n ;
long double xr ;
long double xr_pre = 0.0 ;
long double ea , et ;
long double x_true = 1.191246217 ; // 방정식의 실제근

fp = fopen ( "data.txt" , "w" ); // 텍스트 파일로 저장
fprintf( fp , " iter\t xr\t\t ea\t\t et\t\n " ) ;

printf( " iter\t xr\t\t ea\t\t et\t\n " ) ;
printf("━━━━━━━━━━━━━━━━━━━━━━━━━\n");


while( iter < 30) // 반복구문 사용
{
xr = (xl + xu ) / 2 ; // bisection method에서의 xr
ea = ( (xr - xr_pre) / xr ) * 100; // 근사오차

컴파일 실행환경

Visual studio2010

참고자료

· 없음