C语言实现亚当姆斯方法

#include<stdio.h>
#include<stdlib.h>
double fx(double x, double y)
{
return(y - 2 * x / y);
}
int main() {

double x1, K1, K2, K3, K4, y1;
double x0, y0, h,yx1,yx0,y2,y3,y4,x2,x3,x4,yp;
int N;
printf("请输入x0，y0，h，N:");
scanf_s("%lf%lf%lf%d", &x0, &y0, &h, &N);
yx0 = y0;
//printf("%lf%lf%lf%d", x0, y0, h, N);
for (int i = 0; i < 3; i++) {x1 = x0 + h;K1 = fx(x0, y0);K2 = fx(x0 + h / 2, y0 + h * K1 / 2);K3 = fx(x0 + h / 2, y0 + h * K2 / 2);K4 = fx(x1, y0 + h * K3);y1 = y0 + h * (K1 + 2 * K2 + 2 * K3 + K4) / 6;if (i ==0)yx1 = y1;if (i ==1)y2 = y1;if (i ==2)y3 = y1;printf("%lf\t%lf\n", x1, y1);x0 = x1;y0 = y1;
}
x0 = x0 - 3 * h;
x1 = x0 + h;
x2 = x0  + 2 * h;
x3 = x0  + 3 * h;
for (int n = 4; n <= N; n++)
{x4 = x3 + h;yp = y3 + h * (55 * fx(x3, y3) - 59 * fx(x2, y2) + 37 * fx(x1, yx1) - 9 * fx(x0, yx0)) / 24;y4 = y3 + h * (9 * fx(x4,yp) + 19 * fx(x3, y3) - 5 * fx(x2, y2) + fx(x1, yx1)) / 24;printf("%lf\t%lf\n", x4, y4);       x0 = x1, yx0 = yx1;x1 = x2, yx1 = y2;x2 = x3, y2 = y3;x3 = x4, y3 = y4;
}
system("pause");
return 0;

}// 为4阶龙格库塔修改得极为复杂仅供参考

C语言实现亚当姆斯方法相关推荐

matlab双端测距算法,一种使用亚当姆斯法的双端故障测距算法
高压输电线路是电力系统的重要组成部分,随着电力系统规模的扩展,高压远距离输电线路日益增多,而高压输电线路故障对电力系统.工农业生产和人们日常生活带来一定影响.高压输电线路的准确故障测距是从技术上保证电 ...
【C/C++】龙格库塔+亚当姆斯求解数值微分初值问题
1 /* 2 解数值微分初值问题: 3 龙格-库塔法求前k个初值 + 亚当姆斯法 4 */ 5 #include<bits/stdc++.h> 6 using namespace std; ...
【数值分析】python实现四阶亚当姆斯预测矫正方法
python代码 import numpy as np import matplotlib.pyplot as pltDT = np.array(input("份数:").spli ...
四阶亚当姆斯预报-校正系统
#include <iostream> #include <array> #include <vector> using namespace std;double ...
阿当姆斯matlab,数值计算方法与MATLAB应用
第1篇数值计算的基本方法和概念第1章算法与误差 1.1 算法 1.2 误差第2章方程求解 2.1 引言 2.2 二分法 2.3 迭代法 2.4 牛顿法 2.5 弦截法 2.6 解非线性方程组 ...
3.2.4乔姆斯基的语言观2
自然语言的使用需要说出或写下符号才算完成.人类婴幼儿对母语的学习,是从父母.周围其他人的讲话,不断接收实际例子,以熟悉不同场景下字.词.复合词.词组等素材的使用,以及各种组合方式的使用:然后自己不断实 ...
语言学中不同视角的抽象与具体——谈索绪尔、乔姆斯基、韩礼德的语言观
任何领域的学者或许都无法避免区分一种具体的现象和一种抽象的本质的问题,即现象-本质/具体-抽象的问题.这是因为研究者们似乎都不会仅仅满足于自己能观测到的东西--他们更想知道为什么自己会观测到这些东西. ...
C语言复习之判断Armstrong数(阿姆斯壮数)
首先我们来了解一下什么是Armstrong数(阿姆斯壮数) 在百度百科上阿姆斯壮数是这样说的:Armstrong 数,就是n位数的各位数的n次方之和等于该数比如: 153=1^3+5^3+3^3 1 ...
人类语言能力的自然演化：乔姆斯基对阵达尔文|《达尔文的危险思想》
导语查尔斯·达尔文用他的<物种起源>为生命的多姿多彩提供了一种解释:是演化和自然选择造就了这个星球上缤纷的生命.然而自诞生的那一天起,达尔文自然选择和演化的理论就引发过很多争议.很多人认 ...

C语言实现亚当姆斯方法

C语言实现亚当姆斯方法相关推荐

最新文章

热门文章