|
4th order RK method
#include <iostream.h>
#include <math.h>
double g=9.81, m=90, Cd=0.225;
double dvdt(double v, double t)
{
return g-Cd*v*v/m;
}
double dsdt(double v, double t)
{
return g*t-Cd*v*v*t/m;
}
void main()
{
double t0, h;
double s,v;
t0=0;
h=0.25;
v=0;s=0;
cout << "Use 4
|
|
|
Method\n");
printf(" h = 0.25 \t\t error\t\t h = 0.25\t\t error\n");
for(i=0 ; i<=4 ; i++)
{
t = 0.25*i;
printf("w(%1.1f) = %e\t%e\tw(%1.1f) = %f\t%e\n",t ,w[i] ,absol(y(t)-w[i]) ,t ,v[i], absol(y(t)-v[i]));
}
}
double f(double t, double w)
{
return 5*exp(5*t)*(w-t)*(w-t)+1;
}
double absol(double
|
|
|
오일러 방식과 비교해서 정상적으로 출력되어야하는 값과 좀 더 일치해지는 것을 확인할 수 있었습니다. 지금껏 배운 선형 방정식을 풀기위한 방법들을 사용해서 앞으로 전공으로 진행할 프로젝트에서도 같이 사용한다면 좋을 것 같다는 생
|
|
|
method of this study, I will study this theme according to chronological order. This research will refer to primary sources first, and it will also consult secondary sources related to the theme. Regarding the scope of my research, I will examine mainly various great monks who contributed to monasti
|
|
메뉴펼치기
