1. Составить программу расчета.
2. Произвести расчет интегральной кривой, начиная с x0 ( начальная точка ) до x0+1, с шагом h = 1/n, n = 2,4,8,16,…,N, где N – число, при котором разница между численным решением и аналитическим не превышает 0,0001.

Система уравнений
dy/dx=z-3y
dz/dx=z-y
Начальные условия x=1, y=z=1

Аналитическое решение: у=(A-B+Bx)e^-2x
z=(-A-Bx)e^-2x
A= -1; B= -2.