function B9
clc;
clear;
close all;

%Phần đầu (Khai báo biến, hàm...)------------------------------------------
%Khai báo biến symbolic t
syms t;

%Khai báo các hàm x và y theo t
x = 3 * t.^2 - 4/3 * t.^3;
y = 8 * t;

%Tìm vận tốc theo phương x và y = đạo hàm của x, y theo t
vx = diff(x);
vy = diff(y);

%Tìm gia tốc theo phương x và y = đạo hàm của vx, vy theo t
ax = diff(vx);
ay = diff(vy);

%Phần sau (Thực hiện yêu cầu các câu hỏi)----------------------------------

%Nhập khoảng thời gian cần vẽ đồ thị:
disp('Nhap khoang thoi gian can ve do thi');
t1 = input('Moc thoi gian dau tien: t1 = ');
t2 = input('Moc thoi gian thu hai: t2 = ');

%----------------Hết phần nhập dữ liệu vẽ đồ thị câu a--------------------%

%Câu a: Vẽ đồ thị x, y theo t với giá trị t từ 0 tới 5

%Đồ thị 1: Đồ thị x, y ----------------------------------------------------
subplot(3,1,1)
fplot(x, y, [t1 t2]);

%Ký hiệu các trục x, y và tiêu đề
title('Quy dao cua vat trong khoang thoi gian tu t1 den t2');
xlabel('x');
ylabel('y');

%Hiển thị lưới
grid on;

%Vẽ 2 trục tọa độ
xline(0);
yline(0);

%Đồ thị 2: Đồ thị x, t ----------------------------------------------------
subplot(3,1,2);
fplot(t, x, [t1 t2]);

%Ký hiệu các trục và tiêu đề
title('Hoanh do x cua vat trong khoang thoi gian tu t1 den t2');
xlabel('t');
ylabel('x');

%Hiển thị lưới
grid on;

%Vẽ 2 trục tọa độ
xline(0);
yline(0);

%Đồ thị 3: Đồ thị y, t ----------------------------------------------------

subplot(3,1,3);
fplot(t, y, [t1 t2]);

%Ký hiệu các trục và tiêu đề
title('Tung do y cua vat trong khoang thoi gian tu t1 den t2');
xlabel('t');
ylabel('y');

%Hiển thị lưới
grid on;

%Vẽ 2 trục tọa độ
xline(0);
yline(0);

%----------------------------Hết câu a------------------------------------%

%Nhập thời điểm cần xét cho các câu b, c, d
time = input('Nhap thoi diem can xet, t = ');

%----------------Hết phần nhập dữ liệu cho câu b, c, d--------------------%

%Câu b: Tính độ lớn vận tốc tại t = 1s
disp('Cau b: Do lon van toc chat diem tai t la:');
v = sqrt(subs(vx, t, time)^2 + subs(vy, t, time)^2);
disp(v);

%----------------------------Hết câu b------------------------------------%

%Câu c: Tính độ lớn gia tốc tại t = 1s
disp('Cau c: Do lon gia toc chat diem tai t la:');
a = sqrt(subs(ax, t, time)^2 + subs(ay, t, time)^2);
disp(a);

%----------------------------Hết câu c------------------------------------%

%Câu d:
%Tính bán kính cong của quỹ đạo tại t = 1s
disp('Cau d: Ban kinh cong quy dao tai t la:');

%Vận tốc tại thời điểm t
v = sqrt(vx^2+vy^2);

%Tính gia tốc tiếp tuyến
at = diff(v);

%Tính gia tốc hướng tâm
an = sqrt(a^2 - at^2);

%Tính bán kính quỹ đạo
R = subs(abs(v^2/an), t, time);
disp(R);
end