Табулирование кусочно непрерывной функции.

#include<stdio.h>
#include<stdlib.h>
#include<conio.h>
#include<math.h>
#include<dos.h>
#define C 60
#define M 2
double f(double x)
{ double y;
    if(x<=0.5) y=sqrt(1-x);
    else if(x<2) y=exp(pow(-x,2)) ;
    else y=x*cos(x);
    return y;
 
}
void main(void) 
 
 {unsigned N, c, r;
    double x1, x2, dx, x, y, s;
    clrscr();
    printf("Program tabulates function from x1 to x2, step dx\n");
    printf("Enter x1, x2, dx\n");
    scanf("%lf%lf%lf", &x1, &x2, &dx);
    printf("Program waits 5 seconds and clears screen every N points. Enter N");
    scanf("%u", &N);
    if(dx<=0)
 
    { printf("Error. dx must be>0 dx*=-1;");
        exit(1);
    }
    clrscr();
    x=x1;
    c=C;r=1;s=0;
    do{y=f(x);
         s+=y*dx;
         gotoxy(c,r);
         printf("f(%5.2lf)=%5.2lf", x, y);
         if(r>N)
            { delay(5000);
                clrscr();
                c=C;
                r=1;
 
            }
         else{ r++; c-=M;}
         x+=dx;
 
        }
        while(x<=x2);
        getch();
 
}

do{y=f(x);
         s+=y*dx;
         gotoxy(c,r);
         printf("f(%5.2lf)=%5.2lf", x, y);
         if(r>N)
            { delay(5000);
                clrscr();
                c=C;
                r=1;
 
            }
         else{ r++; c-=M;}
         x+=dx;
 
        }
        while(x<=x2);
        getch();