Based on the fun of:
Imagine we had a function of which its nth derivative at $x=n$ was some function:
$$\frac{d^{x}y(x)}{dx^x}=y^{(x)}(x)=ax+b\implies y(0)=b\implies y’(1)=a+b,\int y(x)dx\bigg|_{x=-1}=a-b$$
Now let’s use the fractional integral as the inverse of the xth derivative as $x$ is any constant. Therefore we can use Cauchy’s Repeated Integral formula:
$$y^{(x)}(x)=ax+b\implies y(x)+c_1=\text I^{x}_x (ax+b)$$
We can also generalize the lower integration bound as there may be no analog for the indefinite integral:
$$\,_a\text I^n_x g(x)\mathop=^\text{def}\frac 1{Γ(n)}\int_a^xg(t)(x-t)^{n-1}dt \implies y(x)+c_1=\,_{c_0}\text I^x_x (ax+b)= \frac 1{Γ(x)}\int_{c_0}^x(at+b)(x-t)^{x-1}dt \mathop=^{x\ne 0,-1}\\\implies y(x)=^?\frac{(x-c_0)^x(ac_0x+ax+bx+b)}{(x+1)!}+c_1$$
Here is an interactive graph of the conjectured answer.
The problem is that this formula only works for the trivial $a=b=0$ case.
Here is my graphical attempt for $y^{(x)}(x)=1$: $$y(0)=1,y’(1)=1,y’’(1)=1$$
maybe we can assume a function looking like the following with all constants of integration being $0$ for simplicity. Here is the graph with these properties:
Notice how this function matches up with our conjectured result:
$$y^{(x)}(x)=1\implies y(x)=\frac{x^x}{x!}$$
Actually we can generalize this process by doing:
$$y^{f(x)}(x)=g(x)\implies y(x)=\text I_x^n g(x),n=f(x)$$
given that $f\ne y$ so that we can solve the equation.
Then just use a definition of the $n$th integral or derivative and to finish.
I have also found this solution $n\in\Bbb N$ derivatives for $ax+b=1$ with many initial conditions:
$$y’’’(x)=1,y(0)=y’(1)=y’’(2)=y’’’(3)=1\implies y^{(n)}(x)=1,y^{(n)}(n)=y^{(n-1)}(n-1)=1=…=y’’(2)=1,y’(1)=1,y(0)=1:$$
which would give the following polynomial for $n=6$:
$$y(x)= \frac1{720} (720 - 4866 x + 4560 x^2 - 1500 x^3 + 270 x^4 - 24 x^5 + x^6), y^{(6)}(6)= y^{(5)}(5)= y^{(4)}(4)= y^{(3)}(3)= y^{(2)}(2) = y’(1)=y(0)=1 $$
Therefore our conjectured solution for $n\in\Bbb N$ is:
$$y=\frac1{n!}\sum_{k=0}^n c_k(-1)^k x^k,$$
where $c_k$ are the coefficients of our solution polynomial which may relate to these Polynomial function.
A self answer found the general solution assuming natural number order derivatives.
Therefore, how can $y^{(x)}(x)=ax+b$ be solved for $y(x)$ using fractional derivatives? Please correct me and give me feedback!