#include <iostream>
using namespace std;

#include <math.h>
#define TORAD 3.1415926/180.0

float ruth(float,float,float,float);

int main()
{
  float z1,z2,ecm;
  float tmin,tmax,dt;
  
  cout << "Qual o Z1? ";
  cin >> z1;
  cout << "Qual o Z2? ";
  cin >> z2;
  cout << "Qual o Ecm? ";
  cin >> ecm;
  cout << "Qual o ThetaMin ThetaMax e DeltaTheta? ";
  cin >> tmin >> tmax >> dt;
  
  cout <<endl;
  cout <<"Usando for para o loop"<<endl;
  
  for(float t = tmin; t<=tmax; t = t+dt)
  {
    float R = ruth(z1,z2,ecm,t);
    cout <<"Theta = "<<t<<"  Sig_Ruth = "<<R<<endl;
  }
  
  cout <<endl;
  cout <<"Usando while para o loop"<<endl;
  
  float t = tmin;
  while(t<=tmax)
  {
    float R = ruth(z1,z2,ecm,t);
    cout <<"Theta = "<<t<<"  Sig_Ruth = "<<R<<endl;
    t = t + dt;
  }
  
  return 0;
}


float ruth(float Z1, float Z2, float ECM, float THCM)
{
  float ZSQ=(Z1*Z2)*(Z1*Z2);
  float CRUT=1.296*ZSQ;
  float A = 1.0/sin(THCM*TORAD/2.);
  float CSC4=A*A*A*A;
  float SIGRUT=CRUT*CSC4/(ECM*ECM);
  return SIGRUT;
}