Ошибка в решении нелинейной системы уравнений

syms x y a b c double
Solution = solve(exp(6*a + b)+ c - 0.25,exp(3*a + b)+ c - 1.25,exp(2*a + b)+ c - 1.75,a,b,c);
A = double(Solution.a);
B = double(Solution.b);
C = double(Solution.c);

clear, clc
% x(1)->a, x(2)->b, x(3)->c
F = @(x) [exp(6*x(1) + x(2)) + x(3) - 0.25
    exp(3*x(1) + x(2)) + x(3) - 1.25
    exp(2*x(1) + x(2)) + x(3) - 1.75];
x0 = [1 0 0]; % начальное предположение
X = fsolve(F,x0);
a = X(1)
b = X(2)
c = X(3)

T=[0 0.6]
Q=[6.9*T + 1.24]
plot(T,Q);
axis([0 24 0 6])
 
T=[0.6 12];
%W=[exp(-0.6844*T + 2.0454) + 0.1326];
plot(T,(exp(-0.6844*T + 2.0454) + 0.1326));
axis([0 24 0 6])
 
T=[12 13]
E=[15.8667*T - 191.3667]
plot(T,E);
 
T=[13 24]
R=[exp(-0.0695*T - 2.1279) + 1.7908]
plot(T,R);
axis([0 24 0 6])

clear, clc
 
T=0:0.01:0.6;
Q=6.9*T + 1.24;
plot(T,Q);
hold on
 
T=0.6:0.01:12;
W=exp(-0.6844*T + 2.0454) + 0.1326;
plot(T,W);
 
 
T=12:0.01:13;
E=15.8667*T - 191.3667;
plot(T,E);
 
T=13:0.01:24;
R=exp(-0.0695*T - 2.1279) + 1.7908;
plot(T,R);
axis([-1 24 0 6])

clear, clc
 
t = 0:0.01:10;
s = double(sin(t)>0); % прямоугольный сигнал
h = exp(-2*t); % переходная функция
S = abs(diff(s)); % модуль производной сигнала
Y = conv(S,h); % делаем свертку
plot(t,s,':b',t,Y(1:length(t)),'r')
axis([-1 10 -0.1 1.1])

clear, clc
 
T=0:0.01:0.6;
Q=6.9*T + 1.24;
plot(T,Q);
hold on
 
T=0.6:0.01:12;
W=exp(-0.6844*T + 2.0454) + 0.1326;
plot(T,W);
 
T=12:0.01:13;
E=5.1261*T - 61.3789 ; % причесал
plot(T,E);
 
T=13:0.01:24;
R=exp(-0.6844*(T-13+0.6) + 2.0454) + 0.1326; 
% или если раскрыть скобки:
% R=exp(-0.6844*T+10.5324) + 0.1326; 
plot(T,R);
axis([-1 24 0 6])

clear ,clc
 
T1 = 0:0.01:24;
T2 = 0:0.4:24;
Y = @(t) (6.9*t+1.24).*((t>=0)&(t<=0.6)) + ... % я ленивая попа :3
    (exp(-0.6844*t+2.0454)+0.1326).*((t>0.6)&(t<=12)) + ...
    (5.1261*t-61.3789).*((t>12)&(t<=13)) + ...
    (exp(-0.6844*t+10.5324)+0.1326).*(t>13);
 
plot(T1,Y(T1),'b', T2,Y(T2),'or');
axis([-1 24 0 6])