1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
| #include<iostream>
#include<math.h>
using namespace std;
void metodSympsona(double a,double b)
{
double x[41],y[41],iks[41],igrik[41],m,h,integral,s,sum=0;int i,j,k;
cout<<"vvedite kolvo x m";
cin>>k;
switch(k)
{
case 1:
m=10;
break;
case 2:
m=20;
break;
case 3:
m=40;
break;
default:
cout<<"Vibrano nedopustimoe zna4enie";
}
h=(b-a)/m;
for(i=1;i<m;i++)
{
x[i]=a+(i-1)*h;
y[i]=0.1*pow(x[i],3)+pow(x[i],2)-10*sin(x[i]);
iks[i]=(x[i]+x[i+1])/2;
igrik[i]=0.1*pow(iks[i],3)+pow(iks[i],2)-10*sin(iks[i]);
}
for(i=1;i<m;i++)
{
s=y[i]+4*igrik[i]+y[i+1];
sum=sum+s;
}
integral=h/6*s;
cout<<integral;
}
int main (void)
{
double a=-4,b=2;
metodSympsona(a,b);
system("pause");
return(0);
} |