Осваиваю матлаб

[pinochet]

Начинает нравиться за неожиданные, но красивые результаты -


А всего-то

>> [T,Y]=ode15s(@sgl,[0,6],[0,0,1]);
>> fre=fft(Y)
>> plot(fre)

где

function M = sgl(t,y)
M = [100*y(2)+2*y(3)*0*y(1)-y(1); -100*y(1)+2*y(3)*1*cos(100*t)-y(2); 1-2*0*y(1)*y(1)-2*1*y(2)*cos(100*t)-y(3)];