原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
主程序外部例程
工程分析程序设计 上机作业（五）
模块化编程
上机目的：掌握内部例程、主程序、外部例程、模块等功能的使用方法。
有一个六边形，求其面积。为求面积，作了3条辅助线。如图所示：（提示，三角形面积=，其中，a、b、c为三个边长）。要求用内部函数来计算每个三角形的面积。
编写一个子例程子程序，SUM(S，T，N1，N2)。把整型数N1到N2进行求和，并把求和的结果放置到S，把整型数N1到N2进行求积，并把求积的结果放置到T。并用这个子程序来计算：
编写函数子程序GDC求两个数的最大公约数。求最大公约数的算法如下：把两个数中大的那个数作为被除数，两数相除得到一个余数。把余数去除除数得到新一轮的余数。不断重复这一过程直到余数为0，这时的除数就是两个数的最大公约数。
调用此函数，求1260，198，72三个数的最大公约数。
编写一个模块程序，提供以下服务：定义出常量?、e。定义出子程序，实现求和、求阶乘n!。并在主程序中计算如下结果：从键盘上输入整数n，求（实型）。
编写外部函数，计算的值，直到最后一项的绝对值小于10-6为止。再编写主程序，从键盘读入x，调用该外部函数并输出sinx的计算结果。注意不能用Fortran的标准函数SIN(X)。
编写程序，用欧拉法数值求解一阶微分方程。
数值解法的基本思想：用差分方程代替微分方程，然后在若干个离散点上逐点求解差分方程，得到各离散点x0, x1, x2, … 等处函数的近似值y0, y1, y2, … 。其中各离散点x0, x1, x2, …之间的距离称为步长h。
欧拉法的基本思想是，一阶微分可用向前差分代替，即。带入微分方程，可得
因此，，其中。
当给定x0, y0和步长h后，即可按下面步骤求得数值解：
请按上述方法编写程序用Euler法求解微分方程，当x0=0时，y0=1.0。取h=0.1，试求出当xi=0.1,0.2,0.3,0.4….,1.0时的yi值。
文档评论（0）