blob: 85e50c8431bb8e47d1c9ec946a9e161a204acb7e (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
|
Tmax=200
Tstart=80
p=poly([1,-2*Tmax/Tstart,1],"smaxT","coeff")
smaxT=roots(p)
smaxT=smaxT(2)
disp(smaxT)
p=poly([1,-4,1],"w","coeff")
w=roots(p)
sfl=smaxT/w(1)
disp(sfl)
ratio=sqrt(((smaxT/sfl)^2+1)/(smaxT^2+1))
disp(ratio)
p=poly([1,-((smaxT/sfl)^2+1)*sfl,smaxT^2],"k","coeff")
k=roots(p)
k=k(2)
disp(k)
|
