参考: Interface(インターフェース) 2017年10月号, pp.158-159
initRad=0.; deltaRad=Pi/8.; sinDeltaRad=Sin[deltaRad]; cosDeltaRad=Cos[deltaRad]; sol=RecurrenceTable[{rad[n]==rad[n-1]+deltaRad, sin[n]==sin[n-1]*cosDeltaRad+cos[n-1]*sinDeltaRad, cos[n]==cos[n-1]*cosDeltaRad-sin[n-1]*sinDeltaRad, rad[0]==initRad, sin[0]==Sin[initRad], cos[0]==Cos[initRad]}, {rad,sin,cos}, {n,0,50}]; radList=sol[[;;,1]]; sinList=sol[[;;,2]]; cosList=sol[[;;,3]]; sinPlot=ListLinePlot[Transpose[{radList,sinList}],PlotStyle->Red]; cosPlot=ListLinePlot[Transpose[{radList,cosList}],PlotStyle->Blue]; Show[sinPlot,cosPlot] sol // Grid