Чтение файла

#include <stdio.h>
#include <stdlib.h>
#include <string.h> 
#include <math.h>
#define eps 0.00000001
 
typedef struct Poly{
    
    double koef;
    
} Poly;
 
int Enter(Poly*p,int n)
{
    int i;
    if(p==NULL){
            
        printf("error(malloc)");    
    }
    if ((n>50) || (n<=0)){
    return printf ("Error size(polinom)");
    }
    printf("Vvedite kofficienti polinoma\n");
    for (i=n;i>=0;i--)
    { 
    printf("k[%d]=",i);
    scanf ("%lf", &p[i].koef);
    }
    printf("Polinom:");
    for (i=n;i>=0;i--)
    { 
        if(i>1)
            printf("(%g)*x^%d + ", p[i].koef, i);
        else
            printf("(%g)%s",p[i].koef, i==1?"*x + ":"\n");
    }   
} 
 
void Output(Poly*p,int n)
{
    int i;
    printf("Polinom:");
    for (i=n;i>=0;i--)
    { 
        if(i>1)
            printf("(%g)*x^%d + ", p[i].koef, i);
        else
            printf("(%g)%s",p[i].koef, i==1?"*x + ":"\n");
    }   
}
 
double Res(Poly*p,int n)
{
    int max_iter=100000;
    double k[50],x;
    double f,df=0;
    int i,q=0;
    printf("Vvedite nachalo intervala\n");
    scanf("%lf",&x);
    while(1)
    {  
       f=0;
       df=0; 
    for (i=n;i>=0;i--)
    { 
        df=df+p[i].koef*i* pow(x,i-1);
        f=f+p[i].koef* pow(x,i);
    }
    q++;
    if (q>max_iter){
    return printf("\nNet kornei\n");
    }
    if (fabs(f)<eps)
        {   
            printf("\nKoreni' x=%g\n",x);
            printf("!!!Esli hotite prervati cicle vvedite ZERO(0)!!!\n");
            printf("Vvedite sledyishee priblizhenie: New x=");
            scanf("%lf",&x);
        }
        else
            x=x-(f/df);
        }
    return 0;
}
 
 
double fResMethod(Poly*p,int n,double x)
{ 
     int i;
     double f=0;
     for (i=n;i>=0;i--)
    { 
      f=f+p[i].koef*pow(x,i);
    }  
          return f;
}
 
double MethodKacatelinix(Poly*p,int n,double x0,double x1) 
{
 double tmp;
 int k=0;
    while(fabs(x1-x0)>eps)   //Âûáðàí îñòàíîâ |x[n]-x[n+1]|<eps
    {
     tmp=x1;
     x1=x1-(x1-x0)*fResMethod(p,n,x1)/(fResMethod(p,n,x1)-fResMethod(p,n,x0));
     x0=tmp;
     k++;
   }
     printf("METHOD KACATELINIH: kol-vo iteraciy k=%d, koren' x=%lf, f(x)=%lf \n",k,x1,fResMethod(p,n,x1));
     return x1;
} 
    
double sign(double x)
{
       return x>0 ? 1.0: x<0 ? -1.0 : 0.0;
}
 
double fx(double x)
{
    return (pow(x,2))-8*x+4;    
}
 
double dfx1(double x)
{
    return 2*x-8;
}
 
double dfx2(double x)
{
    return ((fx(x+eps)-fx(x-eps))/(2.*eps));
}
 
 
double dfx3(double x)
{
    return (1./eps)*(2./3.*(fx(x+eps)-fx(x-eps))-(1./12.)*(fx(x+2.*eps)-fx(x-2.*eps)));     
}
 
double dihotomia(Poly*p,int n,double a, double b)
{   
    double c;
    int k=0;
    while(fabs(b-a)>=eps)
{ 
    c=(a+b)/2.0;
    if ((fResMethod(p,n,c))>=0)  b=c;
    else 
     a=c;
     k++;
}
printf("METHOD DIHOTOMII: kol-vo iteraciy k=%d, koren' x=%lf, f(x)=%lf \n",k,c,fResMethod(p,n,c));
}
 
 
double methodChord(Poly*p,int n,double a, double b)
{
int k=0;
// a, b - ïðåäåëû õîðäû
 
    while(fabs(b - a) >= eps)
    {
        a = b - (b - a) * fResMethod(p,n,b)/(fResMethod(p,n,b) - fResMethod(p,n,a));
        b = a + (a - b) * fResMethod(p,n,a)/(fResMethod(p,n,a) - fResMethod(p,n,b));
        k++;
    }
// a - i-1, b - i-òûé ÷ëåíû
printf("METHOD CHORD: kol-vo iteraciy k=%d, x=%lf, f(x)=%lf\n", k,b,fResMethod(p,n,b));
}
 
double methodSekush(Poly*p,int n,double a, double b)
{
double tmp;
int k=0;
    while(fabs(b-a)>= eps)
     { 
     tmp=b;
     b=b-(b-a)*fResMethod(p,n,b)/(fResMethod(p,n,b)-fResMethod(p,n,a));
     a=tmp;
     k++;
   }
printf("METHOD SEKUSH: kol-vo iteraciy k=%d, x=%lf, f(x)=%lf\n", k,b,fResMethod(p,n,b));
}
 
 
double methodParabol(Poly*p,int n,double a, double b)
{ int k=0; double c,q,l,A,B,C;
      c=(a+b)/2;
      A=((fResMethod(p,n,b)-fResMethod(p,n,c))/(b-c)-(fResMethod(p,n,c)-fResMethod(p,n,a))/(c-a))/(b-a);
      B=((fResMethod(p,n,c)-fResMethod(p,n,a))/(c-a)+A*(a-c));
      C=fResMethod(p,n,a);
      l=a-2*C/(B+sign(sqrt(pow(B,2.)-4.*A*C)));
      k++;
printf("METHOD PARABOL: kol-vo iteraciy k=%d, x=%lf, f(x)=%lf\n", k,l,fResMethod(p,n,l));
}
 
double methodParabolMullera(Poly*p,int n,double a, double b)
{ int k=0; double c,q,l,A,B,C;
      c=(a+b)/2;
      q=(a-c)/(c-b);
      A=q*fResMethod(p,n,a)-q*(1+q)*fResMethod(p,n,c)+q*q*fResMethod(p,n,b);
      B=(2*q+1)*fResMethod(p,n,a)-pow((1+q),2.)*fResMethod(p,n,c)+q*q*fResMethod(p,n,b);
      C=(1+q)*fResMethod(p,n,a);
      l=a-(a-c)*2*C/(B-(sqrt(B*B-4.*A*C)));
      k++;
printf("METHOD PARABOL-MULLERA: kol-vo iteraciy k=%d, x=%lf, f(x)=%lf\n", k,l,fResMethod(p,n,l));
}
 
double methodRiddera(Poly*p,int n,double a, double b)
{ 
   int k=0; double c,d;
      
   while(fabs(b-a)>=eps)
   {
       c=(b+a)/2.0;
       d=c+(c-a)*sign( fResMethod(p,n,a)-fResMethod(p,n,b) )*fResMethod(p,n,c)/
          sqrt(fabs(pow(fResMethod(p,n,c),2.) - fResMethod(p,n,a)*fResMethod(p,n,b)));
       if((fResMethod(p,n,a)*fResMethod(p,n,c))>0)
          b=c;
       else
          a=c;
       k++;
   }                    
printf("METHOD RIDDERA: kol-vo iteraciy k=%d, x=%lf, f(x)=%lf\n", k,d,fResMethod(p,n,d));       
}
 
double methodBrenta(Poly*p,int n,double a, double b)
{ int k=0;double c,d;
   {
       c=(a*fResMethod(p,n,b)-b*fResMethod(p,n,a))/(fResMethod(p,n,b)-fResMethod(p,n,a));
       d=(c*fResMethod(p,n,a)-a*fResMethod(p,n,c))/(fResMethod(p,n,a)-fResMethod(p,n,c));
       if(d>((b+c)/2))  d=(b+c)/2;
       if (d>c+eps) d=c+eps;
       k++;      
   }
printf("METHOD BRENTA: kol-vo iteraciy k=%d, x=%lf, f(x)=%lf\n", k,c,fResMethod(p,n,c));      
}
    
double methodDekkeraBrenta(Poly*p,int n,double a, double b) //http://simenergy.ru/math-analysis/solution-methods/46-method-dekker-brent
{ 
       int k=0; double c,l;
   {
       c=a+(b-a)*0.382;
       l=c+fResMethod(p,n,c)/fResMethod(p,n,a)*(fResMethod(p,n,a)/fResMethod(p,n,a)*(fResMethod(p,n,a)/fResMethod(p,n,b)-fResMethod(p,n,a)/fResMethod(p,n,b))*(b-c)-(1-fResMethod(p,n,c)/fResMethod(p,n,b))*(c-a))
        /((fResMethod(p,n,a)/fResMethod(p,n,b)-1)*(fResMethod(p,n,c)/fResMethod(p,n,b)-1)*(fResMethod(p,n,c)/fResMethod(p,n,a)-1));
       k++;      
   }
printf("METHOD DEKKERA-BRENTA: kol-vo iteraciy k=%d, x=%lf, f(x)=%lf\n", k,l,fResMethod(p,n,l));      
}   
 
 
double methodIteration(Poly*p,int n,double a, double b)
{ 
  int k=0; double c;
  while (fabs(b-a)>=eps)      
  {
    c=a-fResMethod(p,n,a)*(b-a)/(fResMethod(p,n,b)-fResMethod(p,n,a));      
    if(fResMethod(p,n,a)*fResMethod(p,n,c)<0)
      a=c;
    else 
      b=c;
    k++;
  }
printf("METHOD ITERATION: kol-vo iteraciy k=%d, x=%lf, f(x)=%lf\n", k,c,fResMethod(p,n,c));    
}
            
int Read(Poly** p_1,char name[50])
{
    int i=0;
    Poly * p = *p_1;
    printf("Vvedite imya faila bez rasshireniya: ");
    scanf("%s", name, 49);              
    FILE *file = fopen(strcat(name,".txt"), "r");
     if (file == NULL)
        printf("Error opening\n");
    else
    {  
       while(!feof(file))
       {
 
        }
    }
    *p_1 = p;
    return i;
 
}
 
void Write(Poly*p, int n,char name[50],double x0,double x1,double y)
{
    int i=0,z=0,ptr;
    printf("Vvedite imya faila bez rasshireniya: ");
    scanf("%s", name, 49);              
    FILE *file = fopen(strcat(name,".txt"), "w");
    if (file == NULL)
        printf("Error opening\n");
    else
    { 
        for (i=n;i>=0;i--)
        { 
        if(i>1)
            fprintf(file,"(%g)*x^%d + ", p[i].koef, i);
        else
            fprintf(file,"(%g)%s",p[i].koef, i==1?"*x + ":"\n");
        }
        fprintf(file,"Na priblezhenii x1=%lf u x2=%lf->x=%lf",x0,x1,y);
    }
    fclose(file);
}
 
void WriteBin(Poly*p, int n,char name[50])
{   
    printf("Vvedite imya faila bez rasshireniya: ");
    scanf("%s", name, 49);              
    FILE *file = fopen(strcat(name,".txt"), "wb");
    if (file == NULL)
        printf("Error opening\n");
    else
    {
        fwrite (&n, sizeof(n),1, file);
        fwrite (p, sizeof(Poly),n, file);
        fclose(file);
    }
}
 
int ReadBin(Poly ** events_1,char name[50])
{
    int n = 0;
    Poly * p = *events_1;
    printf("Vvedite imya faila bez rasshireniya: ");
    scanf("%s", name, 49);              
    FILE *file = fopen(strcat(name,".txt"), "rb");
    if (file == NULL)
        printf("Error opening\n");
    else
    {
        fread (&n, sizeof(n),1, file);
        p=realloc(p, sizeof(Poly)*n);
        fread (p, sizeof(Poly), n, file);
    }
    *events_1 = p;
    return n;
}
 
 
void WriteKoren(char name[50],double x0,double x1,double y)
{
    int i=0,z=0,ptr;
    printf("Vvedite imya faila bez rasshireniya: ");
    scanf("%s", name, 49);              
    FILE *file = fopen(strcat(name,".txt"), "a");
    if (file == NULL)
        printf("Error opening\n");
    else
    { 
        fprintf(file,"\nNa priblezhenii x1=%lf u x2=%lf->x=%lf",x0,x1,y);
    }
    fclose(file);
}
 
int menu()
{
  int menu;
  printf("Menu:\n\n ");
  printf("Work with a polynomial->");
  printf("          1.  Enter polynomial\n");
  printf("                          2.  View  polynomial\n");
  printf("                          3.  Decide polynomial(cicle)\n");
  printf("                          4.  Decide polynomial(return x)\n\n");
  printf(" Write or read a polynomial->");
  printf("      5.  Write file(txt)\n");
  printf("                      6.  Write file(bin)\n");
  printf("                          7.  Write root to file(txt)\n");
  printf("                          8.  Read file(txt)\n");  
  printf("                          9.  Read file(bin)\n\n"); 
  printf(" Methods for solving equations->");
  printf("  10. List all methods\n");
  printf("                          11. Method Dihotomii\n");
  printf("                          12. Method Chord\n");
  printf("                          13. Method Sekush\n");
  printf("                          14. Method Iteracii\n");
  printf("                          15. Method Parabol\n");
  printf("                          16. Method Parabol-myillera\n");
  printf("                          17. Method Riddera\n");
  printf("                          18. Method Brendta\n");
  printf("                          19. Method Brendta-Dekkera\n\n");
  printf(" We are waiting for you again->"); 
  printf("      20.End program\n");
  scanf("%i", &menu);
  return menu;
}
 
 int main(void)
{
  int n=0;
  double a,b,x0,x1,y;
  char name[50];
  int i=0;  
  struct Poly*p=0;  
  char run=1;
  while(run) {
    switch(menu()) {
      case 1:printf("Vvedite stepen' polinoma\n");
            scanf("%d",&n); 
            p=malloc(sizeof(Poly)*n);
            Enter(p,n);
       break;
     case 2: Output(p,n);
       break;   
     case 3: Res(p,n);
            printf("Nachalo intervala proverki kornei: a=");
            scanf("%lf",&a);
            printf("Konec intervala  proverki kornei : b=");
            scanf("%lf",&b); 
       break;
      case 4: printf("Vvedite priblizhenie x1=");
              scanf("%lf",&x0);    //Ïåðâîå íà÷àëüíîå ïðèáëèæåíèå
              printf("Vvedite priblizhenie x2=");
              scanf("%lf",&x1);    //Âòîðîå íà÷àëüíîå ïðèáëèæåíèå 
              y=MethodKacatelinix(p,n,x0,x1);        
       break;
      case 5: Write(p,n,name,x0,x1,y);
       break;
      case 6: WriteBin(p,n,name);
       break;
      case 7: WriteKoren(name,x0,x1,y);
       break;
      case 8: n=Read(&p,name);
       break;
      case 9: n=ReadBin(&p,name);
       break;   
      case 10: dihotomia(p,n,a,b), methodChord(p,n,a,b),methodSekush(p,n,a,b),methodParabol(p,n,a,b), methodParabolMullera(p,n,a,b),methodRiddera(p,n,a,b),methodBrenta(p,n,a,b),methodDekkeraBrenta(p,n,a,b);
       break;
      case 11: dihotomia(p,n,a,b);
       break;
      case 12: methodChord(p,n,a,b);
       break;
      case 13: methodSekush(p,n,a,b);
      break;
      case 14: methodIteration(p,n,a,b);
      break;
      case 15: methodParabol(p,n,a,b);
      break;
      case 16: methodParabolMullera(p,n,a,b);
      break;
      case 17: methodRiddera(p,n,a,b);
      break;
      case 18: methodBrenta(p,n,a,b);
      break;
      case 19: methodDekkeraBrenta(p,n,a,b);
      break;
      case 20:
       return run=0;
      break;
    }
  }
  getch();
  return 0; 
}