数值分析欧拉法c语言程序

欧拉法c 语言程序

#include "stdio.h"

#include "math.h"

float f(a,b)

float a,b;

{

return (b-(2*a)/b);

}

main()

{ float x0,y0,x1,y1,h,yc,yp;

int N;

int n=1;

printf("please input x0,y0,h,N:");

scanf("%f,%f,%f,%d",&x0,&y0,&h,&N);

sign:x1=x0+h;

yp=y0+h*f(x0,y0);

yc=y0+h*f(x1,yp);

y1=(yp+yc)/2;

while(n!=N)

{printf("%8.7f, %8.7f, %d\n",x0,y0,n); n=n+1;

x0=x1;

y0=y1;

goto sign;

}

printf("output\nx1=%10.9f\ny1=%10.9f",x1,y1); getch();

}

欧拉法c 语言程序

#include "stdio.h"

#include "math.h"

float f(a,b)

float a,b;

{

return (b-(2*a)/b);

}

main()

{ float x0,y0,x1,y1,h,yc,yp;

int N;

int n=1;

printf("please input x0,y0,h,N:");

scanf("%f,%f,%f,%d",&x0,&y0,&h,&N);

sign:x1=x0+h;

yp=y0+h*f(x0,y0);

yc=y0+h*f(x1,yp);

y1=(yp+yc)/2;

while(n!=N)

{printf("%8.7f, %8.7f, %d\n",x0,y0,n); n=n+1;

x0=x1;

y0=y1;

goto sign;

}

printf("output\nx1=%10.9f\ny1=%10.9f",x1,y1); getch();

}


相关文章

  • 数值分析习题大作业

    • 一. [设计题一] .............................................................................................................. ...查看


  • 欧拉法与龙格库塔法比较分析

    • 解微分方程的欧拉法,龙格-库塔法简单实例比较 欧拉方法(Euler method)用以对给定初值的常微分方程(即初值问题)求解分为前EULER法.后退EULER法.改进的EULER法. 缺点: 欧拉法简单地取切线的端点作为下一步的起点进行计 ...查看


  • 数值分析Matlab作业龙格库塔欧拉方法解二阶微分方程

    • Matlab 应用 使用Euler和Rungkutta方法解臂状摆的能量方程 背景 单摆是常见的物理模型,为了得到摆角θ的关于时间的函数,来描述单 摆运动.由角动量定理我们知道 M=Jε 化简得到 dθg +sinθ=02 dtl 在一般的 ...查看


  • 多种微分方程数值计算方法分析

    • 2010年8月 第4期 城 市 勘 测 Urban Geotechnical Investigation& Surveying Aug. 2010 No. 4 文章编号:1672-8262(2010)04-117-03中图分类号:O ...查看


  • CAE有限元分析软件的比较

    • CAE有限元分析软件的比较 随着现代科学技术的发展,人们正在不断建造更为快速的交通工具.更大规模的建筑物.更大跨度的桥梁.更大功率的发电机组和更为精密的机械设备.这一切都要求工程师在设计阶段就能精确地预测出产品和工程的技术性能,需要对结构的 ...查看


  • 微分方程求解数值方法

    • 重庆三峡学院数值计算方法实验课测验 微分方程求解的后退欧拉法(后退梯形法).龙格库塔 法(三阶.四阶公式) 院 系 数学与统计学院 专 业 信息与计算科学 姓 名 庞欢 年 级 2009级 学 号 [1**********]2 指导教师 刘 ...查看


  • 三维弹塑性流体力学欧拉方法动态并行研究

    • 第四届全国计算爆炸力学会议论文集 三维弹塑性流体力学欧拉方法动态并行研究 李征 廖振民 梁仙红 刘军 冯其京 何长江 (北京应用物理与计算数学研究所,北京 100088) 摘 要:欧拉方法是流体力学数值模拟的主要方法,适于计算多介质大变形的 ...查看


  • 电力系统稳定计算

    • 目录 一.需求分析 . ............................................................................................................ ...查看


  • 实验四 常微分方程初值问题数值解法

    • 实 验 报 告 课程名称 数值分析 实验项目 常微分方程问题初值问题数值解法 一. 实验目的 1.理解如何在计算机上实现用Euler 法.改进Euler 法.Runge -Kutta 算法求一阶常微分方程初值问题 ⎧y '(x ) =f ( ...查看


热门内容