From a7dd19f779aa1f1f785a144d641a974484eb648c Mon Sep 17 00:00:00 2001
From: margaret-k <60312043+margaret-k@users.noreply.github.com>
Date: Mon, 25 May 2020 23:57:54 +0300
Subject: [PATCH] Create dz15

---
 dz15 | 136 +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
 1 file changed, 136 insertions(+)
 create mode 100644 dz15

diff --git a/dz15 b/dz15
new file mode 100644
index 0000000..f573dd3
--- /dev/null
+++ b/dz15
@@ -0,0 +1,136 @@
+#include <stdio.h>
+#include <stdlib.h>
+#include <ctype.h>
+
+float* read(char *name1){ 
+  FILE *f1;
+  f1 = fopen(name1, "rb");
+  float *r1 = (float *) malloc(sizeof(float));
+  int n1 = 0;
+  while(!feof(f1))
+  {
+    float c;
+    fread(&c, sizeof(float), 1, f1);
+    //printf("%f", c);
+    r1[n1] = c;
+    n1++;
+    r1 = (float *) realloc(r1, sizeof(float));}
+    fclose(f1);
+return r1;
+}
+void input(char *name){ 
+  FILE *f = fopen(name, "wb");
+  int n;
+  printf("n = ");
+  scanf("%d", &n);
+  for(int i = 0; i <= n; i++){
+    printf("Koef pry x^%d = ", i);
+    float g;
+    scanf("%f",&g);
+    float x[2] = {g, (float) i};
+    fwrite(&x, sizeof(float), 2, f);
+  }
+  fclose(f);
+}
+
+void output(char *name){ 
+  float *r = read(name);
+  int n = sizeof(r)/sizeof(float);
+    for(int i = n - 2; i > 3; --i){
+      if(i % 2 == 1){printf("%f*x^%d + ",r[i - 1],(int) r[i]);}
+    }
+    printf("%f*x + %f = 0", r[2] ,r[0]);
+  free(r);
+}
+void add(char *name1, char *name2, char *res){ 
+  float *r1 = read(name1);
+  float *r2 = read(name2);
+  int n1 = sizeof(r1)/sizeof(float);
+  int n2 = sizeof(r2)/sizeof(float);
+  int n = n1<n2?n1:n2;
+  FILE *re = fopen(res, "wb");
+  for(int i = 0; i <= n; ++i){
+    if (i % 2 == 0){
+    float x[2] = {r1[i] + r2[i], (float)i + 1};
+    printf("%f %f\n",r1[i] + r2[i], (float) i+1  );
+    fwrite(&x, sizeof(float), 2, re);}
+  }
+  if(n1 > n2){
+  for(int i = n + 1; i < n1; ++i ){
+    if (i % 2 == 0){
+    float x[2] = {r1[i], (float) i + 1};
+    fwrite(&x, sizeof(float), 2, re);}
+  }
+}
+else{
+  for(int i = n + 1; i < n1; ++i ){
+    if (i % 2 == 0){
+    float x[2] = {r2[i],(float) i + 1};
+    fwrite(&x, sizeof(float), 2, re);}
+}
+}
+fclose(re);
+free(r1); free(r2);
+}
+
+void sub(char *name1, char *name2, char *res){ 
+  float *r1 = read(name1);
+  float *r2 = read(name2);
+  int n1 = sizeof(r1)/sizeof(float);
+  int n2 = sizeof(r2)/sizeof(float);
+  int n = n1<n2?n1:n2;
+  FILE *re = fopen(res, "wb");
+  for(int i = 0; i <= n; ++i){
+    if (i % 2 == 0){
+    float x[2] = {r1[i] - r2[i], (float)i + 1};
+    printf("%f %f\n",r1[i] - r2[i], (float) i+1  );
+    fwrite(&x, sizeof(float), 2, re);}
+  }
+  if(n1 > n2){
+  for(int i = n + 1; i < n1; ++i ){
+    if (i % 2 == 0){
+    float x[2] = {r1[i], (float) i + 1};
+    fwrite(&x, sizeof(float), 2, re);}
+  }
+}
+else{
+  for(int i = n + 1; i < n1; ++i ){
+    if (i % 2 == 0){
+    float x[2] = {-r2[i],(float) i + 1};
+    fwrite(&x, sizeof(float), 2, re);}
+}
+}
+fclose(re);
+free(r1); free(r2);
+}
+
+void mult(char *name1, char *name2, char *res){ 
+  float *r1 = read(name1);
+  float *r2 = read(name2);
+  int n1 = sizeof(r1)/sizeof(float);
+  int n2 = sizeof(r2)/sizeof(float);
+  float r[2*n1*n2];
+  FILE *re = fopen(res, "wb");
+  for(int i = 0; i <=n1; ++i){
+    for(int j = 0; j <= n2; ++j){
+      if(i%2 == 0 && j%2 == 0){
+        r[i+j] += r1[i] * r2[j];
+        r[i+j+1] = i+j;
+      }
+    }
+  }
+  for(int i = 0; i <= 2*n2*n1; ++i ){
+    if (i % 2 == 0){
+    float x[2] = {r[i],(float) i + 1};
+    fwrite(&x, sizeof(float), 2, re);}
+}
+fclose(res); free(r1); free(r2); free(r3);
+}
+int main(){
+  char *name1 = "1.dat";
+  char *name2 = "2.dat";
+  char *res = "res.dat";
+  input(name1);
+  output(name1);
+
+}