使用c語(yǔ)言解常微分方程-c-ode-資料下載頁(yè)
2025-07-22 00:08本頁(yè)面
【正文】 0] / b[0]。 for (i = 0。 i M 2。 i++) { d[i + 1] = b[i + 1] a[i] * u[i]。 if (0 == d[i + 1]) { printf(error!!!\n)。 return 0。 } u[i + 1] = c[i + 1] / d[i + 1]。 } d[M 1] =b[M 1]a[M 2] * u[M 2]。 //回代 y[0] = r[0] / d[0]。 for (i = 1。 i M。 i++) { y[i] = (r[i] a[i 1] * y[i 1]) / d[i]。 } x[M 1] = y[M 1]。 for (i = M 2。 i = 0。 i) { x[i] = y[i] u[i] * y[i + 1]。 } return x。}*/double* BVP_FD(double pfun(double), double qfun(double), double rfun(double), double t0, double alpha, double tn, double beta, int M){ int i = 0。 double h = (tn t0) / (M1)。 double t = t0。 int m = M 2。 double *a, *b, *c, *d, *u, *r, *y, *x。 a = (double *)malloc(m * sizeof(double))。 b = (double *)malloc(m * sizeof(double))。 c = (double *)malloc(m * sizeof(double))。 d = (double *)malloc(m * sizeof(double))。 u = (double *)malloc(m * sizeof(double))。 r = (double *)malloc(m * sizeof(double))。 y = (double *)malloc(m * sizeof(double))。 x = (double *)malloc(M * sizeof(double))。 t = t0 + h。 for (i = 0。 i m。 i++) { b[i] = 2 + h*h*qfun(t)。 r[i] = h*h*rfun(t)。 a[i] = h / 2 * pfun(t) 1。 c[i] = h / 2 * pfun(t) 1。 t = t + h。 } r[0] = r[0] ((h / 2. )*( pfun(t0+h)) 1.)*alpha。 r[m 1] = r[m 1] (h / 2 * pfun(tn h) 1)*beta。 d[0] = b[0]。 u[0] = c[0] / b[0]。 for (i = 1。 i m 1。 i++) { d[i] = b[i] a[i] * u[i 1]。 if (d[i] == 0) { return 0。 } u[i] = c[i] / d[i]。 } d[m 1] = b[m 1] a[m 1] * u[m 2]。 y[0] = r[0] / d[0]。 for (i = 1。 i m。 i++) { y[i] = (r[i] a[i] * y[i 1]) / d[i]。 } x[m] = y[m 1]。 for (i = m 2。 i = 0。i ) { x[i + 1] = y[i] u[i] * x[i + 2]。 } x[0] = alpha。 x[M 1] = beta。 return x。}double* BVP_Lshooting( double pfun(double), double qfun(double), double rfun(double), double t0,double a, double tn,double b, int M ){ return 0。}ifndef ODE_HHHHdefine ODE_HHHH// 2nd order RungeKutta method for solving initial value problemdouble* RungeKutta_Heum( double fun(double,double), double t0, double y0, double tn, int M)。// 4th order RungeKutta method for solving initial value problemdouble* RungeKutta4th( double fun(double,double), double t0, double y0, double tn, int M)。// ODE systems and high order ODE solverdouble** RKSystem4th( void fun(double,double[],double []), int N, double t0, double y0[],double tn, int M)。// ODE boundary value problem using shooting methoddouble** BVP_Shooting( void fun(double,double[],double[]), double t0,double a, double tn,double b, int M )。// BVP using finite difference methoddouble* BVP_FD( double pfun(double), double qfun(double), double rfun(double), double t0,double a, double tn,double b, int M )。endifVincent. Wang [24]201。quipe de Neutrons Dosim233。trie
環(huán)評(píng)公示相關(guān)推薦
鄂ICP備17016276號(hào)-1