Mam następującą funkcję (krzywa Viviani):

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

Po prostu sprawdź, czy jest ważny:

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

Jak wyprowadzić funkcjęPhi (może wiele razy), co reprezentuje krzywą Viviani w punkciet gdziet idzie od0 do2*pi? Czy zdefiniowałemPhi odpowiedni dla takiej pochodnej? próbowałemdiff, ale to nie zachowałoPhi jakbym tego potrzebował.

Gdyby druga pochodna byłaPhi_d2, Muszę uzyskać jego wartość (na przykład wt = 0).

Jak mogę to osiągnąć?