\[\left(\frac{e^{-uk(t-D)}}{1+e^{-uk(t-D)}}\right) \left(\frac{1-e^{-dk(t-L-D)}}{1+e^{-dk(t-L-D)}}\right)^{pw}\]
c1 <- makepar_F_sharkfin()
C1<- make_function(c1)
tt <- seq(0, 565, by=5)
plot(tt, C1(tt), type ="l", xlab = "Time (in Days)", ylab = expression(Fs1(t)))
c2a <- makepar_F_sharkfin(L = 90, dk = 1/110)
c2b <- makepar_F_sharkfin(L = 180, dk = 1/40)
c2 <- makepar_F_sharkfin(L = c(90, 180), dk = c(1/110, 1/40), pw=c(2,1), N=2)
C2<- make_function(c2)
C2t <- C2(tt)
plot(tt, C2t[1,], type ="l", xlab = "Time (in Days)", ylab = expression(Fs1(t)))
lines(tt, C2t[2,])