Answered
Plot between Nu vs Ra
Replace Ra_values = 10:100:100000; by Ra_values = (10:100:100000).';

1 year ago | 0

Answered
Resetting only a section of my state during event driven ode45 simulation
After the event has been triggered, control is returned to the calling program. If you want to continue integration with the e...

1 year ago | 1

Answered
Index exceeds the number of array elements. Index must not exceed 1.
path_length_BS_IRS_reflect_UE(i) does not exist for i > 1 because path_length_BS_IRS_reflect_UE is a scalar value, namely 223.60...

1 year ago | 0

Answered
Dot indexing is not supported for variables of this type.
Test before your loop % Loop through each entry in the list and extract data for i = 1:numEntries dt(i) = dateti...

1 year ago | 0

| accepted

Answered
How do I know what method Matlab is using when plotting ltiview?
I don't know Itiview nor risetime, but isn't there a literature link in the documentation for the function you are using ? If ...

1 year ago | 0

Answered
How to solve algebraic equations for different values of variables ?
% Solve the equation eqn = Ac.*I.*eta-m.*c.*(Te-Ti)==0 for i=1:numel(eqn) s = solve(eqn(i),Te); S(i) = vpa(s(s>=0)); en...

2 years ago | 0

Answered
function file used with 2 Matrices as input
You are not allowed to use the same name for the function and an output array. Thus rename one of them. And ln_AOU as output pa...

2 years ago | 0

| accepted

Answered
Solving a linear programming problem uising linprog
A = -[110 205 160 160 420 260;... 4 32 13 8 4 14;... 2 12 54 285 22 80]; b = -[2000; 55; 800]; Aeq = []; beq = ...

2 years ago | 1

| accepted

Answered
BVP Error, help with explaination
The solver is no longer able to capture the steep gradient of your solution at x=0. Look at sol.x in the last step of the cont...

2 years ago | 0

Answered
Data type symfun (code problem)
If your MATLAB release throws an error (R2024a does not as you can see above), try [n(x) d(x)] = numden(df(x)) instead of [n ...

2 years ago | 1

Answered
The estimation error is strangely obtained from Simpson's 1/3 rule...
syms a b c d e real syms t real f(t) = a*t^5 + b*t^4 + c*t^3 + d*t^2 + e*t; s = 0; i = -1/2; while i < 959.5 i = i + ...

2 years ago | 2

| accepted

Answered
V2H Optimization: No feasible solution found
Use energyStorage = optimconstr(T); energyStorage(1) = BSS_SOC(1) == 0; energyStorage(2:T) = BSS_SOC(2:T) == BSS_SOC(1:T-1) -...

2 years ago | 1

| accepted

Answered
Need to solve the following equation with three knowns and 2 unknowns.
H_BM = rand(100,1); g_mean = rand(100,1); Wi = rand(100,1); Hi = H_BM; gi = g_mean; x = Hi.*gi; y = Wi + Hi.*gi; fitl...

2 years ago | 1

Answered
Why does Matlab not recognize fieldnames function?
You overwrite the "fieldnames" function by setting it to a variable value: fieldnames = fieldnames(results.digitalMap);

2 years ago | 0

| accepted

Answered
Simulate Fick's 1st and 2nd laws of mass diffusion.
The correct update for inner grid points is C_new(i,j) = C_old(i,j) + dt * (D_x*(C_old(i-1,j)-2*C_old(i,j)+C_old(i+1,j))/dx^2 +...

2 years ago | 0

| accepted

Answered
Indices shown in plot are off
My guess is that the x-axis has a resolution of 10.

2 years ago | 0

Answered
issue in calculation of wave speed
The below code will at least give you the correct update of the solution. You don't initialize "half_max_positions" ; that's wh...

2 years ago | 0

Answered
Minimize difference between two data sets
Use [output] = [method](@(AB) my_func(AB(1),AB(2)), [A B]) instead of [output] = [method](@(A,B) my_func(A,B), [A B])

2 years ago | 0

| accepted

Answered
1D advection diffusion with pdepe
You mean something like this ? X = linspace(0, 1, 100); W = linspace(0, 2, 30); m = 0; sol = pdepe(m, @pdefun, @ic, @bc, X, ...

2 years ago | 0

| accepted

Answered
traveling waves of fisher kpp equation
Your equation reads u_{t} = - (u_{xx}+u_{yy}) + u(1-u) You miss the minus sign in your code (at least if the equation you wrot...

2 years ago | 0

Answered
Summing results at the end of a loop
x = linspace(0, 90.*(10.^-3), 100); w = [7.5.*(10.^-3),7.5.*(10.^-3),7.5.*(10.^-3)]; u = [(8.883).*(10.^-3),(8.883).*(10.^-3),...

2 years ago | 0

Answered
calculate the normal of a 3D plane
I would like to know if it is correct to calculate the normal of this plane in the following way: N = cross(P(1,:) - P(2,:), P(...

2 years ago | 0

| accepted

Answered
"Empty sym" when solving a system of linear equations with different conditions
Your system only has a (nontrivial) solution if i = 1. In this case, the solution is given as below. syms Z1A Z2A Z3A Z4A Z1B Z...

2 years ago | 0

| accepted

Answered
Temperature profile in a pipe flow with Matlab pdepe solver
I corrected as much as I could, but the computation of Q looks different in your code in comparison to the article. %% Main cod...

2 years ago | 0

Answered
Substitute partial differential into symfun
syms G(h) h(t) t Q = G(h) dGdt = diff(Q,t) % -> D(G)(h(t)) occurs s = children(dGdt) % Now define G(h) and derive: G ...

2 years ago | 1

| accepted

Answered
Using a loop to add rows to a struct, but the counter overwrites the rows.
I don't know how "trial" and "blockInfo" are organized. Here is one suggestion: for ifile = 1:3 for s = 1:nTrials struc...

2 years ago | 0

Answered
I am trying to perform a simple operation FOR loop. I think I am not able to get the indexing right.
%% why does not the following operation work Because the numbers 1,2 and 3 in t1, t2 and t3 are not indices you can loop over, ...

2 years ago | 0

| accepted

Answered
Kelvin functions ker kei functions evaluation
The value kei(0) = -pi/4 only exists in the limit : syms x a = exp(pi*1i/4); limit(imag(besselk(0,a*x)),x,0,"right") limit(i...

2 years ago | 0

Answered
constrained optimization with Matlab, symbolic function
If your function f is symbolic, you first have to convert it to a numerical function by using "matlabFunction". k_tot = ...; a...

2 years ago | 0

Answered
I am unable to find the minimum value of function f in this optimization problem. The solver always returns 0 for all variables which is not possible. I am unable to correct.
The solution is correct: A*x <= b is satisfied, and since all elements in the f-vector are positive, 0 is optimal. Maybe you wa...

2 years ago | 0