м. Симпсона, не выводит значение интеграла на экран( работы на пару сек.

[maksimka95]

здравствуйте, столкнулся с проблемой...
вот код на метод Симпсона для вычисления интеграла, но я не могу понять почему не выводит значение интеграла.. а выводит просто "Otvet I=" и все..
_____________
задание:
подынтегральная ф-я: 1/( sin^2(x/3) )
шаг: П/16
точность: 0.001
отрезок: [ П/4;3П/4 ]
число разбиений = 8
______________

#include <iostream> 
#include <math.h> 
using namespace std;
double F(double x) 
{ return (1/(pow(sin(x/3),2)));
} 
int main() 
{double I0,I1,S1,S2,h,eps,a,b,x; 
int n,j; 
cout<<"Vvefite a i b"<<endl; 
cin>>a>>b; 
cout<<"Vvedite cislo razbieniy [a,b]"<<endl; 
cin>>n; 
cout<<"Vvedite tochnost"<<endl; 
cin>>eps; 
I0=1E+10; 
I1=0; 
while(fabs(I1-I0)>eps) 
{ I0=I1; 
h=(b-a)/n; 
S1=0; 
for(j=1;j<=n-1;j+2) 
{x=a+j*h; 
S1=S1+F(x); 
} 
S2=0; 
for(j=2;j<n;j+2) 
{x=a+j*h; 
S2=S2+F(x); 
} 
I1=(h/3)*(F(a)+F(b)+4*S1+2*S2); 
n=2*n; 
} 
cout<<"Otvet I="<<I1<<endl; 
system(“pause”);

return 0;
}

____________________

вот не могу понять, бесконечный цикл?
укажите пожалуйста на ошибку. Спасибо!

[alexander13]

Как минимум, в циклах for надо заменить

for(j=1;j<=n-1;j+2)

на

for(j=1;j<=n-1;j+=2)

Всю прогу не смотрел.

[maksimka95]

блииин, я так ступил жестко(((( это ж я знал( что-то даже и не обратил внимания на это.. привык уже "j++"..
спасибо огромное) программа рабочая..))