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