differentiate v(t) = v0 * (1-(t/t).^4).^2 when t = 20 t = 50 v0 = 9.7 in matlab

main.m
% Given parameters
v0 = 9.7;
t = linspace(20, 50); % create a vector of t values between 20 and 50

% Define the function
v = @(t) v0.*(1 - (t./t).^4).^2; % v(t) function

% Differentiate the function
dv = @(t) -16*v0*(t.^3).*(1 - (t./t).^4); % v'(t) function

% Plot the function and its derivative
hold on % to plot both functions in same figure
plot(t, v(t), 'LineWidth', 2) % plot v(t)
plot(t, dv(t), 'LineWidth', 2) % plot v'(t)

% Add labels and legend
xlabel('t')
ylabel('v(t), v''(t)')
legend('v(t)', 'v''(t)')
519 chars
20 lines