二阶龙格-库塔(共5页)_第1页
二阶龙格-库塔(共5页)_第2页
二阶龙格-库塔(共5页)_第3页
二阶龙格-库塔(共5页)_第4页
二阶龙格-库塔(共5页)_第5页
全文预览已结束

下载本文档

版权说明：本文档由用户提供并上传，收益归属内容提供方，若内容存在侵权，请进行举报或认领

文档简介

1、精选优质文档-倾情为你奉上#include<stdlib.h>#include<stdio.h> /*n表示几等分，n+1表示他输出的个数*/ int RungeKutta(double y0,double a,double b,int n,double *x,double *y,double (*function)(double,double) double h=(b-a)/n,k1,k2;int i; x0=a; y0=y0; for(i=0;i<n;i+) xi+1=xi+h; k1=function(xi,yi); k2=function(xi+h/2,y

2、i+h*k1/2);yi+1=yi+h*k2; return 1; double function(double x,double y) return y-2*x/y; int main() int i; double x6,y6; printf("用二阶龙格-库塔方法n"); RungeKutta(1,0,1,5,x,y,function); for( i=0;i<6;i+) printf("x%d=%f,y%d=%fn",i,xi,i,yi); return 1; 另一个#include<stdio.h>#include<ma

3、th.h> double fxy(double xi,double yi) /*定义函数fxy*/ double y; y=yi-2*xi/yi; return(y); void main() double x0,y0,h,xi,yi,yi_1,xk2,yk2,xk3,yk3,xk4,yk4,k1,k2,k3,k4; int i; x0=0; /*赋初始值*/ y0=1; h=0.1; xi=x0; yi=y0; for(i=1;i<=10;i+) /*循环开始*/ k1=h*fxy(xi,yi); /*求解k1值*/ xk2=xi+0.5*h; /*求解k2的值*/ yk2=yi

4、+0.5*k1; k2=h*fxy(xk2,yk2); xk3=xi+0.5*h; /*求解k3的值*/ yk3=yi+0.5*k2; k3=h*fxy(xk3,yk3); xk4=xi+h; /*求解k4的值*/ yk4=yi+k3; k4=h*fxy(xk4,yk4); yi_1=yi+(k1+2*k2+2*k3+k4)/6; /*求解yi+1的值*/ yi=yi_1; if(i=1) printf("输出函数yi的近似值:n"); /*输出所有的yi值*/ printf("y0 =%.10f",y0); printf("y%d =%.10f",i,yi); if(i+1)%1=0) /*每

人人文库> 全部分类> 教育资料 > 辅导培训

温馨提示

  • 1. 本站所有资源如无特殊说明，都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等，如果需要附件，请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸，网页内容里面会有图纸预览，若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间，仅对用户上传内容的表现方式做保护处理，对用户上传分享的文档内容本身不做任何修改或编辑，并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容，请与我们联系，我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

最新文档

评论

0/150

提交评论