我正在编写一个C ++程序来寻找大学作业的一阶微分方程的解。程序启动然后一旦我输入迭代次数,我就会收到错误消息" Euler的method.exe已停止工作"。这是我的代码:
#include <functional>
#include <vector>
using namespace std;
double f_r(double x, double r) {
return r;
}
double f_s(double x, double s) {
return -x/s;
}
double eulerstep(const function<double(double,double)>& f, double xsub0, double ysub0, double h) {
double ysub1 = ysub0+ h * f(xsub0,ysub0);
return ysub1;
}
double euler(const function<double(double,double)>& f, double xsub0, double ysub0, double h, int n) {
vector<double> xsub;
vector<double> ysub;
xsub[0] = xsub0;
ysub[0] = ysub0;
n = ysub.size();
for (int i=1; i<n; i++){
xsub[i+1] = xsub[i] + h;
ysub[i+1] = ysub[i] + h * f(xsub[i],ysub[i]);
cout << xsub[i] << " , " << ysub[i] << endl;
}
return ysub[n];
}
int main() {
int nsteps = 0;
cout << "Number of steps?" << endl;
cin >> nsteps;
double h = 1.0 / nsteps;
double r = euler(f_r,0,1,h,nsteps);
double s = euler(f_s,0,1,h,nsteps);
return 0;
}
&#13;
我很确定问题在于我如何定义我的向量,但我不习惯使用它们,因此无法查看我出错的地方。如果有人能指出我方法中的错误,我将非常感激
由于
答案 0 :(得分:1)
这是一个更新的功能。您必须设置大小,而不是询问空向量的大小。
double euler(const function<double(double, double)>& f, double xsub0, double ysub0, double h, int n)
{
vector<double> xsub;
vector<double> ysub;
xsub.resize(n+1); // so we can access [n], it must be size n+1
ysub.resize(n+1);
xsub[0] = xsub0;
ysub[0] = ysub0;
for (int i = 1; i<n; i++) {
xsub[i + 1] = xsub[i] + h;
ysub[i + 1] = ysub[i] + h * f(xsub[i], ysub[i]);
cout << xsub[i] << " , " << ysub[i] << endl;
}
return ysub[n];
}
答案 1 :(得分:0)
${rvm-prompt}
您正在实例化一对向量。它们最初是空的。
vector<double> xsub;
vector<double> ysub;
然后继续为向量的内容赋值。这是崩溃的地方,因为向量是空的,并且不包含任何值。
xsub[0] = xsub0;
ysub[0] = ysub0;
引用向量中的现有元素。向量必须已经包含至少vector[x]
个元素,但是没有任何向量中包含任何内容。它们没有元素0,也没有元素1,也没有任何元素。
您的代码看起来希望每个向量都包含x+1
个元素,因此在尝试使用每个向量之前,应该显式调用每个向量的n+1
方法。