Разложение ln(1-x) в ряд Тейлора

#include <stdio.h>
#include <math.h>
 
// ln(1-x)
 
float ln1(float x, float e)
{
float r;
int n;
int sn;
int f;
float xn;
 
r = 0.0;
xn = x;
 
for(n = 0, sn = -1, f = 1; fabs((xn * (float) sn ) / (float) f) > e; n++, xn = xn * x , f++ )
r = r + (xn * (float) sn ) / (float) f;
 
return r;
}
 
void main()
{
printf("%f\n", ln1(0.2,0.001));
printf("%f\n", log(0.8));
}

#include <iostream>
#include <cmath>
 
using namespace std;
 
//факториал n
int fact(int n)
{
    return ( (n < 2) ? 1 : fact(n-1)*n);
}
 
// производная n-го порядка функции ln(x-1)
double diff (int n, double x)
{
    return ((n == 0) ? log(x-1) : pow(-1, (n+1)) * fact(n-1) / pow(x-1,n));
}
 
// разложение функции в ряд Тейлора в окрестности точки a=0 (т.е. в ряд Маклорена)
void taylor (int n)
{
    cout << "i*pi" << " ";
    for(int i = 1; i<=n; i++)
        cout << diff(i,0)/fact(i) << "*x^" << i << " ";
}
 
int main()
{
    taylor(7);
}