Eu tenho a seguinte função (curva de Viviani):

Phi     = @(t)[ cos(t)^2, cos(t)*sin(t), sin(t) ]

Apenas verifique se é válido:

s = linspace(0,T,1000);
plot3(cos(s).^2, cos(s).*sin(s), sin(s));

Como derivar a funçãoPhi (talvez várias vezes), o que representa a curva de Viviani em um pontot Ondet vai de0 para2*pi? Eu definiPhi adequado para tal derivado? eu tenteidiff, mas não manteve oPhi como eu precisaria disso.

Se a segunda derivada seriaPhi_d2, Preciso obter seu valor (por exemplo, emt = 0).

Como posso conseguir isso?