/***************************************************************************
 *   Copyright (C) 2007 by Félix Robles Elvira   *
 *   f@ubuntu   *
 *                                                                         *
 *   This program is free software; you can redistribute it and/or modify  *
 *   it under the terms of the GNU General Public License as published by  *
 *   the Free Software Foundation; either version 2 of the License, or     *
 *   (at your option) any later version.                                   *
 *                                                                         *
 *   This program is distributed in the hope that it will be useful,       *
 *   but WITHOUT ANY WARRANTY; without even the implied warranty of        *
 *   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the         *
 *   GNU General Public License for more details.                          *
 *                                                                         *
 *   You should have received a copy of the GNU General Public License     *
 *   along with this program; if not, write to the                         *
 *   Free Software Foundation, Inc.,                                       *
 *   59 Temple Place - Suite 330, Boston, MA  02111-1307, USA.             *
 ***************************************************************************/


#ifdef HAVE_CONFIG_H
#include <config.h>
#endif

#include <iostream>
#include <cstdlib>
#include <iostream>
#include <math.h>
#define e0 8.8541878176e-12
#define pi 3.1415926535897932384626433832795028841971693993751058209749445923078164062862089986280348253421170679f


using namespace std;

class miclase{
      public:
      double sumq[2];
      
};

inline float p (int i){
	return (float)(i%2);
}
double bi[1000],qi[1000];
miclase suma(double r1, double r2, double v1, double v2, double d, bool k);
miclase sacarq(double r1, double r2, double v1, double v2, double d);
int main (int argc, char *argv[])
{
	while(1){
		double r1, r2, v1,  v2,  d;
		cout << "Pi es 3 exactamente: double r1, r2, v1,  v2,  d" << endl;
		cin >> r1 >> r2 >> v1 >> v2 >> d;
		miclase s = sacarq(r1, r2, v1,  v2,  d);
		double c1 = 4*pi*e0*r1;
		double c2  = 4*pi*e0*r2;
		double c12 = (s.sumq[0]-c1*v1)/v2;
		double c21 = (s.sumq[1]-c2*v2)/v1;
		cout << "c12 c21 c1 c2: " << c12 << " " << c21 << " " << c1 << " " << c2 << endl;
		cout << "q1 q2: " << s.sumq[0] << " " << s.sumq[1]  << endl;
	}
  
  system("pause");
  return 0;
}
//Saca las q
miclase sacarq(double r1, double r2, double v1, double v2, double d){
	miclase f,s;
        f = suma(r1, r2, v1,  v2,  d,0);
        s = suma(r1, r2, v1,  v2,  d,1);
	cout << "f q1 q2: " << f.sumq[0] << " " << f.sumq[1]  << endl;
	cout << "s q1 q2: " << s.sumq[1] << " " << s.sumq[0]  << endl;
        f.sumq[0]+=s.sumq[1]; f.sumq[1]+=s.sumq[0];
        return f;
}
miclase suma(double r1, double r2, double v1, double v2, double d, bool k)
{
        miclase s;
       if(k)
       {    double foo= r1;
            r1=r2;
            r2=foo;
            foo = v1;
            v1=v2;
            v2=foo;
       }
       qi[0]= 4*pi*e0*r1*v1;
       bi[0]=0;//(r1*r1)/d;
       s.sumq[0]=s.sumq[1]=0;
	cout << endl << "b0" << ": " << bi[0] << endl;
       for(int i =1; i< 30; i++)
       {         
                 bi[i]=( (r1*r1)*((i+1)%2)+(r2*r2)*(i%2)) / (d-bi[i-1]);
                /* if(bi[i]<0)
                            bi[i]= -bi[i]; */
                 cout << endl << "b" << i << ": " << bi[i] << endl;
       }
       cout << endl << "q0" << ": " << qi[0] << endl;
       for(int i =1; i< 30; i++)
       {         
                           qi[i]=qi[i-1]*( bi[i]-r1*p(i)-r2*p(i+1) )
                /(d-bi[i-1]-r1*p(i)-r2*p(i+1));
		cout << endl << "q" << i << "= " << qi[i]  << " = ";
		cout   <<qi[i-1]<< " *( "<< bi[i] << " - "<< r1 << " * "<<  p(i) << " - "<< r2 << " * "<< p(i+1) << " ) / ( "<<  d<< " - "<< bi[i-1] << " - "<< +r1*p(i)+r2*p(i+1) << " ) "<< endl;
               /* if(qi[i]<0)
                           qi[i]= -qi[i];*/
                
       }
       for(int i =0; i< 30; i++)
       {
               //  cout << "s.sumq[0]: " <<  s.sumq[0] << " s.sumq[1]: " <<  s.sumq[1] << endl;
               if(i%2 == 0)
               {        
                      s.sumq[0]+=qi[i];
               }
               else
               {      
                      s.sumq[1]+=qi[i];
               }
       }
       return s;
}