p> t1 = t1 + h/2;
y1 = y + ((h/2) * z1);
z1 = z1 + ((h/2) * (e * pow ((n0/h), 1/2)-b1 * z-b2 * fabs (z) * z-c1 * y-c3 * pow (y, 3 )));
mas1 [i] = y1;
i = i +1;
}
/* close the file */
// fclose (stream);
y_rez = 0;
y_max = mas1 [0];
for (i = 0; i <= 119; i + +)
{
y_rez = y_rez + (pow ((mas1 [i]-mas [i]), 2));
if (Mas1 [i]> y_max)
y_max = mas1 [i];
}
eps = pow (y_rez/n, 0.5)/fabs (y_max);
printf ("% f epsilon for our operations n ", eps);
fprintf (stream, "% f epsilon for our operations n", eps);
fclose (stream);
В
return 0;
}
br/>