Табулирование функции через массив - C++ - Ответ 7455191

#include <stdio.h>
#include <iostream>
#include <cmath>
#include <iomanip>
#include <valarray>
using namespace std;
double f(double u) {
return sqrt(2*u)+2*u-6.3;
}
int fz(double d,double k,double z) {
return d*pow(z, 4)+k*log(pow(z, 2)+pow(k,2));
}
int main() {
setlocale(LC_ALL, "Russian");
double l, r, c, eps,b, a, xn, xk, dx, x[100], y[100];
int j,N;
cout << "Введите l:\n"; //решение нелинейного уравнения, начало
    cin >> l;
    cout << "Введите r:\n";
    cin >> r;
    cout << "Введите eps: \n";
    cin >> eps;
while( fabs(f(r) - f(l)) > eps ) {
c = ( l + r ) / 2;
if( f(c) * f(r) < 0 )
l = c;
else
r = c;
}
b=((l + r)/2);
cout << fixed << setprecision(2) << "Корень=" << b << "\n\n"; //решение уравнение, конец. работает отлично.
cout << "Введите a:\n"; //табулирование, начало. работает на получение НЕверных результатов.
    cin >> a;
    cout << "Введите xn:\n";
    cin >> xn;
    cout << "Введите xk:\n";
    cin >> xk;
    cout << "Введите dx:\n";
    cin >> dx;
    j=1;x[j]=xn;
    do
    {
    fz(a,b,x[j]);
    j++;
    x[j]=x[j-1]+dx;
    }
    while (x[j]<=xk);
    cout << "Результаты: \n";
    for (int j = 1; j < N; j++)
    cout<< fixed << setprecision(3) << "x[" << j <<"]=" << x[j] <<' ' << "y[" << j <<"]=" << y[j] <<"\n";
cin.get();
}