方程的改进欧拉公式_第1页
方程的改进欧拉公式_第2页
方程的改进欧拉公式_第3页
方程的改进欧拉公式_第4页
方程的改进欧拉公式_第5页
已阅读5页，还剩3页未读， 继续免费阅读

下载本文档

版权说明：本文档由用户提供并上传，收益归属内容提供方，若内容存在侵权，请进行举报或认领

文档简介

1、2013-2014（1）专业课程实践论文题目：方程的改进欧拉公式 一、算法理论按照计算的数值解时候，如果每次只迭代一次，相当于将欧拉公式与梯形格式结合使用：先用欧拉公式求出，称为预测值，然后用梯形格式校正求得近似值，即预测： ；校正： 将上式称为由欧拉公式和梯形格式得到的预测校正系统，也叫改进的欧拉公式。它是显式的单步法。为了便于编制程序上机，上式常改写成 改进的欧拉方法算法如下：（1） 输入区间等分数N，初始值。（2） 输出在的N个点处的近似值。（3） 置。（4） 计算，置输出。（5） 若，置，转（4）；否则，停机.计算框图见二。 二、算法框图 三、算法程序#include <stdi

2、o.h> #include <math.h> #include <stdlib.h> double fun(double x,double y) return x*y; void Euler(double (*fun)(double ,double), double h,double x0,double y0 ,double low,double up) int n =0,i=0,j=0; double x1 = x0,y1=y0; n = (up-low)/h; printf("第n个变量 变量的值xn 函数值ynn"); printf(&q

3、uot;%lf %lf n",x0,y0); for ( i = 1 ; i<=n ; i+) x1 = x0 + i*h; y1 = y1+ h*fun(x1,y1); printf("%d %lf %lfn",i,x1,y1); void main() double h,x0,y0,low,up; printf("*输入初始值x0="); scanf("%lf",&x0); printf("*输入初始值y0="); scanf("%lf",&y0); prin

4、tf("*请输入利用“改进的欧拉公式”使用变量的范围n"); printf("*输入区间下界low="); scanf("%lf",&low); printf("*输入区间上界up="); scanf("%lf",&up); printf("*输入步长h="); scanf("%lf",&h); Euler(fun,h,x0,y0,low,up);return 0; 四、算法实现例1用改进的欧拉公式，求解常微分方程初值问题的解。，解：（1）输入的数值； （2）输入初始值的数值； （3）输入利用改进的欧拉公式使用的变量下界的数值； （4）输入利用改进的欧拉公式使用的变量上界的数值； （5）输入步长的数值； （6）得到结果。 例2. 用改进的欧拉公式，求解常微分方程的初值问题的解。，解：（1）输入的数值； （2）输入初始值的数值； （

人人文库> 全部分类> 教育资料 > 课件下载

温馨提示

  • 1. 本站所有资源如无特殊说明，都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等，如果需要附件，请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸，网页内容里面会有图纸预览，若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间，仅对用户上传内容的表现方式做保护处理，对用户上传分享的文档内容本身不做任何修改或编辑，并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容，请与我们联系，我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

最新文档

评论

0/150

提交评论