
clc;
clear all;

Ts = 600/3600;
E = 100;


soc_0 = @(eta, P) 0.5 + (1/eta - eta)*P*Ts/2/E;

P = 100;
eta = 0.5;
s0 = soc_0(eta, P);

% verification
1 - (s0 + eta * P * Ts/E)
s0 - 1/eta * P * Ts/E

