在R中绘制一个水平集

Question

知道如何在R中绘制这个等式吗？我被要求这样做，不知道如何绘制这个没有给出数字。

p(x, y) = xy - max(.008x, .00075xy)

约束：

x = exp(-y)/(1+exp(-y))  

任何帮助都会很棒，谢谢！

Answer 1

你可以定义两个这样的函数：

#first you create the second equation i.e. x
f_x <- function(x) exp(-x)/(1+exp(-x))

#and then you rewrite the first function where you replace x with exp(-y)/(1+exp(-y))  
f_y <-  function(y) f_x(y) * y - max(.008 * f_x(y), .00075 * y * f_x(y))

然后使用curve：

绘制它

curve(f_y)

输出：