If you use the editor, you can set a stop as normal, right-click on it, select “set/modify condition” and enter the condition (the stop will turn from red to yellow). From command line, you can use dbstop in file if expression dbstop in file at location if expression e.g. dbstop in myFile at 200 if … Read more
You want make F an argument of your derivative function and pass the right anonymous function to ode45: [t,N] = ode45(@(t,y) rateEquations(t,y,F), timeSpan, initialConditions)
Introduction & Test Setup To perform memory tests to inquire about the points raised in the question, let’s define the inputs A and B: A = rand(M,N) B = rand(1,N) Here, M and N are the size parameters and are kept as really large numbers. I would be using repmat for comparisons as that seems … Read more
From the documentation on parfor: The body of a parfor-loop cannot contain global or persistent variable declarations. In the context of your problem, i.e., calling a function within the parfor that in turn references a global, this translates into: “parfor will probably not give expected or meaningful results”. This makes perfect sense. Consider the following … Read more
Fortran doesn’t support anonymous functions. The simple work around is to write a function that has a name. There are then two possible approaches in modern Fortran for capturing the value of any additional parameters required for the function beyond the variable being minimised: The procedure to be minimised is expressed as a deferred binding … Read more
Use bwconncomp instead since it returns the coordinate indexes for region in a separate cell, where the size of each is easily discerned: >> BW = [1 0 0; 0 0 0; 0 1 1]; % two regions >> CC = bwconncomp(BW) CC = Connectivity: 8 ImageSize: [3 3] NumObjects: 2 PixelIdxList: {[1] [2×1 double]} … Read more
You can take a look at the code here. I think the function is similar to what you are trying to solve. The function extracts an arbitrary plane from a volume given the size of the plane, the center point of the plane, and the plane normal, i.e. [A,B,C]. It also outputs the volumetric index … Read more
It seems like it might be the best to use datetick. Usage: datetick(‘x’) or datetick(‘x’, dateformat) with the available formats as explained in the documentation.
Here is an example using scatteredInterpolant: %# get some 2D matrix, and plot as surface A = peaks(15); subplot(121), surf(A) %# create interpolant [X,Y] = meshgrid(1:size(A,2), 1:size(A,1)); F = scatteredInterpolant(X(:), Y(:), A(:), ‘linear’); %# interpolate over a finer grid [U,V] = meshgrid(linspace(1,size(A,2),50), linspace(1,size(A,1),50)); subplot(122), surf(U,V, F(U,V)) Note that you can evaluate the interpolant object at … Read more
To round to d significant digits: >> d = 3; %// number of digits >> x = 5.237234; %// example number >> D = 10^(d-ceil(log10(x))); >> y = round(x*D)/D y = 5.2400 To round to d decimal digits: >> d = 3; %// number of digits >> x = 5.237234; %// example number >> D … Read more