вычисление производных первого и второго порядка, используя численные методы (С++ или другой язык)

[mari2853]

за ночь надо сделать задачу на С++, хотя в принципе это не особо важно, хоть на чем. Суть- вычисление производных первого и второго порядка, используя численные методы(теория есть). переделываю проект с другой задачей на свою. может кто-то сталкивался с подобным???
навояла за ночь чисто по логике задачу. конечно код не очень, но нужно было очень срочно...
1. проблема- надо все как-то привести к общему виду и запуску
2. проблема с формулой, где в одной операции необходим цикл, как это можно сделать??

#include <stdio.h>
#include <math.h>
#include <iomanip>
#include <iostream>  //потоковый ввод-вывод (cin,cout)
#include <conio.h>     //работа с консолью (clrscr,kbhit)
#include <math.h>  
using namespace std;
int n; int a;
double a[10];
double y[10][10];
double f(double x)
   
void main()

{
  cout<<"Program 1";
  cout<<"\n\nSource data:\n";
  char s[10];
 { 
    FILE *f1=fopen("data.txt","r");    //в этой части читаю первую строку с файла
     fscanf(f1,"%s",n);
       n=atof(s);
       cout<<""<<n<<endl;

  }
   
  double a[20];
  for (int i=0; i<=n; i++)
  {
     FILE *f1=fopen("data.txt","r+");   ////считываю опять все сначала и заношу все в одномерную матрицу
      fscanf(f1,"%s",a);
      a=atof(s);
      cout<<""<<a<<endl;
      fclose(f1);
  }

  double y[20];   int k=1; int m=n;
  {do
       {
       for (int i=1; i<n; i++)              ////преобразования с матрицей  
          for (int j=m+1; j<=m; j++)
         y[k]=a[j+1]-a[j];
       a[j]=y[k];

         }
     n--;
     k++
  }
   while (n>1)     //вот с этими формулами проблема, она должна быть зациклена на увеличение , пока проставила цифрами, именно к эт примеру подходит
      }
int i=1; int j=1;
pr1= 1/(a[i+1]-a)*((y[j])-1/2*a[j+1]+1/3*a[j+2]-1/4*a[j+3];
pr2= 1/((a[i+1]-a)*(a[i+1]-a))*((a[j+1]-a[j+2]+11/12*a[j+3];
    

  printf("\nTable konechnych raznostej:\n");    //здесь бы все аккуратненько вывести в файл, но нифига не помню, как это сделать
  printf("%5s %10s %10s  %10s %10s\n","y", "y^2","y^3","y^4");
  printf("%5d %10.6f %10.6f\n",y[j]);

  fprintf(f2,"\nTable proizvodnoj:\n");
  fprintf(f2,"%5s %10s %10s\n","1","2");
  fprintf(f2,"%5d %10.6f %10.6f\n",pr1, pr2);

  cout<<"\nResults:\n";
  cout<<"\nInt"<<Int;
  cout<<"   ("<<n<<" iterations)";
    //записываем результат в файл result.txt
  fprintf(f2,"Result\n");
  fprintf(f2,"Int(trap)=%f   (%d iterations)\n",Int,n);
  fclose(f2);
  cout<<"\n\nPress ENTER for exit";
    //ждем нажатия любой клавиши для завершения программы
  }

входной файл из 11 строк, в каждой по одному числу. ПОМОГИТЕ!!!!!


___________
1) Название темы должно адекватно отражать суть решаемой задачи/проблемы.

2) Код нужно оформлять по правилам:
тегом [CODE]..[/СODE] (это кнопочка с решёточкой #)

На первый раз я исправил.

В дальнейшем темы с подобным названием будут закрываться/удаляться,
а автор такой темы получать штрафы.

Прошу учесть на будущее...

Модератор.