I have taken c++ but it has been about a year ago. I am taking mechanics of materials now and need to solve a bending moment problem. I am trying to get values of Sigma when a is between 0 and 20mm in 2mm increments. I have wrote a little bit of code and need to figure out if this would the correct approach? Also why is this loop infinite?
I do not have a copy of my text anymore so any help would be appreciated.

The equation I am solving for is Sigma = MC / I
Remember I need 16 values for Sigma when a is varying.

Code:
#include<iostream>
#include<cmath>
using namespace std;
 
int main()
{
 
    double n,M,C,I,Sigma,b,h,a,x,i;
 
    n = 2.66;
    C = .020;
    M = 1500;
 
    for(a=0;a<=.020;a+.002)
    {
 
           Sigma = (1500 * .020)/ ((1/12)*(.060)*pow(.040-2*a,3));
           cout<<Sigma<<endl;
     };
 
 
 
 
    system("pause");
    return 0;
 
}