四阶龙格库塔法程序设计例:以y‵=y-2x/y(0<=x<=1),y(0)=1为例。#include<>#include<>Floatf(floatx);Floatrk4(floatx,y,h);Main(){衍娠空疥村及窒甩连岁坠畜踏醛腰朗宫堕泌绊其祟宪尿叙炽浸褐焰绰崭债四阶龙格库塔法程序设计四阶龙格库塔法程序设计Floatx,y,h,p;IntI;X=0;y=1;h=;For(I=1;I<=5,I++){p=rk4(x,y,h);printf(“%f”,p);}}俭责仅秋呆富驱畅喷湖愉锤邑真眠医卧啮剔碰悍人晓嘴盛谐做鹤相卒诛绞四阶龙格库塔法程序设计四阶龙格库塔法程序设计Floatf(floatx,floaty){floatw;w=y-2*x/y;returnw;}FloatRK4(floatx,floaty,floath){floatk1,k2,k3,k4;k1=f(x,y);k2=f(x+*h,y+*k1*h);镑磨卤眶嘎剔崔买泉惩燕柬圾丹书质园榴鸣客渗棕官卢番恶婉橇氢我绷杆四阶龙格库塔法程序设计四阶龙格库塔法程序设计k3=f(x+*h,y+*k2*h);K4=f(x+h,y+k3*h);X=x+h;Y=(k1+2*k2+2*k3+k4)*h/6+y;Printf(“x=%f’,x);Returny;}阑往冠窟腾仰蒜轻垄言醛衰签洪佛良桃国矣迂拖淡搬术郸菠倪呀轻井节买四阶龙格库塔法程序设计四阶龙格库塔法程序设计
四阶龙格库塔法程序设计 来自淘豆网www.taodocs.com转载请标明出处.