你知道r是什么精度吗

R默认的double精度正常情况下应该够用了，感觉是你算法有问题。

  myfun <- function(P, y){
    ropt <- P[1]
    topt <- P[2]
    tmin <- P[3]
    tmax <- P[4]
    fun0 <- function(x){
      x[x < tmin] <- tmin
      x[x > tmax] <- tmax
      return(x)
    }
    y <- fun0(y)
    temp1 <- (ropt*(-(y-tmin)^3*(2*tmax+tmin-3*topt)^3+3*(y-tmin)^2*(2*tmax+tmin-3*topt)^2*(tmax-topt)^2-
      6*(y-tmin)*(2*tmax+tmin-3*topt)*(tmax+tmin-2*topt)*(tmax-topt)^2*(tmin-topt)
      +6*(tmax+tmin-2*topt)^2*(tmax-topt)^2*(tmin-topt)^2*log(y*(2*tmax+tmin-3*topt)+2*topt^2
      -tmax*(tmin+topt))))/( 3*(2*tmax+tmin-3*topt)^4*(-tmin+topt) )
    z     <- tmin
    temp2 <- (ropt*(-(z-tmin)^3*(2*tmax+tmin-3*topt)^3+3*(z-tmin)^2*(2*tmax+tmin-3*topt)^2*(tmax-topt)^2-
      6*(z-tmin)*(2*tmax+tmin-3*topt)*(tmax+tmin-2*topt)*(tmax-topt)^2*(tmin-topt)
      +6*(tmax+tmin-2*topt)^2*(tmax-topt)^2*(tmin-topt)^2*log(z*(2*tmax+tmin-3*topt)+2*topt^2
      -tmax*(tmin+topt))))/( 3*(2*tmax+tmin-3*topt)^4*(-tmin+topt) )
    temp3 <- temp1-temp2
    return(temp3)

  }


x <- seq(0, 100, len=2000)

par <- c( 0.8004191, 59.2451616, 24.4516643, 76.6365741)

plot(x, myfun(par, x), type="l")


谢谢上述两位的讨论，我附上一个函数，这个函数有四个参数。方程式较为复杂，我猜想就是因为这个复杂方程影响了计算的精度，导致曲线不光滑，按照道理，上述曲线应该是光滑的。

谢谢你们继续讨论！

你们稍微更改一下par中的任何一个值，稍微变化一下，就会出现一个光滑函数，为什么呢？

请将函数myfun 最后的return(temp3),改为return(data.frame(t1=temp1,t2=temp2,t3=temp3))。。。瞧瞧t1,t2,t3各是啥数？？我的计算结果是：

x <- seq(0, 100, len=2000)

par <- c( 0.8004191, 59.2451616, 24.4516643, 76.6365741)
yda<-myfun(par,x)
head(yda)
            t1           t2                     t3
1 2.519266e+15 2.519266e+15  0
2 2.519266e+15 2.519266e+15  0
3 2.519266e+15 2.519266e+15  0
4 2.519266e+15 2.519266e+15  0
5 2.519266e+15 2.519266e+15  0
6 2.519266e+15 2.519266e+15  0
tail(yda)
               t1           t2                      t3
1995 2.519266e+15 2.519266e+15 23.5
1996 2.519266e+15 2.519266e+15 23.5
1997 2.519266e+15 2.519266e+15 23.5
1998 2.519266e+15 2.519266e+15 23.5
1999 2.519266e+15 2.519266e+15 23.5
2000 2.519266e+15 2.519266e+15 23.5
summary(yda)
       t1                  t2                               t3      
Min.   :2.519e+15   Min.   :2.519e+15   Min.   :-0.50  
1st Qu.:2.519e+15   1st Qu.:2.519e+15   1st Qu.: 0.00  
Median :2.519e+15   Median :2.519e+15   Median : 7.00  
Mean   :2.519e+15   Mean   :2.519e+15   Mean   :10.47  
3rd Qu.:2.519e+15   3rd Qu.:2.519e+15   3rd Qu.:23.50  
Max.   :2.519e+15   Max.   :2.519e+15   Max.   :24.00