I will be highly grateful for your help. I am learning how to use *geom\_function* in r. Following is my function:

```
x0 <- 0.5
x1 <- 1
x2 <- 2
x3 <- 3

x <- c(x0, x1, x2, x3)

myfn <- function(w, b, a, x){
    w^(1-b)/(1-b)-a*w-(w>100 & w<=200)*x[3]*(w-100)-(w>200)*x[3]*100-x[4]*(w-200)
}

```

My objective is to plot above function using *geom\_function* to see how this function behaves with different values of arguments *a* and *b* and following is my code:

```
y=seq(0,1000,5)

ggplot()+
  xlim(c(0,1000))+
  geom_function(fun=myfn(w=y, b=-4, a=0.5, x=x))

```

Problem: I feel my logic is correct but when I execute above code, I get nothing. I will be highly grateful for the help. Thank you very much for the help in advance. Any help or direction will be highly appreciated.