Проблемма с функцией,подскажите что не так

#include "stdafx.h"
#include<iostream>
#include <clocale>
#include <cmath>
using namespace std;
 
 
 
 
void input (double x, double e);
//функция для ввода значений
 
double solut (double &result,double &tochnost,double &e,double &x);
//решение 
 
void output(double result,double tochnost);
//функция для вывода результата
 
 
void main()
{
    setlocale(LC_ALL, "Russian");
double x=0,e=0,tochnost,result;
input(x,e);
solut(tochnost,result,x,e);
result=solut(x,e,tochnost,result);
tochnost=solut(x,e,tochnost,result);
output (result,tochnost);
 
 
 
}
 
 
void input (double x, double e)
{double alf=0;
cout<<"Введите значение альфа: ";
    cin>>alf;
    cout << "x = ";
    cin >> x;   
    cout << "Точное e = ";
    cin >> e; // 0.01, 0.001, 0.0001 
}
 
double solut (double &result,double &tochnost,double &e,double &x)
{
    double n=0, fact=0, j,alf=0; // n - степень
double a = 1, slagaemoe,k=n;
slagaemoe = cos(alf);
    while(fabs(slagaemoe)>e)
    {                       
    result = result + slagaemoe;    
    n++; 
    k=n;
    j=1; // счётчик факториала
    fact=1;
    while(j<=k) 
        {       
            fact=fact*j;
            j++;            
    //}     
        a = a*x;        
 
        slagaemoe = a*cos(k*3.14/2+alf)/fact;
    }
        return (tochnost);
        
 
}
}
void output (double result,double tochnost)
{
cout << endl << result << endl << tochnost << endl;
 
 
}

void input (double& x, double& e);

#include <iostream>
#include <clocale>
#include <cmath>
 
using namespace std;
 
void input(double& alpha, double& x, double& e);
double solut(double alpha, double x, double e);
 
int main()
{
   double alpha, x, e;
 
   input(alpha, x, e);
 
   cout << "result: " << solut(alpha, x, e) << endl;
 
   return 0;
}
 
void input(double& alpha, double& x, double& e)
{
   cout << "Введите значение альфа: ";
   cin >> alpha;
   cout << "x = ";
   cin >> x;
   cout << "Точное e = ";
   cin >> e;
}
 
double solut(double alpha, double x, double e)
{
   double result = 0.0;
 
   unsigned fact = 1;
   double localX = 1;
   unsigned n = 1;
 
   double summand = cos(alpha);
 
   while (fabs(summand) > e)
   {
      result += summand;
      summand = localX * cos(alpha + n * 3.14 / 2) / fact;
 
      fact *= ++n;
      localX *= x;
   }
 
   return result;
}