TAILIEUCHUNG - Integration of Ordinary Differential Equations part 3

dv=vector(1,nvar); for (i=1;i | 714 Chapter 16. Integration of Ordinary Differential Equations dv vector 1 nvar for i 1 i nvar i Load starting values. v i vstart i y i 1 v i xx 1 x1 x x1 h x2-x1 nstep for k 1 k nstep k Take nstep steps. derivs x v dv rk4 v dv nvar x h vout derivs if float x h x nrerror Step size too small in routine rkdumb x h xx k 1 x Store intermediate steps. for i 1 i nvar i v i vout i y i k 1 v i free_vector dv 1 nvar free_vector vout 1 nvar free_vector v 1 nvar CITED REFERENCES AND FURTHER READING Abramowitz M. and Stegun . 1964 Handbook of Mathematical Functions Applied Mathematics Series Volume 55 Washington National Bureau of Standards reprinted 1968 by Dover Publications New York . 1 Gear . 1971 Numerical Initial Value Problems in Ordinary Differential Equations Englewood Cliffs NJ Prentice-Hall Chapter 2. 2 Shampine . and Watts . 1977 in Mathematical Software III . Rice ed. New York Academic Press pp. 257-275 1979 Applied Mathematics and Computation vol. 5 pp. 93-121. 3 Rice . 1983 Numerical Methods Software andAnalysis New York McGraw-Hill . Stepsize Control forRunge-Kutta A good ODE integrator should exert some adaptive control over its own progress making frequent changes in its stepsize. Usually the purpose of this adaptive stepsize control is to achieve some predetermined accuracy in the solution with minimum computational effort. Many small steps should tiptoe through treacherous terrain while a few great strides should speed through smooth uninteresting countryside. The resulting gains in efficiency are not mere tens of percents or factors of two they can sometimes be factors of ten a hundred or more. Sometimes accuracy may be demanded not directly in the solution itself but in some related conserved quantity that can be monitored. Implementation of adaptive stepsize control requires that the stepping algorithm signal information about its performance most important an estimate of its truncation error. In this section we will

• Kỹ thuật lập trình
• tư liệu
• tài liệu lập trình
• thủ thuật lập trình
• ngôn ngữ C#
• ngôn ngữ lập trình
• thành phần cấu tạo vật liệu
• phân loại vật liệu compozit
• công nghệ chế tạo vật liệu compozit
• sản phẩm chế tạo từ vật liệu compozit
• phương hướng phát triển
• Tài liệu điện tử
• Vật liệu điện tử
• Linh kiện điện tử
• Vật liệu linh kiện
• Mạch điện tử
• Đề tài nghiên cứu khoa học
• Tổng hợp
• đặc trưng vật liệu nano
• Vật liệu nano cacbon
• Nano cacbon từ vỏ cua
• Xử lý môi trường
• Bài giảng Y học
• Alcaloid và Dược liệu chứa alcaloid
• Phân lập alc từ dược liệu
• Phân loại alcaloid
• Tính chất lý hóa
• Cấu trúc hóa học
• Tự động học
• Điều khiển học
• Điều khiển tự động
• Lý thuyết hệ thống
• Tài liệu tự động hóa