四阶龙格-库塔法解微分方程(C )一.作业:
用四阶龙格—库塔法求下列方程:
步长
。并与解析解
比较。
二.程序
#include
#include
#include
using namespace std;
double f(double x,double y){
return 10*x*(1-y);
}
int main(){
intn,i;
double h,k1,k2,k3,k4;
cout>n;
double *x=new double [n+1];
double *y=new double [n+1];
double *y1=...
一.作业:
用四阶龙格—库塔法求下列方程:
步长
。并与解析解
比较。
二.程序
#include
#include
#include
using namespace std;
double f(double x,double y){
return 10*x*(1-y);
}
int main(){
intn,i;
double h,k1,k2,k3,k4;
cout<<"Please input the number of intervals:";
cin>>n;
double *x=new double [n+1];
double *y=new double [n+1];
double *y1=new double [n+1];
h=1.0/n;
y[0]=0;
for(i=0;i<=n;i++){
x[i]=i*h;
k1=f(x[i],y[i]);
k2=f(x[i]+h/2,y[i]+k1*h/2);
k3=f(x[i]+h/2,y[i]+k2*h/2);
k4=f(x[i]+h,y[i]+k3*h);
y[i+1]=y[i]+(k1+2*k2+2*k3+k4)*h/6;
y1[i]=1-exp(-5*pow(x[i],2));
cout<
本文档为【四阶龙格-库塔法解微分方程(C )】,请使用软件OFFICE或WPS软件打开。作品中的文字与图均可以修改和编辑,
图片更改请在作品中右键图片并更换,文字修改请直接点击文字进行修改,也可以新增和删除文档中的内容。
该文档来自用户分享,如有侵权行为请发邮件ishare@vip.sina.com联系网站客服,我们会及时删除。
[版权声明] 本站所有资料为用户分享产生,若发现您的权利被侵害,请联系客服邮件isharekefu@iask.cn,我们尽快处理。
本作品所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用。
网站提供的党政主题相关内容(国旗、国徽、党徽..)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。