Не работает функция синус и неправильно выводятся числа

#include <iostream>
#include <stdlib.h>
#include <time.h>
#include <iomanip>
#include <math.h>
 
#define MAX_RAND 0x7fff
 
using namespace std;
 
int main ()
{
    int u1, u3;
    double u2;
        srand(time(0));
        int S = 10;
        double x,y;
        
    
        cout << "S\t" << "U1\t" << "U2\t" << "U3\t" << "X\t" << "Y\t" << "Z\t" << endl;
    for (int i = 1; i < S; i++)
        {
            u1 = 0+rand()%3;
            u2 = -2.5 + (3.5 / (MAX_RAND +  1)) *rand();
            u3 = (6-4)/i;
            x = u1 + u2 + u3;
            y = sqrt(u2) + 0.5*u3;
            //z = 2*u1 + u2 + sin(u3);
    cout << setprecision(1) << i << "\t" << u1 << "\t" << u2 << "\t" << u3 << "\t" << x << "\t" << y << "\t"<< endl;
 
 
        }
    system("pause");
}

#include <iostream>
//#include <cstdlib>
#include <ctime>
#include <iomanip>
#include <cmath>
 
//#define MAX_RAND 0x7fff
 
using namespace std;
 
int main ()
{
    srand((unsigned)time(NULL));
    cout << "S\t" << "U1\t" << "U2\t" << "U3\t" << "X\t" << "Y\t" << "Z\t" << endl;
    const int S = 10;
    for (int i = 1; i < S; i++)
    {
        double u1 = (double)(rand()%3);
        double u2 = -2.5 + (3.5 / (RAND_MAX +  1)) *rand();
        double u3 = (6.0-4.0)/i;
        double x = u1 + u2 + u3;
        double y = sqrt(fabs(u2)) + 0.5*u3;
        double z = 2.0*u1 + u2 + sin(u3);
        cout << setprecision(4) << i << "\t" << u1 << "\t" << u2 << "\t" << u3 << "\t" << x << "\t" << y << "\t";
        cout<<z<<endl;
    }
    system("pause");
    return 0;
}