//龙格-库塔法
#include<stdio.h>
#include<math.h>
#define PI 3.1415
#define h 0.1
#define l -5
#define m 8
#define o -3
//指代f(t,u,v)函数
double funcf(double t,double u,double v){
return v;
}
//指代g(t,u,v)函数
double funcg(double t,double u,double v){
return -(l*u+m*u*u*u+o*u*u*u*u*u);
}
void Runge(){
int n=0;
double k1,k2,k3,k4,l1,l2,l3,l4;
double u[11],v[11],t[10],y[10];
u[0]=sqrt(1+1/sin(1)),v[0]=(-cos(1)/sin(1))*sqrt(1+1/sin(1));
printf("%10s %10s %10s %10s\n","t的值","近似解","精确解","误差");
for(n=0;n<3;n++){
//定义步长
t[n]=h*n;
//精确解
//y[n]=sin(PI*t[n]);
y[n]=sqrt(1+1/sin(2*t[n]+1));
/* printf("%10f\n",y[n]);
printf("%10f %10f %10f %10f\n",t[n],u[n],y[n],fabs(u[n]-y[n]));*/
k1=funcf(t[n],u[n],v[n]);
l1=funcg(t[n],u[n],v[n]);
k2=funcf(t[n]+h/2,u[n]+h*k1/2,v[n]+h*k1/2);
l2=funcg(t[n]+h/2,u[n]+h*l1/2,v[n]+h*l1/2);
k3=funcf(t[n]+h/2,u[n]+h*k2/2,v[n]+h*k2/2);
l3=funcg(t[n]+h/2,u[n]+h*l2/2,v[n]+h*l2/2);
k4=funcf(t[n]+h,u[n]+h*k3,v[n]+h*k3);
l4=funcg(t[n]+h,u[n]+h*l3,v[n]+h*l3);
printf("%10f %10f %10f %10f\n",l1,l2,l3,l4);
u[n+1]=u[n]+(h/6)*(k1+2*k2+2*k3+k4);
v[n+1]=v[n]+(h/6)*(l1+2*l2+2*l3+l4);
}
}
void main(){
Runge();
}
-
等级
本版专家分：0结帖率 50% -
展开阅读全文
-
等级
本版专家分：456858勋章
-
- 状元 2017年 总版技术专家分年内排行榜第一
-
- 榜眼 2014年 总版技术专家分年内排行榜第二
-
- 探花 2013年 总版技术专家分年内排行榜第三
-
-
进士
2018年总版新获得的技术专家分排名前十
2012年 总版技术专家分年内排行榜第七
-
-
代码功能归根结底不是别人帮自己看或讲解或注释出来的；而是被自己静下心来花足够长的时间和精力亲自动手单步或设断点或对执行到某步获得的中间结果显示或写到日志文件中一步一步分析出来的。
提醒：再牛×的老师也无法代替学生自己领悟和上厕所！
单步调试和设断点调试（VS IDE中编译连接通过以后，按F10或F11键单步执行，按Shift+F11退出当前函数；在某行按F9设断点后按F5执行停在该断点处。）是程序员必须掌握的技能之一。
-
等级
本版专家分：6469 -
调试就知道了。
- 《Python程序设计与算法基础教程（第二版）》江红 余青松 全部章节的课后习题，上机实践，课后答案，案例...
1. 这本书对Python的知识点的描述很详细，而且排版看的很舒服. 2. 几个例题：假装自己从零开始学，将一些有代表性、有意思的例题抽取出来. 3. 还有自己对一部分课后复习题，全部课后上机实践题的解题思路
- Python中for循环索引变量的作用域
作者: Eli Bendersky 原文链接：...我从一个小测试开始。这个函数做什么？ def foo(lst): a = 0 for i in lst: a += i b = ...
- python定义一个数据自动去重的函数_110道Python面试题（真题小结）
该文110道面试题全部来自于大家笔试面试时候拍照后发到群里求助的题目，并自己一道一道亲自做了，大部分题目属于巩固基本python知识点的题目，希望对基本知识不熟悉的同学，能认真做一遍，肯定会有不少收获1、一行...
- python程序在运行时、如果发生异常_110道python面试笔试题汇总，你能答对几道？...
该文110道面试题全部来自于大家笔试面试时候拍照后发到群里求助的题目，并自己一道一道亲自做了，大部分题目属于巩固基本python知识点的题目，希望对基本知识不熟悉的同学，能认真做一遍，肯定...
- 对Dataframe中异常值检测求助
对Dataframe中异常值检测求助背景主要需求目前随便用一组数的进展dataframe的进展 背景 工程中检测很多个阀门的开关，并且按照时间顺序对阀门进行排序。 刚用python，纯小白，还请大家不吝赐教~~ 主要需求 挑选出...
- java 读取 nextline_新手求助：.nextLine()可以读取下一行，那读取上一行呢?
很明显应该使用BufferedReader类是最好的."当.nextLine()读取到TMES的时候，让读取的行返回到00000000,4"BufferedReader有个readLine()方法,就是每次都读一行的意思(当然,首先在写readLine()方法之前要写上一个while...
- python123作业第四周_python3学习第四周
(1-8)上周知识知识回顾：1.文件操作文件操作步骤：a.打开文件，得到句柄。b.通过句柄操作文件c.关闭文件文件操作：a.增b.删c.改d...."r+"读追加模式打开文件。常用方法：seek() #位置的跳跃tell() #打印出位置truncate...
- 记一次生成一亿条测试数据
今天接到一个数据组同事的求助，需要至少一亿条数据来测试数据库性能问题。格式只要字段间有空格隔开即可，对文本格式没有要求。 简单思考了一下，Java生成这一亿条数据所准备的工作： 设计数据文本格式：clazz + ...
- Python3-Cookbook总结 - 第二章：字符串和文本
几乎所有有用的程序都会涉及到某些文本处理，不管是解析数据还是产生输出。 这一章将重点关注文本的操作处理，比如提取字符串，搜索，替换以及解析等。 大部分的问题都能简单的调用字符串的内建方法完成。 但是，...
- java较全的面试题
一、Java基础部分 1. JAVA的基本数据类型有哪些 ? String 是不是基本数据类型 ?...的区别. 5. 在JAVA中如何跳出当前的多重嵌套循环？ 6. switch语句能否作用在byte上，能否作用在long上，能否作用在Strin...
- 第 1 部分 从输出流中读取
现在，我们有多个平台和概要及许多新的 API，并且还在制作的差不多有数百种之多。尽管 Java 语言的复杂程度在不断增加，但它对于日常的编程任务而言仍是一个出色的工具。虽然有时您会陷入那些日复一日的编程问题中，
- 【求助】scala运行时报错ArrayIndexOutOfBoundsException
Scala程序在运行时出错 : Exception in thread “main” java.lang.ArrayIndexOutOfBoundsException: 10582 下面是pom文件: &lt;artifactId&gt;Spark&lt;/artifactId&gt; &...
- 20172328《程序设计与数据结构》第七周学习总结
学号20172328《程序设计与数据结构》第七周学习总结 第九章教材学习内容总结 1.继承的概念： 继承是在父类和子类之间建立一种“是”的关系。 用于派生新类的原始类称为父类、超类或基类，被派生出的类就称为子类...
- es7遍历mysql数据_es5、es6、es7知识点的归纳与总结（笔记）
es5、es61.supersuper关键字使得子类调用了父类中的构造函数,父类可以使用子类传进来的属性进行各种运算,super 关键字 还可以调用父类的方法class Father{constructor(x,y){this.x=x;this.y=y;}sum(){console.log...
- 这才是你需要的最基础的.Net基础面试题（通俗易懂，最基础的.Net）
由于小编能力有限，这里的解释大多都是小编自己百度然后自己理解的 可能有的地方会有错误欢迎评论私信我更改 一群热爱技术并...
- 精进：如何成为一个很厉害的人--作者：采铜
文章目录精进：如何成为一个很厉害的人序 用更勇敢的方式去生活01 时间之尺 我们应该怎样对待时间活在“全部的现在” 从当下出发，联结过去与未来那么，一个人应该怎样对待时间呢？不同场合，不同的时间视角由当下向...
- 26种数据挖掘软件比较及介绍
1、数据挖掘软件之ANGOSS软件公司 KnowledgeSTUDIO 平台：NT , Windows 9X 功能：CHAID算法，支持PMML，留有与SAS、S-Plus的接口，能够灵活地导入外部模型和产生规则，包含神经网络建模的能力 优点：响应快，...
- 自动化测试的7个步骤
自动化测试的7个步骤 一 改进软件测试过程 二 定义需求 三 验证概念 四 支持产品的可测试性 五 可延续性的设计（design for sustainability） 六 有计划的部署 七 面对成功的挑战 步骤一：改进软件测试过程 如果你...
- 第五章：并发及并行（Ⅲ）
并发及并行
- 第2章 包装外观（Wrapper Facade）：用于在类中封装函数的结构型模式
OO）类接口来封装低级函数和数据结构。常见的包装外观模式的例子是像 MFC、 ACE和 AWT这样的类库，它们封装本地的 OS C API，比如 socket、 pthreads或 GUI函数。 `_'I 9 ,.a 直接对本地 OS C ...
- 数据可视化
也不希望excel形成垄断，到时候坐地起价。程序员的求知天性，也要求随心所欲得设计，总希望是设置参数而不是选项，是调用函数而不是下拉列表框。自由是要付出代价的，就像到地摊上淘A货，或者自已起炉灶。 总会用...
- python 调用另一个py_为什么我的 Python 程序卡住啦！
所谓“不工作”就是指程序不再响应新的请求，处在了某种自娱自乐的状态，英语有一个很形象但的单词“hung”，但我不知道怎么翻译，姑且称之为“卡住”吧。本人遇到过的有两种情况，一种是卡在系...
- 20172304 2017-2018-2 《程序设计与数据结构》第六周学习总结
20172304 2017-2018-2 《程序设计与数据结构》第六周学习总结 教材学习内容总结 本周学习了数组。 首先是数组元素，数组具有优越性因为它可以声明一个能容纳多个可访问值的变量。数组的数据具有索引而且是从零开始...
- 中文版《Effective C++：改善程序与设计的55个具体做法
第一章 从C转向C++ 对每个人来说，习惯C++需要一些时间，对于已经熟悉C的程序员来说，这个过程尤其令人苦恼。因为C是C++的子集，所有的C的技术都可以继续使用，但很多用起来又不太合适。例如，C++程序员会认为...
- 110道python面试笔试题汇总
p转自：ython爬虫人工智能大数据 数据分析系列教程插播一篇...该文110道面试题全部来自于大家笔试面试时候拍照后发到群里求助的题目，并自己一道一道亲自做了，大部分题目属于巩固基本python知识点的题目，希望...
- @excel注解作用_使用 EasyPOI 优雅导出Excel模板数据（含图片）
EasyPOI功能如同名字Easy，主打的功能就是容易，让一个没接触过POI的人员可以方便的写出Excel导出，Excel模板导出，Excel导入，Word模板导出。通过简单的注解和模板语言(熟悉的表达式语法)，完成以前复杂的写法。...
- 阿里面试可以用python吗_求助！去阿里面试Python，结果被这些高级语法考“焦”了...
Python是一种极易上手的语言，但在工程应用中，由于对效率和维护性的要求，会应用到一些Python高级语法元素，这些导致了进大公司面试做题的“翻车”。考你下，以下这些语法，棒棒的你都知道吗？迭代器(iterator)。...
- 使用 EasyPOI 优雅导出Excel模板数据（含图片）
点击上方蓝色“程序猿DD”，选择“设为星标”回复“资源”获取独家整理的学习资料！作者 |星悬月来源 |blog.csdn.net/u012441819/article/details...
- python魔方程序算法_《Python基础教程》第9章 魔方方法、属性和迭代器
本章会(1)讨论一些重要的魔法方法，(2)然后会讨论两个相关的主题：属性和迭代器，(3)最后介绍一个相关的示例来处理一些稍微复杂的问题。魔法方法或者又称为特殊方法(个人挺不习惯魔法方法的翻译，叫法搞的跟拍科幻片...