qq_42756072的博客

#include<bits/stdc++.h>using namespace std;typedef long long ll;const ll N = 1e3;double f1(double x){ double ans = 1.0 / (1.0+ sin(x) * sin(x) ); return ans;}double f2(double x){ re...

#include<bits/stdc++.h>#define e 2.718281828489#define eps 0.5*1e-5#define eps2 1e-12#define inf 0x3f3f3f3fusing namespace std;typedef long long ll;const ll N = 1e5;double f1(double x)...

比较简单，照着公式算就可以了。#include<bits/stdc++.h>using namespace std;#define inf 0x3f3f3f3f#define eps 1e-6typedef long long ll;const ll N = 1e3;double x[N];double y[N];double f(double x, dou...

最近一段时间再忙期末考试，小学期课程设计的东西，没怎么更新博客....更新一个用四阶龙格库塔方法求解脉冲微分方程，题目来源是一篇论文《Impulsive control of projective synchronization in chaotic systems 》。概述：做法其实也很简单，就是原来求一个因变量的龙格库塔公式变成求五个因变量的龙格库塔公式，当然，变化求...