1.feladat

#include<stdio.h>
#include<math.h>

void kiir (double tomb[], int db){
int i;
for (i=0; i<db; i++)
printf("%lf\n",  tomb[i]);
}
double tavolsag(double PR[], double PRV[],int n){
int i;
double osszeg=0.0;
for (i=0;i<n;++i){
    osszeg+=(PRV[i]-PR[i])*(PRV[i]-PR[i]);
    return sqrt(osszeg);
}

}

int main(void){
double L[4][4]={{0.0,0.0,1.0/3.0,0.0},
{1.0,1.0/2.0,1.0/3.0,1.0},
{0.0,1.0/2.0,0.0,0.0},
{0.0,0.0,1.0/3.0,0.0}
};

double PR[4]={0.0,0.0,0.0,0.0};
double PRV[4]={1.0/4.0,1.0/4.0,1.0/4.0,1.0/4.0};

int i,j;
for(;;){
for(i=0;i<4;i++)
PR[i] = PRV[i];
for (i=0;i<4;i++)
{
double tar=0;
for (j=0;j<4;j++)
tar+=L[i][j]*PR[j];
PRV[i]=tar;
}
    if(tavolsag(PR,PRV,4)<0.000000001)
break;
}
kiir(PR,4);
return 0;
}