본문내용
i] / A[i][i];
for(l=i; l < NUM; l++){
A[j][l] = A[j][l] - ( x * A[i][l] ) ;
if(A[j][l] == -0)
A[j][l] = 0;
}
}
}
solve_cj();
solve_bj();
solve_dj();
}
/* c 값을 구한다. */
void solve_cj(){
int i, j;
double sum;
for(i=NUM-2; i>=0; i--){
sum = 0;
for(j=i+1; j
sum += A[i][j] * c[j+1];
}
c[i+1] = (A[i][ NUM - 1 ] - sum) / A[i][i];
}}
/* b 값을 구한다. */
void solve_bj(){
int j;
for(j=0; j
b[j] = 1 / h(j) * ( a[j+1] - a[j] ) - h(j) / 3 * ( c[j+1] + 2 * c[j]);}
/* d 값을 구한다. */
void solve_dj(){
int j;
for(j=0; j
d[j] = 1 / ( 3 * h(j) ) * ( c[j+1] - c[j] );}
void pivot_A(int x, int y){
int i;
double tmp;
for(i=x; i
tmp = A[x][i];
A[x][i] = A[y][i];
A[y][i] = tmp;
}}
int spline_range(double in){
int i;
for(i=0; i
if(x[i] > in)
return i-1;
return -1;}
for(l=i; l < NUM; l++){
A[j][l] = A[j][l] - ( x * A[i][l] ) ;
if(A[j][l] == -0)
A[j][l] = 0;
}
}
}
solve_cj();
solve_bj();
solve_dj();
}
/* c 값을 구한다. */
void solve_cj(){
int i, j;
double sum;
for(i=NUM-2; i>=0; i--){
sum = 0;
for(j=i+1; j
}
c[i+1] = (A[i][ NUM - 1 ] - sum) / A[i][i];
}}
/* b 값을 구한다. */
void solve_bj(){
int j;
for(j=0; j
/* d 값을 구한다. */
void solve_dj(){
int j;
for(j=0; j
void pivot_A(int x, int y){
int i;
double tmp;
for(i=x; i
A[x][i] = A[y][i];
A[y][i] = tmp;
}}
int spline_range(double in){
int i;
for(i=0; i
return i-1;
return -1;}
소개글