fplot | Nice distribution plots with minimum user input

First their are few things to note here :

In the second part, replace your line

The bounds/range of the plot can be changed by adjusting X_Min, X_Max, Y_Min and Y_Max as you'd like.

There are several problems with the code:

• nchoosek does not vectorize on the second input, that is, it does not accept an array as input. fplot works faster for vectorized functions. Otherwise it can be used, but it issues a warning.
• The result of nchoosek is close to overflowing for such large values of the first input. For example, nchoosek(1000,500) gives 2.702882409454366e+299, and issues a warning.
• nchoosek expects integer inputs. fplot uses in general non-integer values within the specified limits, and so nchoosek issues an error.

You can solve these three issues exploiting the relationship between the factorial and the gamma function and the fact that Matlab has gammaln, which directly computes the logarithm of the gamma function:

Your syntax for calling fplot is the problem, not your function. If you're passing a simple function handle, just use:

Just in case anyone else comes across this problem, the solution is the bias term can be used to find the y-intercept, b = -SVMModel.Bias/betas(2). And the y-intercept is just another point in space [0 b] which can be recovered/unrotated by inverse transforming it through the PCA. This new point can then be used to solve the linear equation y = mx + b (i.e., b = y - mx). So the code should be:

fplot requires an expression that returns a numeric value. f returns a sym (symbolic expression), so it doesn't work, whereas z returns a number - so it does. That's all the difference.

Note that in the working example, you overwrote the y sym on the 2^nd row, which means you provided a "proper" function handle to fplot, "by mistake". In fact, you need much less code for that example to work:

Yes, they are different. For the first code:

I think you just need to swap these steps: