FORTRAN+计算物理学学习日记(8)
第三章 边值问题和本征值问题
3.1numerov算法
例题
program mainimplicit noneinteger(8)::i,nreal(8)::x(600),y(600),a,h,pih=0.01n=500pi=3.14159263open(unit=100,file="xable.csv")open(unit=101,file="yable.csv")do i=1,n+1,1x(i)=(i-1)*hwrite(100,*)x(i)end doclose(100)y(1)=1write(101,*)y(1)y(2)=1-2*pi*pi*h*hwrite(101,*)y(2)a=(3-5*pi*pi*h*h)/(3+pi*pi*h*h)do i=2,n,1y(i+1)=2*a*y(i)-y(i-1)write(101,*)y(i+1)end doclose(101)
end program
绘制图像如下:
FORTRAN+计算物理学学习日记(8)相关推荐
- FORTRAN+计算物理学学习日记(1)
第一周:结合李录的计算物理学学习FORTRAN语言,这周的任务是插值函数,大致编写了四个小时,编写了一个双层循环的插值函数,进行了六次插值计算例题. 例题如下: 编写代码如下: program mai ...
- FORTRAN+计算物理学学习日记(2)
利用Fortran编写数值微分函数,利用向前向后两点式以及五点式求微分. 以上节例题为模板,编写向前向后两点式如下: !在例一的条件上,计算两点一次插值的微商,编写向前两点式,向后两点式 progra ...
- FORTRAN+计算物理学学习日记(5)
2.1常微分方程的简单数值解法 本节编写了四种简单的数值方法去求解常微分方程的初始问题,包括 Euler 方法.Taylor 级数 法.后向 Euler 方法和梯形公式. 注意点:定义格式或者数组定义 ...
- FORTRAN+计算物理学学习日记(6)
2.2Runge-Kutta 方法求解常微分方程 "直接利用 Taylor 级数展开提高算法的阶数有许多困难,特别是要确定函 数 f (x, y) 的导数,这在数值计算中是非常不方便的.为了 ...
- FORTRAN+计算物理学学习日记(7)
2.3多步法求解常微分方程 program mainimplicit noneinteger(8)::n,m,i,kreal(8),external::fa,fbreal(8)::h,x(200),y ...
- FORTRAN+计算物理学学习日记(4)
1.5基本数学运算中的求根 方法一:区间对分法求根 书中例题及编写代码如下 !!利用区间对分法求根 program mainimplicit nonereal(8)::a,b,x,t,ya=2b=3x ...
- Java学习日记-Day01
Java学习日记-Day01 Java语言概述 比特(byte)与字节 内存 Java基础知识图解 人机交互方式 常用的DOS命令 常用快捷键 计算机编程语言介绍 第一代语言 第二代语言 第三代语言 ...
- 深度学习日记 2 - 概率论与信息论基础
深度学习日记 2 - 概率论与信息论基础: 1.随机变量(random variable):是可以随机地取不同值的变量.我们通常用打印机 体的小写字母来表示随机变量本身,而用脚本字体中的小写字母来表示 ...
- GPU(CUDA)学习日记(十一)------ 深入理解CUDA线程层次以及关于设置线程数的思考
GPU(CUDA)学习日记(十一)------ 深入理解CUDA线程层次以及关于设置线程数的思考 标签: cuda存储线程结构网格 2012-12-07 16:30 6298人阅读 评论(4)收藏 举 ...
最新文章
- 深度神经网络中Inception-ResNet模块介绍
- 在浏览器上也能训练神经网络?TensorFlow.js带你玩游戏~
- 理解JVM如何使用Windows和Linux上的本机内存
- 开源 免费 java CMS - FreeCMS1.9 移动APP管理 执行配置
- php webview,Android:控件WebView显示网页 – tinyphp – 博客园
- bios升级工具_雨林木风U盘启动盘装系统制作工具再次升级。。。
- JAVA基础知识总结17(网络编程)
- 基于jsp网上书店设计与实现毕业设计论文
- CAD分图后批量重命名图纸名称
- java libtorrent_基于libtorrent最简单的BT下载程序 | 学步园
- 敏捷开发 建立愿景、使命_使用愿景板创建敏捷产品策略的10条技巧
- 旅游公寓APP开发特点
- 银行外包996,骑驴找马要跳槽
- GitCode 专属开源加速器,帮你实现开源社区运营
- 怎么看股票走势图(转贴)
- 第51次全国计算机等级考试,第51次全国计算机等级考试(NCRE)在我校顺利举行...
- 克莱因瓶计算机模拟,世界未解之谜:神奇的克莱因瓶永远装不满水,它或存在于四维空间?...
- 由Debuggerd导致的Android系统死机问题分析
- 登录亚马逊显示内部服务器出错,亚马逊后台登录提示密码错误有哪些原因?哪些操作会导致账号封号 ......
- 二分法和简单迭代法的优缺点_二分法和牛顿迭代法求解方程的比较.doc
热门文章
- js 遍历formData对象数据
- java流程控制,数组,集合
- lol手游s2服务器维护时间,lol手游s2赛季什么时候结束 赛季结束时间一览
- 150℃下使用的无铅焊锡:为车载基板开辟新路(一)
- 新版发布:大数据应用开发管理集成框架DataSphere Studio 1.1.0
- vue2.x项目中使用$once清除定时器,理解$once
- JavaScript排序算法:归并排序
- php过滤微信昵称的emoj,移除微信昵称中的emoji字符
- 中国电信天翼U盾产品荣获第三届网络安全国家标准优秀应用案例二等奖
- 开放网络发展历史大事记