Необходимо упрощение примитивных вычислений

clear; clc;
syms c p S A;
 
a = 3/1000;
b = 0.2;
d = exp(a/sqrt(3))/4;
n = 0.21;
p0 = 0.29;
 
 
%задаю х от 0 до 0.2 с шагом = 0.002
Xo = 0; Xmax = 0.2; t = 0.002;
x1 = (Xo:t:Xmax)';
 
%пара функций, которые будут в уравнениях
Fi(p) = (1+p/p0)^n;
Psi(p) = p0/(1+n)*(1+p/p0)^(n+1);    
 
%первое уравнение. Задача - получить зависимость S (x1).
B = log (4*S) == 2.*x1 + a.*double(Fi(x1))/sqrt(3);
 
for i = 1:101
D(i,1) = solve(B(i),S);
end
s = double(D);
 
%второе уравнение. Задача - получить зависимость A (S, x1).
C = log(4.*(s-A)).*log(4.*(s-A))/2/a + A./sqrt(s-A) == 3.*x1.*x1/4 + Psi(x1) - Psi(0) + a.*Fi(x1).*Fi(x1)/6 + a.*x1.*Fi(x1)/2;
for i = 1:101
H(i,1) = vpasolve(C(i),A);
end
a1 = double(H);
 
data = [x1, s, a1]
%В результате получаю таблицу: x1   S   A (хорошо бы в double

clear; clc;
%syms c p S A;
a = 3/1000;
b = 0.2;
d = exp(a/sqrt(3))/4;
n = 0.21;
p0 = 0.29;
 
%задаю х от 0 до 0.2 с шагом = 0.002
Xo = 0; Xmax = 0.2; t = 0.002;
x1 = (Xo:t:Xmax)';
 
%пара функций, которые будут в уравнениях
Fi=@(p) ((1+p/p0).^n);
Psi=@(p) ( p0/(1+n)*(1+p/p0).^(n+1));    
 
%первое уравнение. Задача - получить зависимость S (x1).
s=exp(2*x1 + a*Fi(x1)/sqrt(3))/4;
 
syms C A;
Z=3.*x1.*x1/4 + Psi(x1) - Psi(0) + a.*Fi(x1).*Fi(x1)/6 + a.*x1.*Fi(x1)/2;
%второе уравнение. Задача - получить зависимость A (S, x1).
%C = log(4.*(s-A)).*log(4.*(s-A))/2/a + A./sqrt(s-A) == Z;
C = (4.*(s-A)) == exp(sqrt((Z - A./sqrt(s-A))*2*a));
for i = 1:length(x1)
H(i,1) = vpasolve(C(i),A);
end
a1 = double(H);
 
data = [x1, s, a1]
%В результате получаю таблицу: x1   S   A (хорошо бы в double
 
%Проверка
z1=log(4.*(s-a1)).*log(4.*(s-a1))/2/a + a1./sqrt(s-a1);
z2=3.*x1.*x1/4 + Psi(x1) - Psi(0) + a.*Fi(x1).*Fi(x1)/6 + a.*x1.*Fi(x1)/2;
[z1, z2]