我有以下功能（Viviani 曲线）：
Phi = @(t)[ cos(t)^2, cos(t)*sin(t), sin(t) ]
只需检查它是否有效：
s = linspace(0,T,1000);
plot3(cos(s).^2, cos(s).*sin(s), sin(s));
如何衍生功能Phi（可能多次），它代表了一个点维维安尼的曲线，t其中，t从去0到2*pi？我是否定义了Phi适合这样的衍生物？我试过了diff，但它没有Phi像我需要的那样保留它。
如果二阶导数是Phi_d2，我需要得到它的值（例如在 中t = 0）。
我怎样才能做到这一点？