某一地区的病菌传染,三种人员的人数的状态方程,即可能受传染的人数x1,已被传染的病的人数x2,及已治愈的人数x3,并假设他们是因接触而传染。
设:α是x1中单位时间内的传染系数,β是x2中单位时间内被治愈的比例系数,则得以下状态方程:
f`(x1)=-α*x1*x2;
f`(x2)=α*x1*x2-β*x2;
f`(x3)=β*x2;
其中f`(x)是加速度
本题中,α=0.01,β=0.027
初始条件:t=0, x1=620人, x2=10人, x3=70人
选步长:h=0.025, 用四阶RK法模拟。
-------------程序
清单
安全隐患排查清单下载最新工程量清单计量规则下载程序清单下载家私清单下载送货清单下载
-------------------
#include
#include
#define A 0.01
#define B 0.027
void myFun(double *x,double *f){
f[0]=-A*x[0]*x[1];
f[1]=A*x[0]*x[1]-B*x[1];
f[2]=B*x[1] ;
}
int rht4(double *rs,double h,double x[],int n,int Num){
int j,i,ii;
double ps[4]={0,0.5,0.5,1};
double bs[4]={1.0/6,1.0/3,1.0/3,1.0/6};
double K[3];
double y[3];
double f[3];
for(i=0;i
本文档为【c求解四阶龙格库塔法】,请使用软件OFFICE或WPS软件打开。作品中的文字与图均可以修改和编辑,
图片更改请在作品中右键图片并更换,文字修改请直接点击文字进行修改,也可以新增和删除文档中的内容。
该文档来自用户分享,如有侵权行为请发邮件ishare@vip.sina.com联系网站客服,我们会及时删除。
[版权声明] 本站所有资料为用户分享产生,若发现您的权利被侵害,请联系客服邮件isharekefu@iask.cn,我们尽快处理。
本作品所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用。
网站提供的党政主题相关内容(国旗、国徽、党徽..)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。