{$N+}Program M;

CONST
     n=4;
TYPE
    vektor=array[1..n+1] of real;

VAR
   i,j,k:integer;
   a,b,c,x,rho,mu,f:vektor;
   filein:text;
BEGIN

ASSIGN(filein,'data3diag.dat');
RESET(filein);

FOR i:=1 TO n DO
BEGIN
    READ(filein,c[i]);
    READ(filein,a[i]);
    READ(filein,b[i]);
    READ(filein,f[i]);
    rho[i]:=0;
    mu[i]:=0;
    x[i]:=0;
END;

x[n+1]:=0;
rho[n+1]:=0;
mu[n+1]:=0;

FOR i:=1 TO n DO
BEGIN
    mu[i+1]:=-b[i]/(c[i]*mu[i]+a[i]);
    rho[i+1]:=(f[i]-c[i]*rho[i])/(c[i]*mu[i]+a[i]);
END;

FOR i:=n DOWNTO 1 DO
    x[i]:=mu[i+1]*x[i+1]+rho[i+1];

FOR i:=1 TO n DO
    WRITELN(x[i]:7:3);

CLOSE(filein);

READLN;

END.