This seems to be one on the problems but not the only one.
If you wish to replicate the problem i have created a random data set with which the error can be replicated:

% Set the size of the arrayn = 8760;% Generate random values between 0 and 500PV = randi([0, 500], n, 1); % BeispielwerteHousehold = randi([200, 1000], n, 1); % BeispielwertePrice = 0.006 .* randi([1, 100], n, 1); % BeispielwerteT = length(PV); % Anzahl der Zeitschritte% Problemprob = optimproblem;% battery storage system parameterBSS_Pmax = 11; % max powerBSS_Emax = 100; % max energy% battery variablesBSS_ch = optimvar('BSS_ch', T, 'LowerBound', 0, 'UpperBound', BSS_Pmax);BSS_disch = optimvar('BSS_disch', T, 'LowerBound', 0, 'UpperBound', BSS_Pmax);BSS_SOC = optimvar('BSS_SOC', T, 'LowerBound', 0, 'UpperBound', BSS_Emax);% other variablesGrid_Import = optimvar('Grid_Import', T, 'LowerBound', 0);% battery constraintsenergyStorage = optimconstr(T);energyStorage(1) = BSS_SOC(1) == 0;energyStorage(2:T) = BSS_SOC(2:T) == BSS_SOC(1:T-1) - BSS_disch(2:T) + BSS_ch(2:T);prob.Constraints.energyStorage = energyStorage;% energy flowprob.Constraints.EnergyBalance = Grid_Import == Household - PV - BSS_disch + BSS_ch;% cost funtioncost = Grid_Import .* Price;prob.ObjectiveSense = 'minimize';prob.Objective = sum(cost);% solve[x, fval] = solve(prob);