Решение линейного уравнения

clc
clear all
a1 = 9.1; 
a2 = 34.4;
a3 = 25.6;
a4 = 30.2;
P = 14.8;
k = 1:0.01:30;
x1 = length(k);
x2 = length(k);
x3 = length(k);
x4 = length(k);
x1>=1 && 8<=x1; %пределы изменения х1...х4
x2>=12 && 21<=x2;
x3>=14 && 24<=x3;
x4>=18 && 27<=x4;
p1 = a1*x1;
p2 = a2*x2;
p3 = a3*x3;
p4 = a4*x4;
P_m = p1+p2+p3+p4;
 
while (P-P_m)^2<=5 %цикл должен считать до тех пор, пока невязка не станет меньше 5
    x1=x1+0.01;
    x2=x2+0.01;
    x3=x3+0.01;
    x4=x4+0.01;
end
 
x = [x1 x2 x3 x4]

clear, clc
 
a = [11.50; 40.80; 23.50; 24.50];
P = 14.56;
 
N = 1e5; % размер выборки случайных чисел
x(:,1) = 1 + 7*rand(N,1); % от 1 до 8
x(:,2) = 12 + 6*rand(N,1); % от 12 до 18
x(:,3) = 15 + 6*rand(N,1); % от 15 до 21
x(:,4) = 16 + 7*rand(N,1); % от 16 до 23
Y = x*a; % сумма произведений для каждой строки
 
[r i] = min( abs(Y-P) ) ; % индекс ближайшего значения
x(i,:) % значения иксов
Y(i) % значение суммы

clear, clc
 
a = [0.115; 0.408; 0.235; 0.2450];
P = 14.56;
 
F = @(x) abs( x*a-P ); % функция невязки
 
x0 = [4.5, 15, 18, 19.5]; % начальное значение х для поиска
x = fminsearch(F, x0) % находим значения х для минимальной невязки
 
x*a - P % невязка