R1 =  14980
R3 =  14980
R4 =   1440
Rv =   8250
Rr = 100000
UB = 5.10

UqE = R4/(R3+R4)*UB
RE = R3 * R4 / (R3+R4)

R25 = 2000
a = 7.88e-3
b = 1.937e-5


U_a = (0:5)'
U_p = ((RE + Rv) * U_a + Rr * R4 / (R3 + R4) * UB) / (RE + Rv + Rr)
R_T = U_p ./ (UB - U_p) .* R1
T = ( 25 + (sqrt(a .* a - 4 .* b + 4 .* b .* R_T / R25) - a ) ./ ( 2 .* b ) )


U_a = (0:500)';
U_p = ((RE + Rv) .* (U_a ./ 100) + Rr * R4 / (R3 + R4) * UB) / (RE + Rv + Rr);
R_T = U_p ./ (UB - U_p) .* R1;
T = ( 25 + (sqrt(a .* a - 4 .* b + 4 .* b .* R_T / R25) - a ) ./ ( 2 .* b ) );



m=20.987
b=-25.462

Terr = (U_a ./ 100) .* m + b


###
m = 0.047648
b = 1.2132



octave:52> U_a = (0:5)'
U_a =

  0
  1
  2
  3
  4
  5

octave:53> U_a = (0:5)';
octave:54> U_p = ((RE + Rv) * U_a + Rr * R4 / (R3 + R4) * UB) / (RE + Rv + Rr)
U_p =

  0.40822
  0.49551
  0.58280
  0.67009
  0.75737
  0.84466

octave:55> R_T = U_p ./ (UB - U_p) .* R1
R_T =

  1303.4
  1612.1
  1932.7
  2265.9
  2612.6
  2973.5

octave:56> T = ( 25 + (sqrt(a .* a - 4 .* b + 4 .* b .* R_T / R25) - a ) ./ ( 2 .* b ) )
T =

  -25.4620
   -1.3183
   20.6824
   41.2266
   60.7312
   79.4735

octave:57> 

octave:61> I = UB ./ (R1 + R_T )
I =

  0.00031320
  0.00030738
  0.00030155
  0.00029572
  0.00028989
  0.00028407


m=0.047648317
b=1.21322146





R1 =  14980
R3 =  14980
R4 =   1440
Rv =   8250
Rr = 100000
UB = 5.10

UqE = R4/(R3+R4)*UB
RE = R3 * R4 / (R3+R4)

R25 = 2000
a = 7.88e-3
b = 1.937e-5



octave:1> load "daylen.onlynum.dat"
octave:2>         a = 818.869776162182
a = 818.87
octave:3>         b = 284.650966179578
b = 284.65
octave:4>         c = 0.0161339421587499
c = 0.016134
octave:5>         d = 0.33668118331649
d = 0.33668
octave:6> daylen2 = [daylen * [1, 0]', daylen * [0, 1]', a + b .* cos(c * daylen * [1, 0]' + d)]
octave:7> save "daylen.onlynum.withcos.dat" daylen2
octave:8> daylen3 = [daylen * [1, 0]', daylen2 * [0, 1, 0]' - daylen2 * [0, 0, 1]']
octave:11> save "daylen.onlynum.withcos.diff.dat" daylen3

octave:8> daylen4 = [daylen * [1, 0]', (100 .* daylen3 * [0, 1]' ) ./ (daylen2 * [0, 0, 1]')]
octave:11> save "daylen.onlynum.withcos.diff.dat" daylen4