Вывод в файл

#include "stdafx.h" 
#include <iostream> 
#include <cmath> 
#include <string> 
#include <fstream>
 
double f(double x, double y) 
{ 
return (pow(x,3)); 
 
} 
 
void main()
{
    double K1, K2, K3, h, a,b,N, x, y,i,Delta,k;
    N=10;
    printf("vvedite x0=");
    scanf_s ("%lf", &x);
    printf("vvedite y0=");
    scanf_s ("%lf", &y);
    printf("vvedite a=");
    scanf_s ("%lf", &a);
    printf("vvedite b=");
    scanf_s ("%lf", &b);
    h=abs(b-a)/N;
    for (x; x<b-h; x+=h)
    {
        K1 = f(x, y);
        K2 = f(x + h / 3, y + h * K1 / 3 );
        K3 = f(x + h / 2, y + K1*h / 4 + h* K2 / 4 );
 
        y = y + h * (K1/2 - 3 * K2/2 +  2*K3);
        Delta= abs(y-f(x,y)); 
        printf("\n %lf", x);
        printf("\n %lf", y);
        printf("\n %lf", Delta); 
    }
    scanf_s ("%lf", &y);
}

void main()
{
    double K1, K2, K3, h, a,b,N, x, y,i,Delta,k;
    N=10;
    printf("vvedite x0=");
    scanf_s ("%lf", &x);
    printf("vvedite y0=");
    scanf_s ("%lf", &y);
    printf("vvedite a=");
    scanf_s ("%lf", &a);
    printf("vvedite b=");
    scanf_s ("%lf", &b);
    h=abs(b-a)/N;
 
    FILE *f = fopen("data.txt", "a");
 
    for (x; x<b-h; x+=h)
    {
        K1 = f(x, y);
        K2 = f(x + h / 3, y + h * K1 / 3 );
        K3 = f(x + h / 2, y + K1*h / 4 + h* K2 / 4 );
 
        y = y + h * (K1/2 - 3 * K2/2 +  2*K3);
        Delta= abs(y-f(x,y)); 
        printf("\n %lf", x);
        printf("\n %lf", y);
        printf("\n %lf", Delta); 
 
        fprintf(f, "%lf %lf %lf\n", x, y, Delta);
 
    }
 
    fclose(f);
 
    scanf_s ("%lf", &y);
}

fprintf(f, "%f %f %f\n", x, y, Delta);