-<<fig=TRUE,echo=FALSE,results=hide,width=7,height=5>>=
-curve(to.kcal(10^x),from=0,to=max(c(sd(abs(log(c(0.8,1.33)))),
- sd(abs(log(c(1,1.33)))),
- sd(abs(log(c(0.8,1)))))),
- main="Curvature forward",
- xlab="Standard Deviation of Absolute value of the Log of the Curvature of Vesicle",
- ylab="Curvature Forward Adjustment (kcal/mol)")
+<<fig=TRUE,echo=FALSE,results=hide,width=7,height=7>>=
+grid <- expand.grid(x=seq(0,max(c(sd(abs(log(c(1,3)))),
+ sd(abs(log(c(1,0.33)))),sd(abs(log(c(0.33,3)))))),length.out=20),
+ y=seq(0,max(c(mean(log(c(1,3)),
+ mean(log(c(1,0.33))),
+ mean(log(c(0.33,3)))))),length.out=20))
+grid$z <- to.kcal(10^(grid$x*grid$y))
+print(wireframe(z~x*y,grid,cuts=50,
+ drape=TRUE,
+ scales=list(arrows=FALSE),
+ xlab=list("Vesicle stdev log curvature",rot=30),
+ ylab=list("Vesicle average log curvature",rot=-35),
+ zlab=list("Vesicle Curvature Forward (kcal/mol)",rot=93)))
+rm(grid)