第一個(gè)問(wèn)題：
f[x_] := 2 x^2;
xlist = NestList[# - f[#]/f'[#] &,0.5,4]; ylist = f[xlist]; list1 = Transpose[{xlist,ylist}]; list2 =
Transpose[{Drop[xlist,1],Table[0,{Length[xlist] - 1}]}];
list = Riffle[list1,list2];
iterplot = ListLinePlot[list]; funplot = Plot[f[x],{x,0,1},PlotRange -> {0,1},PlotStyle -> {Black}]; Show[{funplot,iterplot}]
第二個(gè)問(wèn)題：
singleTest := Module[{data,dataSorted,exist,howmanysame = 2},
data = RandomInteger[{1,365},90]; dataSorted = Sort[data];
If[Or @@ (Map[dataSorted[[#]] == dataSorted[[# + howmanysame]] &,
Range[Length[data] - howmanysame]]),exist = 1,exist = 0];
exist];
result = Table[singleTest,{200}];
In[203]:= Total[result]/200 // N
Out[203]= 0.505
第三個(gè)問(wèn)題：
x[t_] := t; y[t_] := Sin[t]; r = 0.5;
Manipulate[
Show[
{ParametricPlot[{x[t],y[t]},{t,0,2 \[Pi]},
PlotRange -> {{0,2 \[Pi]},{-2,2}}],
Graphics[{
Circle[{x[a],y[a]} + r*{- y'[a],x'[a]}/Norm[{x'[a],y'[a]}],r],
Point[{x[a],y[a]} + r*{- y'[a],x'[a]}/Norm[{x'[a],y'[a]}]]}]}
],
{a,0,2 \[Pi]}]
第二個(gè)和第三個(gè)粘貼反了.一些意見(jiàn),1里面方程似乎沒(méi)有根；3里面到底是三個(gè)還是三個(gè)以上說(shuō)的不清楚.這個(gè)讓人很不爽.
2里面我只做了圓心,算是hint吧,其他可以自己搞定的
那么1,3解決了
關(guān)于2你說(shuō)的無(wú)摩擦滾動(dòng)=滑動(dòng)+繞圓心一定速率的轉(zhuǎn)動(dòng).依我看只有用弧長(zhǎng)參數(shù)代碼才容易寫(xiě).不過(guò)就是x,y各增加一項(xiàng),琢磨琢磨吧:)