% Given data points t = [0 1 2 3 4 5]; x = [1.0 1.4 2.0 2.5 2.5]; y = [1.0 0.5 1.0 1.5 1.0]; % Fit natural cubic splines xx = linspace(0, 5, 100); % Evaluation points for spline sx = ncspline(t, x, xx); sy = ncspline(t, y, xx); % Plot the curve plot(xx, sx, xx, sy, t, x, 'o', t, y, 'o'); axis equal; grid on; legend('x(t)', 'y(t)', 'Data points for x(t)', 'Data points for y(t)')