作者blackalert (右上角星位)
看板Statistics
标题Re: [问题] r绘图求程式码
时间Tue Mar 12 00:31:02 2013
※ 引述《esther730 (esther)》之铭言:
: f<-function(x)x^4-4*x^3-10
: dx<-deriv(~x^4-4*x^3-10,"x");dx
: d2x<-deriv(~4*x^3-12*x^2,"x");d2x
: f(-1)
: 以上是我做到一半的
: 函数:x^4-4*x^3-10
: 要画函数图和一次微分图,二次微分图,通过(-1,-5)该点切线
: 要本四条线画在一起并用不同颜色标记
: 可是我怎麽试都错误
: Ex:plot(f,add=T)
: 设x,y范围也出现问题
: 求以上程式码
## x轴范围 ##
from = -2
to = 4
## 在范围内切200个点 ##
x = seq(from, to, len=200)
## 在这200个x下的f(x),f'(x),f"(x) ##
y0 = f(x)
y1 = attr(eval(dx), "gradient")
y2 = attr(eval(d2x), "gradient")
## 绘图,也可以用curve跟lines ##
y = cbind(y0, y1, y2)
matplot(x, y, type="l", lty=1)
## 切线,-21跟-16是截距和斜率 ##
abline(-21, -16, col=4)
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 1.168.92.170