f(x,y)a的计算:假设y(xk-1)=y已知 k-1 f(x, y)dx a hf(k-1, Dk-1) 矩形求积公式 k-1 f(,y)≈[/(x-,y-)+f(x,y(x h 梯形求积公式误差为R(Tk) f"(7k,y(7k)) 12 f(x, ydx lf(rk-,Dk-1+4f(k1+o,y( O+f(rk, y(kI h Simpson求积公式误差为R(S) f()(nk,y(nk) 180×16 将以上求积公式代入(1)试式,并加以处理就 可得到相对应的求解公式( , ) : 1 ò 的计算 - k k x x f x y dx ( , ) ( , ) 1 1 1 ò » - - - k k x x f x y dx hf x y k k 矩形求积公式 [ ( , ) ( , ( ))] 2 ( , ) 1 1 1 k k k k x x f x y f x y x h f x y dx k k » - - + ò - 假设y(xk -1 ) = yk -1已知 梯形求积公式,误差为 )) ( , ( ))] 2 , ( 2 [ ( , ) 4 ( 6 ( , ) 1 1 1 1 1 k k k k k k x x f x y x h y x h f x y f x h f x y dx k k » - - + - + - + + ò - Simpson求积公式,误差为 将以上求积公式代入(11)式,并加以处理就 可得到相对应的求解公式 ( , ( )) 12 ( ) 3 k k k f y h R T = - ¢¢ h h ( , ( )) 180 16 ( ) ( 4) 5 k k k f y h R S h h ´ = -