All pastes #17456 Raw Edit

test

public text v1 · immutable
#17456 ·published 2005-07-11 15:44 UTC
rendered paste body
x=1:20000
x=x-10000
x=x*360/10000
 
sinusw1=function(x){
znak=1
x=x*2*pi/360
sinus=znak*x}
sinusw3=function(x){
znak=1
x=x*2*pi/360
sinus=znak*(x-x^3/factorial(3))}
sinusw5=function(x){
znak=1
x=x*2*pi/360
sinus=znak*(x-x^3/factorial(3)+x^5/factorial(5))}
sinusw7=function(x){
znak=1
x=x*2*pi/360
sinus=znak*(x-x^3/factorial(3)+x^5/factorial(5)-x^7/factorial(7))}
sinusw9=function(x){
znak=1
x=x*2*pi/360
sinus=znak*(x-x^3/factorial(3)+x^5/factorial(5)-x^7/factorial(7)+x^9/factorial(9))}
sinw=function(x) sin(2*pi*x/360)
sinww=sinw(x)
sinusw1w=sinusw1(x)
sinusw3w=sinusw3(x)
sinusw5w=sinusw5(x)
sinusw7w=sinusw7(x)
sinusw9w=sinusw9(x)
plot(x,sinww,col="red",lty=1,type="l",ylim=c(-2,2))
abline(h=c(-1,0,1))
lines(x,sinusw1w,col="blue",lty=3)
lines(x,sinusw3w,col="blue",lty=3)
lines(x,sinusw5w,col="blue",lty=3)
lines(x,sinusw7w,col="blue",lty=3)
lines(x,sinusw9w,col="blue",lty=3)
text(c(190,290,200,40,150),c(0.7,0.7,-0.9,1,-0.9),list("5","9","7","1","3"),col="blue")
text(340,-0.75,"sinus",col="red")