проблема в формуле,никак не могу разобраться в чем дело

[ксюшенькаа]

вот листинг, я уже изначально задала матрицу, чтоб ее не вводить, нужен результат, Очень нужно, но мозг отказывается понимать в чем он не прав^^


Const t=4; k=6;
O: array [1..k] of real=
(9.5, 120, 6, 2, 70, 900);
M: array [1..k] of real=
(10.5, 140, 4.5, 15, 110, 650);
S: array [1..t,1..k] of real= ((9.6, 126, 5.8, 5.9, 74, 875), (9.61, 127, 5.81, 5.91, 75, 876), (9.62, 129, 5.82, 5.92, 76, 877), (9.63, 129, 5.83, 5.93, 77, 878));
v: array [1..k] of real= (0.05, 0.3, 0.15, 0.5, 0.2, 0.1);
Type dmyarray=array [1..t,1..k] of real;
Var i,j,n: byte;
x1,x2,D,a,b,c,Ql1,Ql2,Qnl1,Qnl2,Qml 1,Qml2,vn11, vn21, vn31, vn22, vn32, vn33, m33, m23, m13, m03, m22, m12, m02, m11, vn11n, vn21n, vn31n, vn22n, vn32n, vn33n, m33n, m23n, m13n, m03n, m22n, m12n, m02n, m11n, vn11m, vn21m, vn31m, vn22m, vn32m, vn33m, m33m, m23m, m13m, m03m, m22m, m12m, m02m, m11m: real;
E: dmyarray;
Ql: array [1..t] of real;
Qnl: array [1..t] of real;
Qml: array [1..t] of real;
Begin
Ql1:=0;
Ql2:=0;
Qnl1:=0;
Qnl2:=0;
Qml1:=0;
Qml2:=0;
For i:=1 to t do
begin
For j:=1 to k do
begin
E[i,j]:=(M[j]-O[j])/(O[j]-S[i,j]);
Ql1:=Ql1+v[j]*E[i,j];
Ql2:=Ql2+v[j];
Ql[i]:=Ql1/Ql2;
Qnl1 := Qnl1+Exp(1/v[j]*Ln(E[i,j]));
Qnl2:=Qnl2+v[j];
Qnl[i]:=exp(Qnl2*ln(Qnl1/6));
Qml1:=Qml1+Exp((1/v[j])*ln(E[i,j]));
Qml2:=Qml2+v[j];
Qml[i]:=1/6*Qml1/Qml2;
end;
end;
vn11:=Ql[4]-Ql[3];
vn21:=vn11-Ql[2];
vn31:=vn21-Ql[1];
vn22:=vn11-vn21;
vn32:=vn21-vn31;
vn33:=vn22-Ql[1];
m33:=1/6*vn33;
m23:=1/2*vn22+1/2*vn33;
m13:=vn11+1/2*vn22+1/3*vn33;
m03:=Ql[4];
m22:=1/2*vn22;
m12:=vn11+1/2*vn22;
m02:=Ql[4];
m11:=-Ql[4]/vn11;
Writeln ( m33);
Writeln (m23);
Writeln (m13);
Writeln (m03);
Writeln (m22);
Writeln (m12);
Writeln (m02);
Writeln ( m11);
vn11n:=Qnl[4]-Qnl[3];
vn21n:=vn11n-Qnl[2];
vn31n:=vn21n-Qnl[1];
vn22n:=vn11n-vn21n;
vn32n:=vn21n-vn31n;
vn33n:=vn22n-Qnl[1];
m33n:=1/6*vn33n;
m23n:=1/2*vn22n+1/2*vn33n;
m13n:=vn11n+1/2*vn22n+1/3*vn33n;
m03n:=Qnl[4];
m22n:=1/2*vn22n;
m12n:=vn11n+1/2*vn22n;
m02n:=Qnl[4];
m11n:=-Qnl[4]/vn11n;
Writeln (m33n);
Writeln (m23n);
Writeln (m13n);
Writeln (m03n);
Writeln ( m22n);
Writeln (m12n);
Writeln (m02n);
Writeln ( m11n);
vn11m:=Qml[4]-Qml[3];
vn21m:=vn11m-Qml[2];
vn31m:=vn21m-Qml[1];
vn22m:=vn11m-vn21m;
vn32m:=vn21m-vn31m;
vn33m:=vn22m-Qml[1];
m33m:=1/6*vn33m;
m23m:=1/2*vn22m+1/2*vn33m;
m13m:=vn11m+1/2*vn22m+1/3*vn33m;
m03m:=Qml[4];
m22m:=1/2*vn22m;
m12m:=vn11m+1/2*vn22m;
m02m:=Qml[4];
m11m:=-Qml[4]/vn11m;
Writeln (m33m);
Writeln (m23m);
Writeln (m13m);
Writeln (m03m);
Writeln (m22m);
Writeln (m12m);
Writeln (m02m);
Writeln ( m11m);
end.

[ксюшенькаа]

проблема где то здесь, т.к. результат NAN-а это не то что ожидается

Qnl1 := Qnl1+Exp(1/v[j]*Ln(E[i,j]));
Qnl2:=Qnl2+v[j];
Qnl[i]:=exp(Qnl2*ln(Qnl1/6));
Qml1:=Qml1+Exp((1/v[j])*ln(E[i,j]));
Qml2:=Qml2+v[j];
Qml[i]:=1/6*Qml1/Qml2;