C语言-预估校正法求常微分方程

最新推荐文章于 2023-08-27 05:47:37 发布

weixin_33713350

最新推荐文章于 2023-08-27 05:47:37 发布

阅读量935

点赞数

CC 4.0 BY-SA版权

文章标签： c/c++

原文链接：http://www.cnblogs.com/zhubinglong/p/6141934.html

本文提供了一个使用C语言实现的数值解法示例，通过欧拉法和改进的欧拉法求解一阶微分方程组。该程序定义了方程组的初值条件，并迭代计算了一系列近似解。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

#include<stdio.h>
#include<math.h>
#define n 14
int main(){
    double a = 0.0, b = 1.4,h,m=0.0;
    //int n = 14;

    double x[n + 1] = { 0.0 }, y[n + 1] = { 0 };
    h = (b - a) / n;
    x[0] = a;
    y[0] = m;
    for (int i = 1; i <= n; i++)
    {
        y[i] = y[i - 1] + h*(1 + y[i - 1] * y[i - 1]);
        x[i] = a + i*h;
        y[i] = y[i - 1] + h*0.5*((1 + y[i - 1] * y[i - 1]) + (1+y[i]*y[i]));
        printf("x[%d]=%e , y[%d]=%e\n",i, x[i], i, y[i]);
    }
    return 0;
}