adams预测-校正系统

这个程序是大学三年级学习《数值分析》的时候写的程序,计算机以数学作为基础,很好的了解数有助于我们学习程序设计。

#include<iostream.h>
#include<iomanip.h>
#include<math.h>
#include<stdio.h>
#define M 80
double f(double,double);
//void runge4(double,double,double,double,double);
void runge4(double a,double b,double rk[],double h,double x[])
{
    int k,max;
 double k1,k2,k3,k4;
 max=(int)((b-a)/h);
 x[0]=a;
 for(k=1;k<=max;k++)//runge-kutta
 {
        x[k]=x[k-1]+h;
  k1=f(x[k-1],rk[k-1]);
  k2=f((x[k-1]+h*0.5),(rk[k-1]+h*0.5*k1));
  k3=f((x[k-1]+h*0.5),(rk[k-1]+h*0.5*k2));
  k4=f((x[k-1]+h),(rk[k-1]+h*k3));
  rk[k]=rk[k-1]+h*(k1+2*k2+2*k3+k4)/6;
 }
}
void main()
{
 double x[M]={0.0},a,b,h,y[M]={0.0},k[5],y0;
 int n,i,m;
 cout<<"

  • 0
    点赞
  • 3
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值