[기계 수치해석] 7장 1번

※ 미리보기 이미지는 최대 20페이지까지만 지원합니다.

분야
등록일
페이지/형식
구매가격
적립금

레포트 > 공학계열
2011.07.26
11페이지 / docx
1,400원
42원 (구매자료 3% 적립)

다운로드 네이버 로그인

추천자료

소개글

[기계 수치해석] 7장 1번에 대한 자료입니다.

- 풀이
- C program
- 결과 및 고찰

본문내용

3번의 반복만에 원하는 오차율에 해당하는 초기값을 구할 수 있었다. 실제 문제를 풀기 위하여 위 코딩에서 도출하였던 결과값(초기치문제)으로 4계 Runge-Kutta법을 실행하여 온도의 변화와 최고점, 최고점위치에 대해 고찰해 보았다.

원하는 초기값에서의 Runge-Kutta 실행

y1 = A;
y2 = A2;
t=0;
printf("%lf \t %lf\n",t*1000,y1);
fprintf(Project,"%lf \t %lf\n",t*1000,y1);

for (i=0;;i++)
{
t = t + h;
k1_y1 = y1_prime(t,y1,y2);
k1_y2 = y2_prime(t,y1,y2);

k2_y1 = y1_prime(t+0.5*h,y1+0.5*h*k1_y1,y2+0.5*h*k1_y2);
k2_y2 = y2_prime(t+0.5*h,y1+0.5*h*k1_y1,y2+0.5*h*k1_y2);

k3_y1 = y1_prime(t+0.5*h,y1+0.5*h*k2_y1,y2+0.5*h*k2_y2);
k3_y2 = y2_prime(t+0.5*h,y1+0.5*h*k2_y1,y2+0.5*h*k2_y2);

k4_y1 = y1_prime(t+h , y1+h*k3_y1, y2+h*k3_y2);
k4_y2 = y2_prime(t+h , y1+h*k3_y1, y2+h*k3_y2);

y1 = y1 + h/6*(k1_y1 + 2*k2_y1 + 2*k3_y1 + k4_y1);
y2 = y2 + h/6*(k1_y2 + 2*k2_y2 + 2*k3_y2 + k4_y2);

printf("%lf \t %lf\n",t*1000,y1);
fprintf(Project,"%lf \t %lf\n",t*1000,y1);

if (t>=0.0025)
break;
}

}

double y1_prime(double t, double y1,double y2)
{
return y2;
}

double y2_prime(double t, double y1,double y2)
{
return -x;
}

#수치해석 #장 #기계 #수치 #해석