爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

搜索
查看: 3534|回复: 2

[混合编程] IDL编程学习之五点“牛顿-科茨公式”积分

[复制链接]
发表于 2017-8-6 11:43:25 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册

x
题目:
求积分,积分限:[0,0.8],
函数:y=0.2+25.0*x-200.0*x^2.0+675.0*x^3-900.0*x^4+400.0*x^5.0
现在,这个积分以数值形式给出,
X = [0.0, .12, .22, .32, .36, .40, .44, .54, .64, .70, .80]
y = [0.200000, 1.30973, 1.30524, 1.74339, 2.07490, 2.45600, 2.84299, 3.50730, 3.18194, 2.36302, 0.231964]

我们可以先计算下这个积分的精确结果:
x=[0.0,0.8]
dy=0.2*x+25.0*x^2.0/2-200.0*x^3.0/3+675.0*x^4.0/4-900.0*x^5.0/5+400.0*x^6.0/6;(自己写的不定积分)
print,dy

输出:     0.000000      1.64053
可见,精确结果是1.64053


“五点牛顿-科茨公式”积分:
  X = [0.0, .12, .22, .32, .36, .40, .44, .54, .64, .70, .80]
  Y = [0.200000, 1.30973, 1.30524, 1.74339, 2.07490, 2.45600, $
    2.84299, 3.50730, 3.18194, 2.36302, 0.231964]
  print,INT_TABULATED(X, Y)

输出:1.62323


数值积分一向不准
密码修改失败请联系微信:mofangbao
发表于 2017-8-6 14:06:06 | 显示全部楼层
建议自己学的差不多,做点东西,把成果在传上来,大家一起学习,虽然你的学习精神是好的,但是天天这样就是在打扰大家,谢谢,个人观点。
密码修改失败请联系微信:mofangbao
 楼主| 发表于 2017-8-7 08:33:14 | 显示全部楼层
墨墨卡托 发表于 2017-8-6 14:06
建议自己学的差不多,做点东西,把成果在传上来,大家一起学习,虽然你的学习精神是好的,但是天天这样就是 ...

你看IDL的帖子是多是少?根本没人关心好吧
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表