In[]:=
Clear["`*"];n=3;t=0.45;tj={0.35,0.36,0.37};ρs={0.85,0.85,0.85};ρl={0.50,0.45,0.4};ρ=ρsρl;b=1.65;c=0.1;α=0.4;sigma[i_]:=(t-tj[[i]]);ω1={{ω11,ω12,ω13},{ω21,ω22,ω23},{ω31,ω32,ω33}};omega[i_,j_]:=Ifij,0,;ω=Array[omega,{n,n}];ω1//MatrixFormω11=ω[[1]][[1]];ω12=ω[[1]][[2]];ω13=ω[[1]][[3]];ω21=ω[[2]][[1]];ω22=ω[[2]][[2]];ω23=ω[[2]][[3]];ω31=ω[[3]][[1]];ω32=ω[[3]][[2]];ω33=ω[[3]][[3]];ω1//MatrixFormω//MatrixForm
(1-ρ[[i]])ρ[[i]]
;σ=Array[sigma,n];λ=1
α[[i]]
n
∑
i=1
t-tj[[i]]
2
σ
λ(t-tj[[i]])
2(c-λ(t-tj[[i]]))
Out[]//MatrixForm=
ω11 | ω12 | ω13 |
ω21 | ω22 | ω23 |
ω31 | ω32 | ω33 |
Out[]//MatrixForm=
0 | 0.00885475 | 0.00885475 |
0.00795519 | 0 | 0.00795519 |
0.0070588 | 0.0070588 | 0 |
Out[]//MatrixForm=
0 | 0.00885475 | 0.00885475 |
0.00795519 | 0 | 0.00795519 |
0.0070588 | 0.0070588 | 0 |