toyexample1

#include<bits/stdc++.h>
using namespace std;

int main()
{
	ofstream fout("toyexample1.dat");

	double energy_initial = 1.8;
	double energy_final = 3.2;
	double energy_diff = 0.000001;
	int steps = abs((energy_final - energy_initial)/energy_diff);

	double energy = energy_initial;

	double g = 3;
	double m_r = 2.5;
	double q_square = 0;
	double eta = 0.5;
	double M = 0;

	for(int i = 0;i<=steps; ++i)
	{
		double A = g*g*m_r*m_r;
		double C = A*A;
		double D = 1/(6*3.1416*6*3.1416);
		double E = C*D*q_square;
		double F = 1/(energy*energy);
		double G = E*F;
		double B = 2.667; //4/3eta;
		q_square = energy*energy/4 - m_r*m_r;
		double denom = sqrt(m_r*m_r - energy*energy + G);
		double 1bydenom = 1/denom;
		double M = B*A*1bydenom;

		fout<<energy<<'\t'<<M<<endl;

		energy+=energy_diff;
	}

	return 0;
}