bild26.m

function bild26()
% Matrixwerte
m1 = 6.0e3; m2 = 6.0e3; m3 = 1.0e3;
c1 = 3.0e6; c2 = 3.0e6; c3 = 1.0e6;
M = diag([m1 m2 m3]);
C = [c1 + c2, -c2, 0; -c2, c2 + c3, -c3; 0, -c3, c3];

% Anfangsbedingungen
x0 = zeros(3,1);
v0 = [0; 0; 10];
y0 =  [x0; v0];

% DGL lösen
tEnd = 2;
myfunc = @(t, y) schwingNd(t, y, M, C);
[t, y] = ode45(myfunc, [0 tEnd], y0);

% Ergebnis plotten
plot(t, y(:,1), t, y(:,2), t, y(:,3));
legend("x_1", "x_2", "x_3", "Location", "best");

F = getframe(gcf);
imwrite(F.cdata, "bild26.png");