Ask Question

Name:
Title:
Your Question:

Answer Question

Name:
Your Answer:
User Submitted Source Code!


Description:
  jjjj
Language: C/C++
Code:
#include <iostream>
#include <conio.h>
#include <math.h>
const int n = 3;

const int N = 100;

 
using namespace std;

 
double
f (double x) 
{
  
return (2 * pow (x, 3)) / (pow (x, 4));

}


 
double
Gauss (double a, double b) 
{
  
const double Xi[n] = { -0.7745967, 0, 0.7745967 };
  
const double Ci[n] = { 0.5555556, 0.8888889, 0.5555556 };
  
 
double ra = (b - a) / 2;
  
double su = (a + b) / 2;
  
double Q, S = 0.0;
  
for (int i = 0; i < n; i++)
    
    {
      
Q = su + ra * Xi[i];
      
S += Ci[i] * f (Q);
    

return ra * S;

}


 
int
main () 
{
  double a, b;
  
cout << "a = ";
  
cin >> a;
  
cout << "b = ";
  
cin >> b;
  
double s = 0.0;
  
for (int i = 0; i < N; ++i)
    
    {
      
s += Gauss (a + i * (b - a) / N, a + (i + 1) * (b - a) / N);
    

cout << "I = " << s << endl;
  
 
return 0;

}
          
          
Comments: