Friday, December 27, 2019

Evolutes of Curves


F[t_]:={{2t,t^2},{7Sin[t],4Cos[t]+6},
{2t*Cos[t]^2,t^2*Cos[t]^2}};
DF[t_]:=D[F[t],t]; D2F[t_]:=D[DF[t],t];
NDF[t_]:=Table[Sqrt[DF[t][[i,1]]^2+DF[t][[i,2]]^2],{i,3}];
ND2F[t_]:=Table[Sqrt[D2F[t][[i,1]]^2+D2F[t][[i,2]]^2],{i,3}];
FNDF[t_]:=Table[{DF[t][[i,2]]/NDF[t][[i]],
-DF[t][[i,1]]/NDF[t][[i]]},{i,3}] ;
RF[t_]:=Table[NDF[t][[i]]^3/(DF[t][[i,
2]]*D2F[t][[i,1]]-DF[t][[i,1]]*D2F[t][[i,2]]),{i,3}] ;
CF[t_]:=Table[{F[t][[i,1]]+RF[t][[i]]*FNDF[t][[i,1]],
F[t][[i,2]]+RF[t][[i]]*FNDF[t][[i,2]]},{i,3}];
TF=Table[F[t],{t,-8,8,.1}];
TCF=Table[Evaluate[CF[t]],{t,-8,8,.1}];
Table[Show[ParametricPlot[{F[t][[i]],Evaluate[CF[t][[i]]]},
{t,-8,8},PlotStyle->Thick],
Table[ListLinePlot[Style[{TF[[k,i]],TCF[[k,i]]},
LightGray,Dashed]],{k,Length[TF[[All,i]]]}],
PlotRange->{{-8,8},{-2,14}},ImageSize->400],{i,3}]

No comments:

Post a Comment