fork download
copy 
  1. #include <iostream>
  2. #include <ctime>
  3. #include <math.h>
  4. #include <cmath>
  5. #include <cstdlib>
  6. using namespace std;
  7.  
  8. double f(double x)
  9. {
  10. 	return x*x;
  11. }       
  12. int pow(int a, int b) {
  13.     int result = 1;
  14.     for (int i = 0; i < b; i++)
  15.         result *= a; 
  16.     return result;
  17.     }
  18. int main()
  19. {
  20. 	double m;
  21. 	double h;
  22. 	double t[65000];
  23. 	double T_alpha;
  24. 	double delta;
  25. 	double delta2;
  26. 	double N;
  27. 	double x[200];
  28. 	double f_del2[200];
  29. 	double A[200];
  30. 	double alpha;
  31. 	double n;
  32. 	double h1;
  33. 	cout << "m=";
  34. 	cin >> m;
  35. 	cout << endl;
  36. 	cout << "delta=";
  37. 	cin >> delta;
  38. 	cout << endl;
  39. 	cout << "delta2=";
  40. 	cin >> delta2;
  41. 	cout << endl;
  42. 	h = 1 / m;
  43. 	alpha = delta2;
  44. 	n=m;
  45. 	h1=(2*alpha)/m;
  46. 	cout<<h1<<endl;
  47. 	cout << "h=" << h << endl;
  48. 	cout << "alpha=" << alpha << endl;
  49. 	for (int i = 0; i < m; i++)
  50. 	{
  51. 		x[i] = h*i;
  52. 		cout << "x[" << i << "]=" << x[i] << endl;
  53. 	}
  54. 	for(int i = 0; i < n - 1; i++) {
  55.             if (( ((x[i] - alpha) >= 0) && ((x[i] - alpha) <= 1) ) && ( ((x[i] + alpha) >= 0) && ((x[i] + alpha) <= 1) ))
  56. 			{
  57.                   //double j = x[i] - alpha;
  58.                   //int tt = 0;
  59.                   int g=0;
  60. 				  double a=x[i]-alpha;
  61. 				  cout<<"a"<<a<<endl;
  62. 				  double b=x[i]+alpha;
  63. 				  cout<<"b"<<b<<endl;
  64. 				  while (a<=b)
  65. 				  {
  66. 				  t[g]=a;
  67. 				  a+=h1;
  68. 				  g++;
  69. 				  }
  70.                  /* while(j <= x[i] + alpha) {
  71.                           t[tt] = j;
  72.                           tt++;
  73.                           }*/
  74.             } else {
  75.                    cout << "Ne vhodit" << endl;
  76.               }
  77.             }
  78. 	N = sqrt((0, 5 * (m*pow(delta / delta2, 2) - 1)) + 1);
  79. 	cout << "N=" << N << endl;
  80. 	for (int i = 0; i < m; i++)
  81. 	{
  82. 		A[i] = 1 + rand() % 10;
  83. 		cout << "A[" << i << "]=" << A[i] << endl;
  84. 	}
  85. 	for (int i = 0; i < m; i++)
  86. 		f_del2[i] = f(x[i]) +pow(1, i) * A[i] * delta2;
  87. 	for (int i = 4; i < m; i++)
  88. 		f_del2[i] = f(x[i]) + pow(-1, i)*A[i] * delta2;
  89. 	f_del2[2] = f(x[2]) + N*A[2] * delta2;
  90. 	f_del2[3] = f(x[3]) - N*A[3] * delta2;	
  91. 	for (int i = 0; i < m; i++)
  92. 		cout << "f(x[" << i << "])=" << f(x[i]) << endl;
  93. 	for (int i = 0; i < m; i++)
  94. 		cout << "f_del2[" << i << "]=" << f_del2[i] << endl;
  95. 		for(int i = 0;i < n; i++) {
  96.             T_alpha+= 3 / (2 * pow(alpha,3)) * ( (t[i] - x[i]) * f_del2[i]);
  97.             cout << "i = " << i << " T_alpha = " << T_alpha << endl;
  98.             }
  99. 	system("PAUSE");
  100. 	return 0;
  101. }
Success #stdin #stdout #stderr 0s 16456KB
comments (?)
stdin
copy 
10
0.004
0.008
stdout
copy 
m=
delta=
delta2=
0.0016
h=0.1
alpha=0.008
x[0]=0
x[1]=0.1
x[2]=0.2
x[3]=0.3
x[4]=0.4
x[5]=0.5
x[6]=0.6
x[7]=0.7
x[8]=0.8
x[9]=0.9
Ne vhodit
a0.092
b0.108
a0.192
b0.208
a0.292
b0.308
a0.392
b0.408
a0.492
b0.508
a0.592
b0.608
a0.692
b0.708
a0.792
b0.808
N=2.91548
A[0]=4
A[1]=7
A[2]=8
A[3]=6
A[4]=4
A[5]=6
A[6]=7
A[7]=3
A[8]=10
A[9]=2
f(x[0])=0
f(x[1])=0.01
f(x[2])=0.04
f(x[3])=0.09
f(x[4])=0.16
f(x[5])=0.25
f(x[6])=0.36
f(x[7])=0.49
f(x[8])=0.64
f(x[9])=0.81
f_del2[0]=0.032
f_del2[1]=0.066
f_del2[2]=0.22659
f_del2[3]=-0.0499428
f_del2[4]=0.192
f_del2[5]=0.202
f_del2[6]=0.416
f_del2[7]=0.466
f_del2[8]=0.72
f_del2[9]=0.794
i = 0 T_alpha = 74250
i = 1 T_alpha = 208364
i = 2 T_alpha = 603481
i = 3 T_alpha = 530791
i = 4 T_alpha = 754891
i = 5 T_alpha = 932430
i = 6 T_alpha = 1.17813e+06
i = 7 T_alpha = 1.31902e+06
i = 8 T_alpha = 1.32915e+06
i = 9 T_alpha = 1.11142e+06
stderr
copy 
sh: 1: PAUSE: not found
https://ideone.com/LuoaH5
language:
C++ (gcc 8.3)
created:
8 years ago
visibility:
public

Share or Embed source code

Discover > Sphere Engine API

The brand new service which powers Ideone!

Discover > IDE Widget

Widget for compiling and running the source code in a web browser!