c数值分析c17-6.pdfVIP

784 Chapter 17. Two Point Boundary Value Problems where φ(x) is chosen by us. Written in terms of the mesh variable q, this equation is dx ψ = (17.5.7) dq φ(x) Notice that φ(x) should be chosen to be positive definite, so that the density of mesh points is everywhere positive. Otherwise (17.5.7) can have a zero in its denominator. To use automated mesh spacing, you add the three ODEs (17.5.5) and (17.5.7) to your h r P C S t e t a e o a set of equations, i.e., to the array y[j][k]. Now x becomes a dependent variable! Q and ψ p d r p m : m y / a p / i r w b s i l also become new dependent variables. Normally, evaluating φ requires little extra work since w l s g e e i h p it will be composed from pieces of the g ’s that exist anyway. The automated procedure allows w f o t a . i n (