作者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