算法 (1)输入a,b,f(x,y)2N,y(2)置h= b-a 0=a2n=1 (3)计算 =f(rm-l, K,= hf K,=hf( h K K3=hf(e Ka=hf(n-th,y,-+K3) 2 Un=n-1+-(K+2K2+2K,+K4 x=a+nh (4)输出(xn,yn) (5)若n<3,置n+1→n,返回3; 否则,置n+1→n,0→p0,0→c0,转6。0 0 1 1 1 1 1 1 2 1 1 1 , , ( , ), , (2) , , 1 (3) ( , ) ( , ) 2 2 n n n n n n b a a b f x y N y h x a n N f f x y h K K hf K hf x y K − − − − − − − = = = = = = + + 算法 （ ）输入 置 计算 2 3 1 1 4 1 1 3 1 1 2 3 4 0 0 ( , ) ( , ) 2 2 1 ( 2 2 ) 6 (4) ( , ) (5) 3, 1 , 3 1 , 0 , 0 , 6 n n n n n n n n n h K hf x y K hf x h y K y y K K K K x a nh x y n n n n n p c − − − − − = + + = + + = + + + + = +  +  +    输出 若 置 返回 ； 否则，置 转