Вычислить значение y в зависимости от выбранной функции f(x)

#include <stdio.h> //printf scanf
#include <stdlib.h>//system
#include <math.h>  //sqrt, pow
 
double X(double z)
{
    double x = z;
    double ipart;
    double fpart = modf(x,&ipart);
    if(0 <= z)
        x = sin(z);
    return x;
}
 
//Âîò òàê ïðîñòî îïèñûâàåì ôóíêöèþ y = f(x) 
//âîçâðàùàåìîå çíà÷åíèå - çíà÷åíèå ôóíêöèè â òî÷êå õ
double f(double x)
{
    return 2*x;
}
 
int main()
{
    system("chcp 1251");//Ïðîñòàÿ ðóññèôèêàöèÿ
    double x, y, z, a = 0 , b = 0;
    double zn = 0, zk = 0, hz = 0;
 
    //Ââîä à b
    printf("Ââåäèòå a  : ");scanf("%lf",&a);
    printf("Ââåäèòå b  : ");scanf("%lf",&b);
    
    //Ââîäèì øàã è äèàïàçîí çíà÷åíèé àðãóìåíòà
    printf("Ââåäèòå zn : ");scanf("%lf",&zn);
    printf("Ââåäèòå zk : ");scanf("%lf",&zk);
    printf("Ââåäèòå øàã: ");scanf("%lf",&hz);
    //Ïå÷àòàåì øàïêó òàáëèöû
    printf("------------------\n");
    printf("  z  |  x  |  y  |\n");
    printf("------------------\n");
    for(z = zn; z <= zk; z += hz)
    {
        x = X(z);
        y = 2*a*pow(sin(x),2)/3 - 3*b*pow(cos(f(x)),2)/4;
        printf("%.2f | %.2f | %.2f\n",x,y,z);
    }
    system("pause");//Ñòîï-òî÷êà âêîíöå ïðîãðàììû
    return 0;
}