FOURIER(y,t,t1,t2,n) PULSE(t):=STEP(t)-STEP(t-1) PULSE(t) PULSEF(t):=FOURIER(PULSE(t),t,0,4,64) PULSEF(t) FOURIERS(y,t,t1,t2,n):=APPEND([[0,1/(t2-t1)*INT(y,t,t1,t2)]],VECTOR([m_,1/(t2~ -t1)*SQRT(INT(y*COS(2*pi*m_*t/(t2-t1)),t,t1,t2)^2+INT(y*SIN(2*pi*m_*t/(t2-t1)~ ),t,t1,t2)^2)],m_,1,n)) pulsefs:=FOURIERS(PULSE(t),t,0,4,64) pulsefs PULSE1(t):=PULSE(t)*IF(t<1/10,SIN(pi*t*5),IF(t>9/10,SIN(pi*(1-t)*5),1)) PULSE1(t) PULSEF1(t):=APPROX(FOURIER(PULSE1(t),t,0,4,64),6) PULSEF1(t) pulsefs1:=APPROX(FOURIERS(PULSE1(t),t,0,4,64),6) pulsefs1 PULSE2(t):=PULSE(t)*SIN(pi*t) PULSE2(t) PULSEF2(t):=APPROX(FOURIER(PULSE2(t),t,0,4,64),6) PULSEF2(t) pulsefs2:=APPROX(FOURIERS(PULSE2(t),t,0,4,64),6) pulsefs2 PULSE3(t):=PULSE(t)*(1-COS(2*pi*t))/2 PULSE3(t) PULSEF3(t):=APPROX(FOURIER(PULSE3(t),t,0,4,64),6) PULSEF3(t) pulsefs3:=APPROX(FOURIERS(PULSE3(t),t,0,4,64),6) pulsefs3