Rstudioのmanipulate関数でカスプ・カタストロフのポテンシャル関数を描く

カスプ・カタストロフのポテンシャル関数

(1/4)*x^4+(1/2)*u*x^2+v*x

における分裂要因（u）と平常要因（v）をコントロールして最適解x（ここではポテンシャルを局所的に最小化するx。x^4の係数を負にすれば最大化のモデルになる）がどのように変化するかを示す。

library(manipulate)
manipulate(
curve*1,
u = slider(-15,15),
v = slider(-15,15)
)

※7/15 expression関数によるタイトルを追加するなど修正

*1:1/4)*x^4+(1/2)*u*x^2+v*x,-6,6, cex.lab = 1.4, ylab = "F(x)",
main = expression(F(x) == frac(1,4)*x^4 + frac(1,2)*ux^2 + vx