syms t
PHI=[ 1, -t, -t/3 – (2*exp(-3*(-t)))/9 + 2/9, (2*-t)/3 + (2*exp(-3*(-t)))/9 – 2/9;
0, 1, (5*exp(-3*(-t)))/12 – (3*exp(-t))/4 + 1/3, 2/3 – exp(-t)/4 – (5*exp(-3*(-t)))/12;
0, 0, exp(-3*(-t))/4 + (3*exp(-t))/4, exp(-t)/4 – exp(-3*(-t))/4;
0, 0, (3*exp(-t))/4 – (3*exp(-3*(-t)))/4, (3*exp(-3*(-t)))/4 + exp(-t)/4];
PHIT=transpose (PHI);
B=[0;1;2;1];
BT=transpose (B);
GRAMi = PHI*B*BT*PHIT
GRAMfinal=int(GRAMi,t, 0, t)
A= det(GRAMfinal) %% Absolute determinant of the matrix
N = @(t) (A(t));
t = linspace(0, 1, 50);
for k = 1:numel(t)
Nt(k) = N(t(k));
end
figure
plot(t, Nt)
grid
Best Answer