wyznacznik

#include <stdio.h> //pliki nagłówkowe
#include <conio.h>
#include <math.h>


double x[31],y[31],z[31][31];
double det;

void zad1();
void zad2();
void RUN(void);


int main()
{
    RUN();
    return 0;
}

void RUN(void)
{
    zad1();

    zad2();

}

void zad1()
{
    double x1,x2,y1,y2,krokX,krokY;
    int k,i,j;

    do{
        printf("\n Podaj wartosc poczotkowa X:");
        k=scanf("%lf",&x1);
        if (k==0) printf("\n\n Niepoprawny format liczby!\n\n");
        fflush(stdin);
    }while (k==0);

    do{
        printf("\n Podaj wartosc koncowa X:");
        k=scanf("%lf",&x2);
        if (k==0||x2<x1) printf("\n\n Niepoprawny format liczby!\n\n");
        fflush(stdin);
    }while (k==0||x2<x1);

    do{
        printf("\n Podaj wartosc poczotkowa Y:");
        k=scanf("%lf",&y1);
        if (k==0) printf("\n\n Niepoprawny format liczby!\n\n");
        fflush(stdin);
    }while (k==0);

    do{
        printf("\n Podaj wartosc koncowa Y:");
        k=scanf("%lf",&y2);
        if (k==0||y2<y1) printf("\n\n Niepoprawny format liczby!\n\n");
        fflush(stdin);
    }while (k==0|y2<y1);

    krokX=(x2-x1)/31;
    krokY=(y2-y1)/31;
    x[0]=x1;
    y[0]=y1;

    for(i=1;i<31;i++)
    {
        x[i]=x[i-1]+krokX;
        y[i]=y[i-1]+krokY;
    }
    for(i=0;i<31;i++)
    {
        printf("\n[%d]:  %lf   %lf",i,x[i],y[i]);
    }

    for(i=0;i<31;i++)
        for(j=0;j<31;j++)
            z[i][j]=pow((x[i]-15),2)-pow((y[j]-15),2);

    for(i=0;i<31;i++) //dr. Wojciechowski zalecił użyć przekątnej z tym wzorem aby wyznacnzik nie wychodził 0
        z[i][i]=0.05*(i+1);


    FILE* plik=fopen("zad1.txt", "w");
    for(i=0;i<31;i++)
    {
        fprintf(plik,"\n");
        for(j=0;j<31;j++){
            fprintf(plik,"%lf ",z[i][j]);
        }
    }
    fclose(plik);
    printf("\n\n");


}

void zad2()
{
    printf("Transpozycja macierzy\n");
    int i,j;
    FILE* plik2=fopen("zad2.txt", "w");
    for(i=0;i<31;i++)
    {
        fprintf(plik2,"\n");
        for(j=0;j<31;j++){
            fprintf(plik2,"%lf ",z[j][i]);
        }
    }
    fclose(plik2);
    printf("\n\n");

}