kauchy

#include <iostream>
#include <math.h>
#include <fstream>
#include <cmath>

using namespace std;
const double right_ending = 1;
const double left_ending = 0;
const double eintial_cond_func = 1;\\in 0
const double inteial_cond_deliv = 3;\\in 0
const int accuracy = 10000;
int main (){
    double step = (right_ending - left_ending)/accuracy;
    double n = 1,2;
    double tetha;
    double q;
    double xi;
    double solution_1[accuracy]; \\deliv_array: teta
    double solution_2[accuracy]; \\func_array: q
    double variable[accuracy];
    xi = left_ending + 0.001 * step;
    solution_1[0] = intial_cond_func;
    solution_2[0] = intial_cond_deliv;
    variable[0] = xi;
    for (int i = 0, i = accuracy - 1, i++){
        solution_1[i+1] = solution_1[i] + step * solution_2[i];
        solution_2[i+1] = solution_2[i] + step * (-2 / xi * solution_2[i] - solution_1[i]^n);
        xi = xi + step;
        variable[i+1] = xi;
        }
\\ решение диф - solution_1(variable)
    return 0;
}