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