x=c(0.4,0.6,0.8,1,1.2,1.4,1.6,1.8,2,2.2,2.4) y=c(0,0.15,0.4,0.6,0.65,0.75,0.95,1,1,1,1) model= lm(y~ploy(x,3)) plot(x,y) lines(x,fitted(model) )