四阶龙格库塔法程序设计例:以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转载请标明出处.