推导微分方程的解法和定理证明

中值定理证明解微分方程中值定理是微积分中的一个重要定理，它是解微分方程的基础。

本文将介绍中值定理的证明和如何利用它来解微分方程。

一、中值定理的证明中值定理也被称为罗尔定理或拉格朗日中值定理，它的表述如下：设函数 $f(x)$ 在闭区间 $[a,b]$ 上连续，在开区间 $(a,b)$ 内可导，则在 $(a,b)$ 内至少存在一点 $c$，使得 $f'(c)=\dfrac{f(b)-f(a)}{b-a}$。

证明如下：由于 $f(x)$ 在闭区间 $[a,b]$ 上连续，根据最值定理，$f(x)$ 在该区间内必有最大值 $M$ 和最小值 $m$，即 $m\leqf(x)\leq M$，且存在 $x_1,x_2\in[a,b]$，使得 $f(x_1)=M$，$f(x_2)=m$。

当 $f(x)$ 为常数函数时，结论显然成立。

当 $f(x)$ 不为常数函数时，存在 $x_0\in[a,b]$，使得 $f(x_0)\neq f(a)$，$f(x_0)\neq f(b)$。

不失一般性，假设 $f(x_0)>f(a)$。

若 $f(x_0)<f(b)$，则由连续性，$f(x)$ 在 $[a,b]$ 上取到$f(x_0)$ 的值，设为 $d$。

根据介值定理，存在 $[a,x_0]$ 和$[x_0,b]$ 上的某点 $c_1$ 和 $c_2$，使得 $f(c_1)=d$，$f(c_2)=d$。

由于 $f(x)$ 在 $(a,b)$ 内可导，根据导数的定义，有$$ f'(x_0)=\lim\limits_{h\to 0}\dfrac{f(x_0+h)-f(x_0)}{h} $$ 根据极限的性质，可以找到两个数列 $\{x_n\}$ 和 $\{y_n\}$，满足$$ \lim\limits_{n\to\infty}x_n=\lim\limits_{n\to\infty}y_n=x_ 0 $$$$ x_n\in(a,x_0),\ y_n\in(x_0,b) $$$$ f(x_n)<f(x_0),\ f(y_n)>f(x_0) $$于是有$$ f'(x_0)=\lim\limits_{n\to\infty}\dfrac{f(x_n)-f(x_0)}{x_n-x_0}=\lim\limits_{n\to\infty}\dfrac{f(y_n)-f(x_0)}{y_n-x_0} $$根据介值定理，存在 $\alpha\in[c_1,x_0]$ 和$\beta\in[x_0,c_2]$，使得 $f'(\alpha)=f'(\beta)=\dfrac{f(b)-f(a)}{b-a}$。

瑞利-里特尔定理：描述微分方程解的唯一性和稳定性瑞利-里特尔定理是微分方程理论中的一个重要定理，它描述了微分方程解的唯一性和稳定性。

瑞利-里特尔定理的数学表述非常精确，但对于非数学专业的人来说可能有些晦涩难懂。

本文将以通俗易懂的方式解释瑞利-里特尔定理，并给出一些实例来帮助读者更好地理解这个定理的意义和应用。

第一章：引言引言部分将介绍微分方程的基本概念和瑞利-里特尔定理的背景。

首先，我们将简要介绍微分方程的定义和分类，以及为什么微分方程在自然科学和工程学中如此重要。

然后，我们将提到瑞利-里特尔定理的历史背景和发展过程，以及该定理的重要性。

第二章：瑞利-里特尔定理的数学表述这一章节将详细阐述瑞利-里特尔定理的数学表述。

我们将引入微分方程的初值问题的概念，并定义初值问题的解的存在性、唯一性和稳定性。

然后，我们将给出瑞利-里特尔定理的数学表述，包括定理的前提条件和结论。

第三章：瑞利-里特尔定理的证明这一章节将介绍瑞利-里特尔定理的证明思路和主要步骤。

我们将从定理的前提条件出发，逐步推导出结论，并解释每一步的原理和推理过程。

在证明过程中，我们将用到一些基本的数学工具和技巧，如积分、微分等。

第四章：瑞利-里特尔定理的应用这一章节将介绍瑞利-里特尔定理在实际问题中的应用。

我们将给出一些具体的例子，如弹簧振子、电路中的振荡器等，来说明瑞利-里特尔定理对于解析解的存在性和稳定性的重要性。

我们还将讨论定理的局限性和适用范围，并提出一些改进和扩展的可能性。

第五章：总结与展望在这一章节中，我们将对本文进行总结，并展望瑞利-里特尔定理在未来的发展和应用。

我们将强调该定理在微分方程理论研究和实际问题求解中的重要性，并提出一些可能的研究方向和应用领域。

同时，我们也会回顾本文的主要内容和亮点，并给出一些对读者的启示和建议。

通过以上章节的划分和排版，本文将对瑞利-里特尔定理进行全面而深入的解释和阐述，使读者能够更好地理解该定理，并了解它在微分方程理论和实际问题中的重要性和应用。

求数学公式的11种推导方法在数学中，推导公式是一种常见的方法，它可以帮助我们理解数学原理和解决问题。

本文将介绍11种常用的数学公式推导方法。

1. 直接证明法直接证明法是最常见的推导方法之一。

它通过从已知的前提出发，逐步推导出所要证明的结论。

这种方法通常是通过逻辑推理和数学运算来完成的。

2. 反证法反证法是一种通过假设某个结论为假，然后导出逻辑矛盾的方法来推导公式。

如果我们能够证明该假设是错误的，那么所要证明的结论就是对的。

3. 数学归纳法数学归纳法是一种证明递归定义上成立的方法。

它通常分为两个步骤：基础情况的证明和归纳步骤的证明。

4. 同余模运算同余模运算是一种推导数学公式的方法，它基于模运算的性质进行推导。

这种方法通常用于证明数论中的一些定理和公式。

5. 极限和极限运算极限和极限运算是一种通常用于推导数学公式的方法。

通过计算函数的极限，我们可以推导出一些公式，例如泰勒展开式和级数求和公式。

6. 向量分析向量分析是一种用于推导数学公式的方法，它基于向量运算和坐标系的概念。

通过对向量进行运算和变换，我们可以推导出许多与几何和物理相关的公式。

7. 矩阵运算矩阵运算是一种用于推导数学公式的方法，它基于矩阵的性质和运算规则。

通过对矩阵进行运算和变换，我们可以推导出许多与线性代数和线性方程组相关的公式。

8. 微积分微积分是一种用于推导数学公式的方法，它基于导数和积分的概念。

通过对函数进行微分和积分，我们可以推导出许多与曲线，曲面和体积相关的公式。

9. 概率论和统计学推导概率论和统计学是一种用于推导数学公式的方法，它基于概率和统计的概念。

通过对随机变量和概率分布进行分析，我们可以推导出许多与概率和随机过程相关的公式。

10. 微分方程推导微分方程是一种用于推导数学公式的方法，它基于微分方程的性质和解法。

通过对微分方程进行求解和变换，我们可以推导出许多与动力学和振动系统相关的公式。

11. 几何推导几何推导是一种用于推导数学公式的方法，它基于几何的性质和定理。

2.7．微分方程初步2.7.1 概说涉及到量的变化率满足的制约关系，通常是含有导数的方程——微分方程。

涉及到量的变化率满足的制约关系，通常是含有导数的方程——微分方程。

简单例子：（1）放射性物质，在每一时刻t ，衰变的速率d m d t-（由于是减少，因此0d m d t<，速率为标量，是正值）正比于该放射性物质尚存的质量，因此质量应满足一下微分方程。

标量，是正值）正比于该放射性物质尚存的质量，因此质量应满足一下微分方程。

d m km d t-=（2）质量为m 的物体自由落体，取坐标轴沿竖直方向指向地心，下落距离()y y t =应该满足牛顿第二定律F m a =，即，即22d y m g mdt=（3）质量为m 的跳伞员下落，所受空气阻力正比下降的速度，取坐标轴沿竖直方向指向地心，则t 时刻下降距离()y y t =满足满足22dyd y m g kmdtdt-=（1）如下图所示，钢球在以水平光滑杆上，受到弹力而来回整栋，原点位置为O，钢球在t 时刻的坐标()x x t =满足微分方程满足微分方程()22d x kx mdt-=如果钢球还受到一个与速度成正比，方向与速度相反的阻尼力的作用，那么它所满足的微分方程是方程是22dxd x kx h m dtdt--=总结：最简单的一阶微分方程是最简单的一阶微分方程是()d x f t d t=其中t 是自变量，上述方程的一般解应该是是自变量，上述方程的一般解应该是()x f t dt C =+ò最简单的n 阶方程阶方程()nnd xf t dt=它等价于说11n nn d x dt--是()f t 的原函数，即的原函数，即11()n n ndxf t dt C dt --=+ò则再次积分，一直积分下去得到则再次积分，一直积分下去得到111()(1)!n nn n t x f t dt dt C C t C n --=++++-òò2.7.2 一阶线性微分方程考察下面的方程考察下面的方程()()d x a t x b t d t+=方程中有未知函数的一阶导数，且其一阶导数的系数为常数，其余部分未知函数最高层次数为一次，称为线性，上述方程为一阶线性微分方程。

第九章微分方程第一节基本概念一.解释下列名词术语1．微分方程：含有未知函数的导数（或微分）的方程.注意：（1）微分方程的一般形式：，在这个方程中是自变量，是的未知函数，是对的一阶、二阶、n阶导数；（2）方程中未知函数及自变量的记号可以不出现，如：；但未知函数的导数则必须出现.2.微分方程的阶：微分方程中所含的未知函数的最高阶导数的阶数.如：是一阶是二阶是n阶3.微分方程的解：代入微分方程能使方程成为恒等式的函数.例如：是的解.4.微分方程的通解：n阶微分方程的含有n个独立的任意常数的解.例如：是的通解；但是的解，而非通解.注意：这里要说明一下“两个常数独立”的含义----即对于任意给定的不同的的取值，则应得到不同的解，则称两个常数是互相独立的.之所以不是的通解，就是因为不是互相独立的.比如：取或者都可得到解.5.微分方程的初始条件：用来确定通解中的任意常数的一种定解的条件.一阶微分方程的初始条件通常为二阶微分方程的初始条件通常为例如：已知是的通解，可由初始条件通常为。

初始条件的个数与微分方程的阶数相同。

6.微分方程的特解：通解中所含的所有任意常数都确定后的解。

比如：是的满足初始条件的特解。

7.积分曲线：微分方程的解的图形（特解是一条积分曲线；通解是一组积分曲线）二。

用微分方程求解实际问题中的未知函数的步骤：1.建立微分方程和初始条件（难点）；------这通常使一部分同学感到为难，因为它除了需要数学知识之外，还往往要用到力学、物理学、化学、电学、工程技术等方面的知识，甚至还要用到语文的知识。

2.求通解；3.求特解。

我们这一章的重点是：给定一个微分方程，如何求其通解或特解.第二节一阶微分方程一.可分离变量的微分方程求解微分方程有一个特点：就是“对号入座”,什么样的微分方程，就用什么方法去解决,这几乎成了一个固定的格式.因此，判定所给的方程是什么类型就是首要问题。

这是本章的特点.今天，就给大家介绍一种最简单的一阶微分方程：可分离变量的微分方程.1.引例求解解：因为，所以是是的一个原函数。

第二章 基本定理我们在第一章主要学习了初等积分法，掌握了几类常微分方程的解法.但是这些解法只适用于某些特殊的类型，很多其它的常微分方程不能用初等解法进行求解.1841年，法国数学家刘维尔（Liouville ）证明了里卡蒂（Riccati ）方程)0)(()()()(2≠++=x p x r y x q y x p dydx 除了某些特殊的类型外，一般不能用初等积分法求解.例如，很简单的里卡蒂方程22y x dxdy +=就不能用初等积分法求解.自然地，如果一个常微分方程不能用初等积分法求解，那么应该如何处理呢？是否存在解呢？如果存在解，它的解是否唯一呢？解的存在区间是什么呢？初值的微小误差对解有什么影响呢？这些问题在理论的研究和实际应用中，都有着重要的意义.本章将解决这些基本问题. 本章主要介绍解的存在唯一性定理、解的延展定理与比较定理、解对初值的连续依赖性定理以及解对初值的可微性定理，这些定理就回答了我们刚才的疑问，有效的处理解的存在性、唯一性、存在区间、初值对解的影响等问题，为我们使近似解法奠定理论基础，同时这些定理也是常微分方程理论的基础内容，对进一步的学习奠定基础.解的存在唯一性定理对于一般的常微分方程),(y x f dxdy = （） 如果给出了初始条件00)(y x y =，我们就得到了柯西初值问题⎪⎩⎪⎨⎧==00)(),(y x y y x f dx dy （） 这时，在什么样的条件下，柯西初值问题的解存在且唯一呢？解的存在区间是什么呢？我们有如下的解的存在唯一性定理.2.1.1 存在唯一性定理的叙述定理（存在唯一性定理）如果方程（）的右端函数),(y x f 在闭矩形区域b y y b y a x x a x R +≤≤-+≤≤-00002,:上满足如下条件：（1）在2R 上连续；（2）在2R 上关于变量y 满足李普希兹（Lipschitz ）条件，即存在常数N ，使对于2R 上的任何一对点),(y x 和),(y x 有不等式：y y N y x f y x f -≤-),(),(则初值问题（）在区间],[0000h x h x +-上存在唯一解00)(),(y x x y ==ϕϕ 其中),(max ),,min(),(0y x f M M b a h R y x ∈==.在给出定理的证明之前，我们先对定理的条件和结论做些说明：1、在两个条件中，条件（2）,即李普希兹条件比较难于验证，因为李普希兹常数N 难以确定.但是，我们可以将该条件加强，替换为：如果函数),(y x f 在闭矩形区域2R 关于y 的偏导数),(y x f y '存在且有界.这样，可以推出李普希兹条件成立.事实上，因为),(y x f y '有界，故设N y x f y ≤'),(，对2),(),,(R y x y x ∈∀，由拉格朗日中值定理得：y y N y y x f y x f y x f y -≤-'=-),(),(),(ξ我们验证),(y x f y '在闭矩形区域2R 上有界也不容易，可以进一步将条件加强为：),(y x f y '在闭矩形区域2R 上连续.由闭区域上连续函数的性质知：),(y x f y '在闭矩形区域2R 上有界，所以李普希兹条件成立.因此，有如下的关系式：),(y x f y '在2R 上连续⇒),(y x f y '在2R 上存在且有界⇒李普希兹条件2、在定理的结论中，解)(x y ϕ=的存在区间为],[0000h x h x +-，其中),(max ),,min(),(0y x f M M b a h R y x ∈==.为什么解的存在区间不是],[00a x a x +-呢？这是因为我们研究问题的范围为闭矩形区域2R ，方程的解)(x y ϕ=不能超出2R 的范围，又因为),(max ),(y x f M Ry x ∈=，所以M y x f M ≤≤-),( 即 M dxdy M ≤≤- 由⎪⎩⎪⎨⎧=-=00)(y x y M dx dy 和⎪⎩⎪⎨⎧==00)(y x y M dx dy 得：001)()(y x x M x y +--=，002)()(y x x M x y +-= 因此)()()(21x y x y x y ≤=≤ϕ，即)(x y ϕ=夹在)(1x y 与)(2x y 之间.又，)(1x y 与)(2x y 在2R 上的存在区间为],[0000h x h x +-， 故)(x y ϕ=的存在区间也是],[0000h x h x +-.2.1.2 存在性的证明首先，我们给出柯西初值问题（）的等价转化，即求（）的解)(x y ϕ=，等价于求解积分方程 ⎰+=xx d y f y y 0))(,(0ξξξ （） 事实上，如果)(x y ϕ=是初值问题（）的解，即有))(,()(x x f x ϕϕ='且00)(y x =ϕ从0x 到x 积分得：⎰+=xx d f y x 0))(,()(0ξξϕξϕ 即)(x y ϕ=是积分问题（）的解.反过来，如果)(x y ϕ=是积分问题（）的解，即有⎰+=xx d f y x 0))(,()(0ξξϕξϕ 则00)(y x =ϕ且))(,()(x x f x ϕϕ='即)(x y ϕ=是初值问题（）的解.经过等价转化，我们将初值问题（）的求解，转化为积分问题（）的求解. 下面用皮卡（Picard ）逐次逼近来证明积分问题（）的解的存在性，分为三个步骤：1、构造近似函数列{})(x n ϕ任取一个满足初值条件00)(y x y =的函数)(0x y ϕ=作为首项（初始项），并要求在2R 上的存在区间为：],[0000h x h x +-，简单起见，取00)(y x =ϕ，将它代入方程（）的右端，所得到的函数用)(1x ϕ表示，并称为一次近似，即⎰+=xx d f y x 0))(,()(001ξξϕξϕ 再将)(1x ϕ代入方程（）的右端就得到二次近似⎰+=xx d f y x 0))(,()(102ξξϕξϕ 序行此法，可以得到n 次近似⎰-+=xx n n d f y x 0))(,()(10ξξϕξϕ 为了保证上述的逐次逼近过程可以一直进行下去，必须有2))(,(R x x n ∈ϕ，即当],[0000h x h x x +-∈时，有,2,1)(0=≤-n b y x n ϕ 下面用数学归纳法证明b y x n ≤-0)(ϕ.显然，当],[0000h x h x x +-∈时，有b y y y x ≤=-=-0)(0000ϕ假设，当],[0000h x h x x +-∈时，有b y x n ≤--01)(ϕ，那么，对于)(x n ϕ有⎰-=-xx n n d f y x 0))(,()(10ξξϕξϕ从而有b M b M Mh x x M d f y x x x n n =≤≤-≤≤-⎰-00100))(,()(ξξϕξϕ 由数学归纳法知，当],[0000h x h x x +-∈时，有 ,2,1)(0=≤-n b y x n ϕ这样，我们就可以得到一个近似函数列{})(x n ϕ.2、证明近似函数列{})(x n ϕ在区间],[0000h x h x +-上一致收敛.由于无法得到{})(x n ϕ的通项公式，只知道首项和递推关系式，直接证明函数列{})(x n ϕ的收敛性比较困难，为此我们构造函数项级数+-++-+-)]()([)]()([)(1010x x x x x n n ϕϕϕϕϕ （） 它的部分和是)()]()([)]()([)()(10101x x x x x x x S n n n n ϕϕϕϕϕϕ=-++-+=-+因此，证明{})(x n ϕ的收敛性转化为证明级数（）的收敛性，下面我们证明级数（）在区间],[0000h x h x +-上一致收敛.首先研究级数（）的通项)(x n μ⎰=-xx d f x x 0))(,()()(001ξξϕξϕϕ 即⎰=-xx d y f y x 0),()(001ξξϕ 所以00010),()(x x M d y f y x x x -≤≤-⎰ξξϕ 因为⎰+=x x d f y x 0))(,()(001ξξϕξϕ，⎰+=x x d f y x 0))(,()(102ξξϕξϕ，所以 ⎰-≤-x x d f f x x 0))(,())(,()()(0112ξξϕξξϕξϕϕ由李普希兹条件，得 !2)()()()(200011200x x MN d x MN d N x x x x x x -=-≤-≤-⎰⎰ξξξξϕξϕϕϕ 下面用数学归纳法证明!)()(011n x x MN x x nn n n -≤---ϕϕ 显然，2,1=n 的时候，不等式成立（上面已经给出）， 假设!)()(011n x x MN x x n n n n -≤---ϕϕ成立，那么对于1+n 的情形有 )!1(!)()())(,())(,()()(100111000+-=-≤-≤-≤-+--+⎰⎰⎰n x x MN d n x MN d N d f f x x n n x x n n xx n n x x n n n n ξξξξϕξϕξξϕξξϕξϕϕ由数学归纳法知，对一切自然数n ，均有!)()(011n x x MN x x nn n n -≤---ϕϕ又00h x x ≤-，所以级数（）的通项满足：!)(011n h MN v x n n n n -+=≤μ （ ,2,1=n ） 利用比式判别法，可知以n v 为通项的级数收敛，从而以)(x n μ为通项的级数绝对收敛且一致收敛.又，每一个)(x n μ是连续的，所以级数（）的和函数也是连续的，记为)(x ϕ，其存在区间也是],[0000h x h x +-.因此函数列{})(x n ϕ就收敛于)(x ϕ.3、证明)(lim )(x x n n ϕϕ∞→=是积分问题（）的解，从而也是初值问题（）的解. 在⎰-+=x x n n d f y x 0))(,()(10ξξϕξϕ两端取极限，得到 ⎰-∞→∞→+=xx n n n n d f y x 0))(,(lim )(lim 10ξξϕξϕ 即⎰+=xx d f y x 0))(,()(0ξξϕξϕ 所以)(x ϕ是积分问题（）的解，从而也是初值问题（）的解.2.1.3 唯一性的证明下面我们证明解的唯一性.在证明唯一性之前，先介绍一个重要的不等式，即贝尔曼（Bellman ）不等式.贝尔曼引理 设)(x y 为区间],[b a 上的非负连续函数，b x a ≤≤0.若存在,0≥δ 0≥k ，使得)(x y 满足不等式],[,)()(0b a x d y k x y xx ∈+≤⎰ττδ （） 则有],[,)(0b a x e x y x x k ∈≤-δ证明 仅证明0x x ≥的情形，0x x ≤的情形类似.令)(x y 的原函数为⎰=xx d y x R 0)()(ττ，代入（）得 δ≤-')()(x kR x R两边同时乘以积分因子)(0x x k e --，得)()(00)]()([x x k x x k e x kR x R e ----≤-'δ从0x 到x 积分得)()(00)(x x k x x k e e x kR -----≤δδ即)(0)(x x k e x kR -≤+δδ 由（）知，)()(x kR x y +≤δ，所以],[,)(0b a x e x y x x k ∈≤-δ下面证明积分问题（）的解的唯一性.假设积分问题有两个解)(1x y 和)(2x y ，我们只需要证明：)(1x y )(2x y ≡，],[0000h x h x x +-∈事实上，因为 ⎰+=x x d y f y x y 0))(,()(101ξξξ，⎰+=xx d y f y x y 0))(,()(202ξξξ 所以有⎰-≤-x x d y f y f x y x y 0))(,())(,()()(2121ξξξξξ由李普希兹条件知 ⎰-≤-xx d y y N x y x y 0)()()()(2121ξξξ 令N k x y x y x y ==-=,0,)()()(21δ，由贝尔曼引理可知，0)(=x y ，即)(1x y )(2x y ≡.这样，我们就完成了解的存在性与唯一性的证明.2.1.4 三点说明为了更好的理解和掌握解的存在唯一性定理，我们对该定理再做三点说明.1、在存在性的证明过程中，我们利用逐次逼近法构造了近似函数列{})(x n ϕ，其中首项为：00)(y x =ϕ，递推关系式为：⎰-+=xx n n d f y x 0))(,()(10ξξϕξϕ.该方法实际上给出了我们一种求初值问题（）的近似解的方法，当用n 次近似解逼近精确解时，需要给出它的误差估计.事实上，有∑∑∞+=∞=+-≤-≤-101!)()()()(n k kk n k k k n k x x N N M x x x x ϕϕϕϕ 0)!1()(!)!1()(!10001010Nh n k kk n n k k k e n Nh N M k h N n Nh N M k h N NM +=+<≤+∞=+∞+=∑∑2、如果方程（）是线性方程，即)()(x q y x p dxdy +-= 其中)(x p 和)(x q 在区间],[b a 上连续，这时，初值问题（）在带型区域+∞<<-∞≤≤y b x a R ,:2满足定理的条件.事实上，)()(),(x q y x p y x f +-=在2R 上连续，而且)(),(x p y x f y -='在2R 上也连续，所以),(y x f 关于变量y 满足李普希兹条件.这时，初值问题（）的解存在且唯一，存在区间为],[b a .3、定理中的李普希兹条件是保证解唯一的充分条件，那么这个条件是不是必要条件呢？回答是否定的，即李普希兹条件是解唯一的充分非必要条件.下面我们给出一个例子来说明李普希兹条件是解唯一的非必要条件，也就是说，即使李普希兹条件不成立，初值问题（）的解也可能是唯一的.例1 试证方程0,ln ,0≠=⎩⎨⎧=y y y y dx dy 经过xOy 平面上任一点的解都是唯一的.证明 由00,ln ,0≠=⎩⎨⎧=y y y y dx dy 可得：0=y 或x Ce e y ±=. 任给xOy 平面上的一个点),(00y x ，只会对应0=y 或xCe e y ±=中的一个解，也就是说，过xOy 平面上任一点的解都是唯一的.但是，我们有 0ln ln )0,(),(-==-y y y y x f y x f 因为+∞=→y y ln lim 0，所以找不到0>N ，使得 0)0,(),(-≤-y N x f y x f从而方程右端函数在0=y 的任何邻域上不满足李普希兹条件，但是初值问题（）的解却是唯一的，这说明李普希兹条件是非必要条件.习 题1.试判断方程y x dxdy tan =在区域 （1）π≤≤≤≤-y x R 0,11:1；（2）44,11:2ππ≤≤-≤≤-y x R上是否满足定理的条件？2.讨论方程3123y dx dy =在怎样的区域中满足定理的条件.并求通过)0,0(的一切解.3.试用逐次逼近法求方程2y x dxdy -=满足初值条件0)0(=y 的近似解： )(),(),(),(3210x x x x ϕϕϕϕ并在闭矩形区域11,11:2≤≤-≤≤-y x R 给出三次近似的误差估计.4.利用逐次逼近法求方程22x y dxdy -=适合初值条件1)0(=y 的近似解： )(),(),(210x x x ϕϕϕ并在闭矩形区域111,11:2≤-≤-≤≤-y x R 给出二次近似的误差估计.5.试证明定理中的n 次近似解)(x n ϕ与精确解)(x ϕ有如下的误差估计式：10)!1()()(+-+≤-n n n x x n MN x x ϕϕ 6.在条形区域+∞<≤≤y b x a ,内，假设方程（）的所有解都唯一，对其中任意两个解)(),(21x y x y ，如果有)()(0201x y x y <，则必有b x x x y x y ≤≤<021),()(.7.讨论方程323y dx dy = 解的唯一性.延展定理和比较定理由解的存在唯一性定理，我们知道，初值问题（）的解在满足一定条件的情况下存在且唯一，但是解的存在区间不是],[00a x a x +-，而是],[0000h x h x +- 其中),(max ),,min(),(0y x f M M b a h R y x ∈==.如果M 比较大的话，则解的存在区间就非常小，这对我们研究解的性质产生了很大的局限性，只能在很小的范围内有解，当x 超出这个范围时，解的情况就不清楚了.为了解决这个问题，我们有下面的延展定理.2.2.1 延展定理定理（延展定理）如果方程（）的右端函数在区域R R D ⨯⊂上连续，且关于变量y 满足局部的李普希兹条件，即对于D 内的任一闭矩形区域都满足李普希兹条件，则对任何一点D y x ∈),(00，初值问题（）的解)(x y ϕ=可以向左右无限延展，直到))(,(x x ϕ任意接近区域D 的边界.在给出定理的证明之前，先对“))(,(x x ϕ任意接近区域D 的边界”进行说明.当区域D 有界时，积分曲线向左右延展可以任意接近；当区域D 无界时，积分曲线向左、右延展，或者任意接近区域D 的边界（边界存在的话），或者无限远离坐标原点.证明 首先证明区域D 有界的情形.设区域D 的边界为D D L -=（D 为D 的闭包）.对于任意给定的正数ε，记L 的ε邻域为εU ，记L 的2ε邻域为2εU ，记L 的4ε邻域为4εU .则集合22εεU D D -=为闭集，且D D ⊂2ε，所以2εD 有界. 只要证明积分曲线可以到达2εD 的边界2εL ，由ε的任意性知，积分曲线就可以任意接近区域D 的边界L .事实上，以2εD 中的任意一点为中心，以4ε为半径的闭圆区域均包含在区域D 的内部.且在闭区域44εεU D D -=之内.从而，以2εD 中的任意一点为中心，以4221ε=a 为边长的正方形也在闭区域4εD 之内.记 ),(max 4),(1y x f M D y x ε∈=则过2εD 的任意一点),(**y x 的积分曲线，必至少可在区间],[**h x h x +-上存在，其中)82,82min(),min(1111M M a a h εε==. 于是，过点),(00y x 的积分曲线)(x y ϕ=每向左或向右延展一次，其存在区间就伸长一个确定的正数h ，由于2εD 有界，)(x y ϕ=经过有限次延展后一定可以达到2εD的边界2εL .于是也就可以任意接近区域D 的边界L .其次考虑区域D 为无界的情形.这时，我们可以用闭圆区域,2,1},),{(222=≤+=n n y x y x S n与区域D 取交集，令n n S D D =，则 ∞==1n n D D .由于n D 为有界的区域，根据前面的证明，我们可知，过n D 内任一点的积分曲线能够任意接近n D 的边界.因此，过点),(00y x 的积分曲线)(x y ϕ=可以无限接近区域D 的边界.延展定理的证明，关键是第一步证明，也就是区域D 有界的时候，过点),(00y x 的积分曲线)(x y ϕ=向左向右延展的时候，一定要做等速延展，即延展步幅h 是不变的. 例1 试讨论方程2y dxdy=通过点)1,1(的解和通过点)1,3(-的解的存在区间. 解 该题目中研究问题的区域D 为整个坐标平面xOy .方程右端函数满足延展定理的条件.由2y dxdy=可以解得方程的通解为 xC y -=1代入1)1(=y 得：2=C .故通过点)1,1(的解为xy -=21 它可以向左无限延展，而当-→2x 时，+∞→y ，所以通过点)1,1(的解xy -=21的存在区间为)2,(-∞.代入1)3(-=y 得：2=C .故通过点)1,3(-的解为xy -=21 它可以向右无限延展，而当+→2x 时，-∞→y ，所以通过点)1,3(-的解xy -=21的存在区间为),2(+∞.这个例子说明，尽管),(y x f 在整个坐标平面上满足延展定理的条件，解上的点))(,(x x ϕ也能无限接近区域D 的边界，但是延展的方向却不一定是无限向右和向左，可能是向上或向下，从而导致解的存在区间不是),(+∞-∞. 例2 试证明：对任意的0x 及满足条件100<<y 的0y ，方程221)1(y x y y dx dy ++-=的满足条件00)(y x y =的解)(x y y =在),(+∞-∞上存在.证明：令221)1(),(yx y y y x f ++-=，则 222222)1(122),(y x x y y x y y x f y ++--++=' 显然),(),,(y x f y x f y '在xOy 平面上连续，满足解的存在唯一性条件及延展定理的条件，而1,0==y y 是),(y x f dxdy=的解， 因此，满足00)(y x y =，100<<y 的解存在，而且可以无限延展到xOy 平面的边界，且不能穿过1,0==y y ，故只能向左右无限延展，所以，)(x y y =在),(+∞-∞上存在.该例题说明，),(y x f 在整个坐标平面上满足延展定理的条件，当方程的解不能穿过1,0==y y 时，它就不能向上向下无限延展了，只能向左、向右延展，所以解的存在区间就是),(+∞-∞.在这里，1,0==y y 控制了解的延展方向，使它按照我们的要求进行延展，因此就有了下面的比较定理. 2.2.2 比较定理我们在使用延展定理的时候，通常会和比较定理配合使用，从而起到控制延展方向的作用.下面介绍一下比较定理.我们在考察方程（）),(y x f dxdy= 时，通常将右端函数),(y x f 进行放缩的处理，比如),(),(),(21y x F y x f y x F <<这时，我们可以同时考察),(1y x F dx dy =和),(2y x F dxdy = 我们有如下的比较定理：定理 （第一比较定理）设定义在某个区域D 上的函数),(y x f ，),(1y x F 和),(2y x F 满足条件：（1）在D 满足解的存在唯一性定理及延展定理的条件，即在D 上连续，在D 上关于变量y 满足李普希兹条件； （2）在D 上有不等式),(),(),(21y x F y x f y x F <<设初值问题⎪⎩⎪⎨⎧==00)(),(y x y y x f dxdy ，⎪⎩⎪⎨⎧==001)(),(y x y y x F dx dy 和⎪⎩⎪⎨⎧==002)(),(y x y y x F dx dy的解分别为)(x y ϕ=，)(1x y Φ=和)(2x y Φ=，则在它们的共同存在区间上有下列不等式：021),()()(x x x x x >Φ<<Φϕ 021),()()(x x x x x <Φ>>Φϕ证明 仅证当0x x >时，)()(2x x Φ<ϕ，其它的情形相类似. 由比较定理的条件（1），初值问题⎪⎩⎪⎨⎧==00)(),(y x y y x f dxdy 和⎪⎩⎪⎨⎧==002)(),(y x y y x F dx dy的解在0x 的某一邻域内存在且唯一，分别记为)(x y ϕ=和)(2x y Φ=，它们满足0020)()(y x x =Φ=ϕ令)()()(2x x x h ϕ-Φ=，则0)()()(0020=-Φ=x x x h ϕ且0))(,())(,()()()(0002020020>-Φ='-Φ'='x x f x x F x x x h ϕϕ所以函数)(x h 在0x 的某一右邻域内是严格单调增加的.如果在0x x >时，0)(>x h 不是总成立，则至少存在一点01x x >，使得0)(1=x h ，且当10x x x <<时，0)(>x h ，因此在点1x 的左导数0)0(1≤-'x h ，这与0))(,())(,()()()(1112121121>-Φ='-Φ'='x x f x x F x x x h ϕϕ矛盾.因此当0x x >时，0)(>x h 总成立，即)()(2x x Φ<ϕ.比较定理的应用，关键是),(1y x F 和),(2y x F 的选取，因为初值问题⎪⎩⎪⎨⎧==00)(),(y x y y x f dxdy的解)(x y ϕ=的存在区间的延展，受到)(1x y Φ=和)(2x y Φ=的控制，即)(x y ϕ=夹在)(1x y Φ=和)(2x y Φ=之间.因此，我们必须能确定出)(1x y Φ=和)(2x y Φ=的存在区间，这就是我们选取),(1y x F 和),(2y x F 的标准，即⎪⎩⎪⎨⎧==001)(),(y x y y x F dxdy 和⎪⎩⎪⎨⎧==002)(),(y x y y x F dx dy的解)(1x y Φ=和)(2x y Φ=必须能够求得. 下面我们给出第二比较定理.定理 （第二比较定理）设定义在某个区域D 上的函数),(y x f ，),(1y x F 和),(2y x F 满足条件：（1）在D 满足解的存在唯一性定理及延展定理的条件，即在D 上连续，在D 上关于变量y 满足李普希兹条件； （2）在D 上有不等式),(),(),(21y x F y x f y x F ≤≤设初值问题⎪⎩⎪⎨⎧==00)(),(y x y y x f dxdy ，⎪⎩⎪⎨⎧==001)(),(y x y y x F dx dy 和⎪⎩⎪⎨⎧==002)(),(y x y y x F dx dy的解分别为)(x y ϕ=，)(1x y Φ=和)(2x y Φ=，则在它们的共同存在区间上有下列不等式：021),()()(x x x x x >Φ≤≤Φϕ 021),()()(x x x x x <Φ≥≥Φϕ习 题1．设方程为),()(22y x f a y dxdy-= 假设),(y x f 及),(y x f y '在xOy 平面上连续，试证明：对于任意的0x 及a y <0，方程满足00)(y x y =的解都在),(+∞-∞上存在.2.指出方程2)1(2xy e y dxdy-=的每一个解的最大存在区间，以及当x 趋于这个区间的右端点时解的极限.3.讨论方程xx dx dy 1cos 12-= 解的存在区间.4.设),(y x f 在整个平面上连续有界，对y 有连续偏导数，试证明方程),(y x f dxdy=的任一解)(x y ϕ=在区间+∞<<∞-x 上有定义. 5.讨论方程212-=y dx dy 的通过点)0,0(的解，以及通过点)3,2(ln -的解的存在区间.6.在方程)(y f dxdy=中，如果)(y f 在),(+∞-∞上连续可微，且 )0(0)(≠<y y yf ，求证方程满足00)(y x y =的解)(x y 在区间),[0+∞x 上存在，且有0)(lim =+∞→x y x .解对初值的连续依赖性定理和解对初值的可微性定理通过前两节的存在唯一性定理和延展定理，加上比较定理，我们知道了初值问题（）在什么样的条件下，解是存在的，是唯一的，而且存在区间比较小的时候，通过延展定理和比较定理可以将解的存在区间变大，从而在实际问题中可以达到我们的要求.但是，在实际问题中，还有一个问题需要解决，那就是误差问题.我们的初始条件0)(y x y =如果产生了微小的偏差，这个偏差对我们的初值问题（）的解)(x y ϕ=会有什么影响呢？下面我们来解决这个问题. 我们在研究初值问题（）的时候，习惯上把x 和y 当作常数来看待，这样初值问题（）的解)(x y ϕ=被看作x 的函数.实际上，如果0x ，y 变化，初值问题（）的解)(x y ϕ=也会发生变化.例如方程x y dx dy =经过点),(00y x 的解为x x y y 0=，可以看作0,,y x x 的函数.对于一般的情形，初值问题（）的解也可以看作00,,y x x 的函数，记为),,(00y x x y ϕ=，代入0)(y x y =得：0000),,(y y x x =ϕ.如果我们的初始条件00)(y x y =发生了微小的误差，变为了*0*0)(y x y =，初值问题（）的解也变化不大的话，称解连续依赖于初值.下面我们给出连续依赖性的严格定义.定义 设初值问题⎪⎩⎪⎨⎧==*0*0)(),(y x y y x f dxdy的解),,(*0*0y x x y ϕ=在区间],[b a 上存在，如果对于任意给定的正数ε，存在正数δ(δ的选取与,ε*0*0,y x 有关)，使得对于满足δδ<-<-*00*00,y y x x （）的解),,(00y x x y ϕ=都在],[b a 上存在，且有],,[,),,(),,(*0*000b a x y x x y x x ∈<-εϕϕ则称初值问题（）的解),,(00y x x y ϕ=在点),(*0*0y x 连续依赖于初值,0x 0y .定理 （解对初值的连续依赖性定理）设),(y x f 在区域D 内连续，且关于变量y 满足李普希兹条件.如果D y x ∈),(*0*0，初值问题（）有解),,(*0*0y x x y ϕ=,且当b x a ≤≤时，D y x x x ∈)),,(,(*0*0ϕ，则对任意的正数ε，存在0>δ，使对于满足δδ<-<-*00*00,y y x x的任意),(00y x ，初值问题⎪⎩⎪⎨⎧==00)(),(y x y y x f dx dy的解),,(00y x x y ϕ=也在区间],[b a 上存在，且有εϕϕ<-),,(),,(*0*000y x x y x x证明 对于任意给定的正数ε，取εδ<<10，使得闭区域}),,(,),{(1*0*0δϕ≤-≤≤=y x x y b x a y x U整个含在区域D 内，这是可以做到的，因为区域D 是开区域，且当b x a ≤≤时，Dy x x x ∈)),,(,(*0*0ϕ，所以，只要1δ的选取足够小，以曲线),,(*0*0y x x y ϕ=为中线，宽度为12δ的带形开区域U 就整个包含在区域D 内， 选取δ满足)(110a b N e M--+<<δδ其中N 为李普希兹常数，),(max ),(y x f M Uy x ∈=，同时还要求δ的选取，必须保证闭正方形δδ≤-≤-*0*02,:y y x x R含于带形开区域U 内.由存在唯一性定理知，对于任一200),(R y x ∈，初值问题（）在0x 的某邻域上存在唯一解),,(00y x x y ϕ=，而且),,(00y x x y ϕ=在x 的该邻域上可以表示为ττϕτϕd y x f y y x x xx )),,(,(),,(000000⎰+=而),,(*0*0y x x y ϕ=可以表示为ττϕτϕd y x f y y x x xx )),,(,(),,(*0*0*0*0*0*⎰+=对上述两式做差得：ττϕτττϕτϕϕd y x f d y x f y y y x x y x x xx xx )),,(,()),,(,(),,(),,(*0*000*00*0*000*⎰⎰-+-=-所以ττϕτττϕτϕϕd y x f d y x f y y y x x y x x xx x x )),,(,()),,(,(),,(),,(*0*000*00*0*000*⎰⎰-+-≤-ττϕτττϕττϕτd y x f d y x f y x f y y x x xx |)),,(,(||)),,(,()),,(,(|0000*0*0*00*0*0⎰⎰+-+-≤δττϕττϕτδM d y x f y x f xx +-+≤⎰|)),,(,()),,(,(|00*0*0*0ττϕτϕδd y x y x N M xx |),,(),,(|)1(00*0*0*0-++≤⎰由贝尔曼引理，得εδδδϕϕ<<+≤+≤---1)(*0*000)1()1(),,(),,(*a b N x x N e M eM y x x y x x因此，只要在),,(00y x x y ϕ=有定义的区间上，就有εϕϕ<-),,(),,(*0*000y x x y x x .下面我们证明：),,(00y x x y ϕ=在区间],[b a 上有定义.事实上，因为εϕϕ<-),,(),,(*0*000y x x y x x即解),,(00y x x y ϕ=夹在εϕ+=),,(*0*0y x x y 和εϕ-=),,(*0*0y x x y 之间，而且，初值问题（）满足延展定理的条件，所以，解),,(00y x x y ϕ=可以向左向右无限延展，直到无限接近区域D 的边界，于是，它在延展的时候，必须由直线a x =和直线b x =穿出区域U ，从而),,(00y x x y ϕ=在区间],[b a 上有定义. 解对初值的连续依赖性说明，初值),(00y x 无法准确得到，但是我们能得到测量数据),(*0*0y x ，只要误差比较小，即 δδ<-<-*00*00,y y x x .我们就可以用),(*0*0y x 代替),(00y x 去计算，得到初值问题的解),,(*0*0y x x y ϕ=，这个解可以非常接近真实解),,(00y x x y ϕ=，即εϕϕ<-),,(),,(*0*000y x x y x x .同理，如果方程的右端函数),(y x f 不能准确得到，只能得到),(y x f 的近似函数),(~y x f ，即)),((,),(),(~D y x y x f y x f ∈<-δ我们就可以用),(~y x f 代替),(y x f 去计算，得到初值问题⎪⎩⎪⎨⎧==00~)(),(y x y y x f dx dy的解),,(00~y x x y ϕ=，那么),,(00~y x x y ϕ=能否代替),,(00y x x y ϕ=呢？我们有下面的解的连续依赖性定理.定理 （解对被积函数的连续依赖性定理）在区域D 上，),(y x f 和),(~y x f 都连续，而且关于变量y 满足李普希兹条件， 若初值问题⎪⎩⎪⎨⎧==00~)(),(y x y y x f dx dy在b x a ≤≤上有解),,(00~y x x y ϕ=，则对任意给定的正数ε，存在0>δ，只要),(y x f 满足)),((,),(),(~D y x y x f y x f ∈<-δ则初值问题（）的解),,(00y x x y ϕ=在b x a ≤≤上存在，且有εϕϕ<-),,(),,(00~00y x x y x x .证明 由解的存在唯一性定理知，初值问题⎪⎩⎪⎨⎧==00~)(),(y x y y x f dx dy的解),,(00~y x x y ϕ=存在，设其存在区间为],[b a ，且有⎰+=xx d y x f y y x x 0))],,(,([),,(00~~000~ξξϕξϕ而初值问题⎪⎩⎪⎨⎧==00)(),(y x y y x f dx dy的解),,(00y x x y ϕ=也存在，且可以表示为⎰+=xx d y x f y y x x 0))],,(,([),,(00000ξξϕξϕ则⎰⎰-=-xx xx d y x f d y x f y x x y x x 0))],,(,([))],,(,([),,(),,(0000~~0000~ξξϕξξξϕξϕϕ从而有⎰-≤-xx d y x f y x f y x x y x x 0|)),,(,()),,(,(|),,(),,(0000~~0000~ξξϕξξϕξϕϕ⎰-+-=xx d y x f y x f y x f y x f 0|)),,(,()),,(,()),,(,()),,(,(|0000~00~00~~ξξϕξξϕξξϕξξϕξ⎰-+-≤xx d y x f y x f y x f y x f 0|)),,(,()),,(,(||)),,(,()),,(,(|0000~00~00~~ξξϕξξϕξξϕξξϕξ⎰+-≤xx d y x y x N 0)|),,(),,((|0000~ξδξϕξϕ⎰-+-≤xx d y x y x N a b 0|),,(),,(|)(0000~ξξϕξϕδ由贝尔曼引理，得)(0000~)(),,(),,(a b N e a b y x x y x x --≤-δϕϕ取)(a b N e ab ---<εδ，则εϕϕ<-),,(),,(0000~y x x y x x .且解),,(00y x x y ϕ=在b x a ≤≤上存在.例1 考虑方程0,ln ,0≠=⎩⎨⎧-=y y y y dx dy解的情况.解 显然1,1,0-===y y y 是方程的解，当1,1,0-≠≠≠y y y 时，有y y dx dyln -=这时解得上半平面的通解为xCe ey -=，下半平面的通解为xCe ey --=.可以看到，对于Ox 轴上的初值)0,(0x ，在任意有限闭区间上解对初值连续依赖，但是，在),0[+∞上，无论),(00y x ，00≠y 如何接近)0,(0x ，只要x 充分大，过),(00y x 的积分曲线就不能与过)0,(0x 的积分曲线（即0=y ）任意接近了.这个例子说明，解在有限闭区间上对初值连续依赖，不能推广到无限区间，即，在无限区间上解对初值的连续依赖定理就不成立了. 我们有时不仅要求解对初值连续依赖，而且还要知道解),,(00y x x y ϕ=对初值00,y x 的偏导数0,y x ∂∂∂∂ϕϕ是否存在.下面给出解对初值的可微性定理.定理 （解对初值的可微性定理）如果函数),(y x f 以及),(y x f y '在区域D 内连续，则初值问题⎪⎩⎪⎨⎧==00)(),(y x y y x f dx dy的解),,(00y x x y ϕ=在它有定义的区间上有连续偏导数0,y x ∂∂∂∂ϕϕ.并且有⎰-=∂∂'xx y d y x f ey x f x y x x 000)),,(,(00000),(),,(ττϕτϕ及⎰=∂∂'xx y d y x f ey y x x 000)),,(,(000),,(ττϕτϕ习 题1.若函数),(y x f ，),(y x R 在区域D 内连续且满足李普希兹条件，设初值问题⎪⎩⎪⎨⎧=+=*0*0)(),(),(y x y y x R y x f dxdy的解为),,(*0*0~y x x y ϕ=，存在区间为],[b a .对任意的正数ε，存在0>δ，使对于满足)),((,),(D y x y x R ∈<δ的),(y x R ，以及满足δδ<-<-*00*00,y y x x的任意),(00y x ，初值问题⎪⎩⎪⎨⎧==00)(),(y x y y x f dx dy的解),,(00y x x y ϕ=也在区间],[b a 上存在，且有εϕϕ<-),,(),,(*0*0~00y x x y x x2.已知方程)sin(xy dxdy= 试求000000),,(==⎥⎦⎤⎢⎣⎡∂∂y x x y x x y 和0000000),,(==⎥⎦⎤⎢⎣⎡∂∂y x y y x x y3.设),,(00y x x ϕ是初值问题⎪⎩⎪⎨⎧==00)(),(y x y y x f dxdy的解，试证明0),(),,(),,(00000000=∂∂+∂∂y x f y y x x x y x x ϕϕ欧拉折线法在第一章，我们介绍了方程的初等解法，即用微积分的知识求得常微分方程的函数解.但是绝大多数的方程不能用初等方法求解，在第二章的前三节中，我们给出了柯西初值问题⎪⎩⎪⎨⎧==00)(),(y x y y x f dxdy在什么样的条件下，解存在且唯一；在什么条件下，解的存在区间可以延展；在什么条件下连续依赖于初值；在什么条件下，解对初值是可微的.有了这些准备，我们就可以研究柯西初值问题的近似解.下面我们介绍求近似解的方法，欧拉折线法.假定函数),(y x f 在区域：+∞<<-∞≤≤y b x a ,上连续，且关于变量y 满足李普希兹条件，求柯西初值问题⎪⎩⎪⎨⎧==00)(),(y x y y x f dxdy在区间],[0b x 上的近似解，我们采用的方法是：（1）等分区间],[0b x ，分点为nk kh x x k ,,1,0,0 =+=；小区间长度n x b h 0-=，（2）第一个小区间上用切线段逼近曲线：))(,(0000x x y x f y y -+=，（3）求出1x 所对应的纵坐标))(,(010001x x y x f y y -+=，（4）依次重复（2），（3）得到每个小区间上的线段，从而得到欧拉折线. 这样，我们就用欧拉折线作为柯西初值问题在区间],[0b x 近似解.欧拉折线法的前提是：柯西初值问题的解存在且唯一，而且解的存在区间是],[0b x .例1试用欧拉折线法，取步长1.0=h ，求初值问题⎪⎩⎪⎨⎧=+=1)1(22y y x dxdy的解在4.1=x 时的近似值.解 令22),(y x y x f +=，2)1,1(=f ，这时12-=x y ，代入1.11=x 得：2.11=y ，65.2)2.1,1.1(=f ，这时2.1)1.1(65.2+-=x y ， 代入2.12=x 得：465.12=y ，586225.3)465.1,2.1(=f ， 这时465.1)2.1(586225.3+-=x y ，代入3.13=x 得：8236225.13=y ，0155990225.5)8236225.1,3.1(=f ， 这时8236225.1)3.1(0155990225.5+-=x y ，代入4.14=x 得：53251824022.24=y习 题1. 试用欧拉折线法，取步长1.0=h ，求初值问题⎪⎩⎪⎨⎧=-=0)0(22y y x dxdy 的解在5.1=x 时的近似值.2．试用欧拉折线法，取步长1.0=h ，求初值问题⎪⎩⎪⎨⎧=+=2)1(22y y x dxdy 在区间]4.1,1[上的近似解.。