Расчет выбранной функции, вид которой определяется в функции вывода Out_Rez

#include <vcl.h>
#pragma hdrstop
#include <math.h>
#include <stdio.h>
#include <iostream.h>
double fun1 (double);
double fun2 (double);
void Out_Rez (int, double, double, double);
//---------------------------------------------------------------------------
 
#pragma argsused
void main()
{  double a, b, h;
   cout<<"Input a, b, h"<<endl;
   cin>>a>>b>>h;
   cout<<"\n\t Function Y(x)"<<endl;
   Out_Rez (1, a, b, h);
   cout<<"\n\t Function S(x)"<<endl;
   Out_Rez (2, a, b, h);
   cout<<"\n Press any key..."<<endl;
      }
       double fun1 (double f, double s, double k, double x, double n) {
         s=0;
         for(k=1;k<=n;k++) {
         f=pow(x,k)*cos(k*M_PI/3)/k;
         s+=f; }
         return s;
         }
       double fun2 (double x) {
         return -0.5*log(1-2*x*cos(M_PI/3)+pow(x,2));
         }
void Out_Rez (int kod, double xn, double xk,double h) {
   double (*fun)(double);
   switch (kod) {
     case 1:
     fun = fun1;
     break;
     case 2:
     fun = fun2;
     break;  }
   for(double x=xn; x<=xk; x+=h)
     printf("x=%5.2lf \t y=%8.4lf%\n", x, fun(x));
 }