summaryrefslogtreecommitdiff
path: root/packages/gsl/src/Numeric
diff options
context:
space:
mode:
Diffstat (limited to 'packages/gsl/src/Numeric')
-rw-r--r--packages/gsl/src/Numeric/GSL/gsl-ode.c12
1 files changed, 6 insertions, 6 deletions
diff --git a/packages/gsl/src/Numeric/GSL/gsl-ode.c b/packages/gsl/src/Numeric/GSL/gsl-ode.c
index 72c8617..4e24cb6 100644
--- a/packages/gsl/src/Numeric/GSL/gsl-ode.c
+++ b/packages/gsl/src/Numeric/GSL/gsl-ode.c
@@ -40,9 +40,9 @@ int ode(int method, int control, double h,
40 case 5 : {T = gsl_odeiv_step_rk2imp; break; } 40 case 5 : {T = gsl_odeiv_step_rk2imp; break; }
41 case 6 : {T = gsl_odeiv_step_rk4imp; break; } 41 case 6 : {T = gsl_odeiv_step_rk4imp; break; }
42 case 7 : {T = gsl_odeiv_step_bsimp; break; } 42 case 7 : {T = gsl_odeiv_step_bsimp; break; }
43 case 8 : { printf("Sorry: ODE rk1imp not available in this GSL version\n"); exit(0); } 43 case 8 : { fprintf(stderr,"Sorry: ODE rk1imp not available in this GSL version\n"); exit(0); }
44 case 9 : { printf("Sorry: ODE msadams not available in this GSL version\n"); exit(0); } 44 case 9 : { fprintf(stderr,"Sorry: ODE msadams not available in this GSL version\n"); exit(0); }
45 case 10: { printf("Sorry: ODE msbdf not available in this GSL version\n"); exit(0); } 45 case 10: { fprintf(stderr,"Sorry: ODE msbdf not available in this GSL version\n"); exit(0); }
46 default: ERROR(BAD_CODE); 46 default: ERROR(BAD_CODE);
47 } 47 }
48 48
@@ -180,11 +180,11 @@ int ode(int method, int control, double h,
180 if (status != GSL_SUCCESS) { 180 if (status != GSL_SUCCESS) {
181 int k; 181 int k;
182 printf ("error in ode, return value=%d\n", status); 182 printf ("error in ode, return value=%d\n", status);
183 printf("last successful values are:\n"); 183 fprintf(stderr,"last successful values are:\n");
184 printf("t = %.5e\n", t); 184 fprintf(stderr,"t = %.5e\n", t);
185 for (k=0; k < xin; k++) 185 for (k=0; k < xin; k++)
186 { 186 {
187 printf("y[%d] = %.5e\n", k, y[k]); 187 fprintf(stderr,"y[%d] = %.5e\n", k, y[k]);
188 } 188 }
189 break; 189 break;
190 } 190 }