Метод Эйлера // ошибка NAN(IND)

[Nurik1509]

#include "stdafx.h"
#include <iostream>
#include <fstream>
#include <math.h>
#define n 200
using namespace std;

int main()
{
	ofstream fout("a2.txt", ios::out);

	long double x1[n], x2[n], x3[n], u[n], dt = 0.01, t;
	long double m = 729.0, l2;
	int i;
	x1[0] = 2.0;
	x2[0] = 5.0;
	x3[0] = 3.0;
	l2 = 3.0;

	for (i = 0; i < n-1; i++)
	{
		t = i*dt;
		x1[i + 1] = x1[i] + dt*(x2[i]* x2[i]*x2[i]+(x1[i]*x1[i]*x2[i])/(1.0+ x2[i]* x2[i])*log(1.0+ x1[i]* x1[i]));
		u[i] = -pow(m, 4 / 3)*pow((10.0*(x3[i] + l2*(x1[i] - pow(sin(t), 3)))) + 16.0 * (x1[i] - pow(sin(t), 3)), 3);
		x2[i + 1] = x2[i] + dt*u[i];
		x3[i + 1] = x3[i] + dt*(u[i]/pow(m,1/3)-m*l2*pow((x3[i]+l2*(x1[i] - pow(sin(t), 3))), 3));
	}

	for (i = 0; i <= n-1; i++)
	{
		fout << i*dt << " " << x1[i] << endl;
	}
	
	return 0;
}

выдает ошибку nan (ind) // не могу понять // то ли слишком большая цифра или в коде ошибка //

[alexzk]

Вот тут
4 / 3 = 1
1/3 = 0
Нужно

4./ 3.
1./3.

ВОбщем у вас там кругом автоматом в целое преобразовыается, попробуйте верно числа записать для начала.

[Nurik1509]

#include "stdafx.h"
#include <iostream>
#include <fstream>
#include <math.h>
#define n 200
using namespace std;

int main()
{
	ofstream fout("a2.txt", ios::out);

	long double x1[n], x2[n], x3[n], u[n], dt = 0.01, t;
	long double m = 729.0, l2;
	int i;
	x1[0] = 2.0;
	x2[0] = 5.0;
	x3[0] = 3.0;
	l2 = 3.0;

	for (i = 0; i < n-1; i++)
	{
		t = i*dt;
		x1[i + 1] = x1[i] + dt*(x2[i]* x2[i]*x2[i]+(x1[i]*x1[i]*x2[i])/(1.0+ x2[i]* x2[i])*log(1.0+ x1[i]* x1[i]));
		u[i] = -pow(m, 4.0 / 3.0)*pow((10.0*(x3[i] + l2*(x1[i] - pow(sin(t), 3)))) + 16.0 * (x1[i] - pow(sin(t), 3)), 3);
		x2[i + 1] = x2[i] + dt*u[i];
		x3[i + 1] = x3[i] + dt*(u[i]/pow(m,1.0/3.0)-m*l2*pow((x3[i]+l2*(x1[i] - pow(sin(t), 3))), 3));
	}
	for (i = 0; i <= n-1; i++)
	{
		fout << i*dt << " " << x1[i] << endl;
	}
	return 0;
}

попробовал так // все равно результат так же ///
первые 5-6 раза с начальными значениями вычисляет /// а потом nan-ы пойдет


0 2
0.01 3.26238
0.02 -1.69102e+22
0.03 3.56347e+76
0.04 -3.25201e+239
0.05 -nan(ind)
0.06 -nan(ind)
0.07 -nan(ind)
0.08 -nan(ind)
0.09 -nan(ind)
0.1 -nan(ind)
0.11 -nan(ind)
0.12 -nan(ind)
0.13 -nan(ind)
0.14 -nan(ind)
0.15 -nan(ind)
0.16 -nan(ind)
0.17 -nan(ind)
0.18 -nan(ind)
0.19 -nan(ind)
0.2 -nan(ind)
0.21 -nan(ind)


а вообще это задача слежения // есть опорный график (синус куб) и моя система (график) должна повторять за ним

здесь только второй график хотел нарисовать /// не знаю может в коде ошибка

[Nurik1509]

может задачу не понял

[Аватар]

Проверь формулу для u[i], в проге оно должно немерено расти и растет ))

[p51x]

Вы хоть представляете, что за число "-3e+239"?.. А что будет если его умножить или на него поделить?

[Nurik1509]

Сообщение от p51x

Вы хоть представляете, что за число "-3e+239"?.. А что будет если его умножить или на него поделить?

вообще по идее оно не должно иметь очень большое число /// проверил уже по формуле // вот не могу найти ошибку

[Nurik1509]

Сообщение от Nurik1509

вообще по идее оно не должно иметь очень большое число /// проверил уже по формуле // вот не могу найти ошибку

формулу можете посмотреть? и если сможете показать на ошибку в коде