Risposto

Crank-Nicholson method

I suggest the following code: xstart = 0; xend = 1; tstart = 0; tend = 1; nx = 10; nt = 100; x = linspace(xstart,xend,nx)...

Crank-Nicholson method

I suggest the following code: xstart = 0; xend = 1; tstart = 0; tend = 1; nx = 10; nt = 100; x = linspace(xstart,xend,nx)...

19 giorni fa | 0

| accettato

Risposto

Calculate total heat loss by conduction given temperature and depth profile vectors

Here is another way to determine the temporal change of heat content in the substratum. rho*cp*dT/dt = d/dx (k*dT/dx) Integrat...

Calculate total heat loss by conduction given temperature and depth profile vectors

Here is another way to determine the temporal change of heat content in the substratum. rho*cp*dT/dt = d/dx (k*dT/dx) Integrat...

21 giorni fa | 0

| accettato

Risposto

Calculate total heat loss by conduction given temperature and depth profile vectors

How can I estimate the total heat loss (w/m2) at each time interval (i.e., the heat loss for every profile)? Total heat flow in...

Calculate total heat loss by conduction given temperature and depth profile vectors

How can I estimate the total heat loss (w/m2) at each time interval (i.e., the heat loss for every profile)? Total heat flow in...

21 giorni fa | 0

Risposto

data not generating and error in plotting

And what is the "correct" result ? According to your mathematical description, I get this: % Define initial parameters lambda...

data not generating and error in plotting

And what is the "correct" result ? According to your mathematical description, I get this: % Define initial parameters lambda...

22 giorni fa | 0

| accettato

Risposto

Problem with the fitting

% Define the data dr_data = [2.5453, 0.042123; 5.0907, 0.075326; 7.636, 0.059506; 10.1813, 0.071553; 12.7267, 0.071365; 15.272,...

Problem with the fitting

% Define the data dr_data = [2.5453, 0.042123; 5.0907, 0.075326; 7.636, 0.059506; 10.1813, 0.071553; 12.7267, 0.071365; 15.272,...

22 giorni fa | 0

Risposto

Second order coupled differential equation

X0 = [0; 0; 0]; % initial conditions r = linspace(1, 50, 100); % resistance range tspan = linspace(0,100, 100); freq_vector...

Second order coupled differential equation

X0 = [0; 0; 0]; % initial conditions r = linspace(1, 50, 100); % resistance range tspan = linspace(0,100, 100); freq_vector...

22 giorni fa | 0

| accettato

Risposto

Use of Fmincon for stress relax data

I suggest you read the examples for the different suitable solvers. They will show you exactly how to proceed for your case. Ma...

Use of Fmincon for stress relax data

I suggest you read the examples for the different suitable solvers. They will show you exactly how to proceed for your case. Ma...

22 giorni fa | 0

Risposto

How to use Matlab to generate a discrete uniformly distributed decimal with a value range of {0, 0.05} and a step size of 0.01?

n = 10; x = randi(6,n,1); x = (x-1)*0.01

How to use Matlab to generate a discrete uniformly distributed decimal with a value range of {0, 0.05} and a step size of 0.01?

n = 10; x = randi(6,n,1); x = (x-1)*0.01

23 giorni fa | 0

Risposto

use of pdepe for a space-dependent diffusivity

First: Don't name the vector "diff" since "diff" is an internal MATLAB function. Name it D, e.g. Second: To get the correct val...

use of pdepe for a space-dependent diffusivity

First: Don't name the vector "diff" since "diff" is an internal MATLAB function. Name it D, e.g. Second: To get the correct val...

23 giorni fa | 0

| accettato

Risposto

If I have an array and I need to find an element in the array how do I go about it?

Compute 100*states(:,1)+10*states(:,2)+states(:,3) and use "find" to compare with C.

If I have an array and I need to find an element in the array how do I go about it?

Compute 100*states(:,1)+10*states(:,2)+states(:,3) and use "find" to compare with C.

24 giorni fa | 1

Risposto

Solve a differential equation and determine when it gets close to steady state

You can use an event function for the ODE integrator and stop integration when abs(dy/dt) (in your case abs(0.1 * (c-20) * (23-c...

Solve a differential equation and determine when it gets close to steady state

You can use an event function for the ODE integrator and stop integration when abs(dy/dt) (in your case abs(0.1 * (c-20) * (23-c...

24 giorni fa | 0

Risposto

the degradation of ozone in aqueous solutions (water). When I plot the all three components, the plot stays at the 0 value for all components.

There seems to be a singularity in the solution at t approximately 320. ProjectOzone() function ProjectOzone clc ki = 7e-3;...

the degradation of ozone in aqueous solutions (water). When I plot the all three components, the plot stays at the 0 value for all components.

There seems to be a singularity in the solution at t approximately 320. ProjectOzone() function ProjectOzone clc ki = 7e-3;...

25 giorni fa | 0

Risposto

The code below runs nonstop, and the results are not shown. why is that?

The code below runs nonstop, and the results are not shown. why is that? "solve" won't succeed for your system of equations bec...

The code below runs nonstop, and the results are not shown. why is that?

The code below runs nonstop, and the results are not shown. why is that? "solve" won't succeed for your system of equations bec...

25 giorni fa | 0

Risposto

Create number of for loops depending on size of N

The N columns of the resulting C-matrix contain i_first, i_second,... imax = 8; N = 4; i_average = (sum(0:N-1)/N+sum(imax:-1:...

Create number of for loops depending on size of N

The N columns of the resulting C-matrix contain i_first, i_second,... imax = 8; N = 4; i_average = (sum(0:N-1)/N+sum(imax:-1:...

26 giorni fa | 0

Risposto

How to make Matlab recognize that an expression rewritten by using subexpr() is equal to the original expression?

You have to resubstitute sigma by (a+b+c+d+e)^(1/2) in x_2 to compare with x: syms a b c d e f x x = a + (a+b+c+d+e)^(1/2) + c...

How to make Matlab recognize that an expression rewritten by using subexpr() is equal to the original expression?

You have to resubstitute sigma by (a+b+c+d+e)^(1/2) in x_2 to compare with x: syms a b c d e f x x = a + (a+b+c+d+e)^(1/2) + c...

26 giorni fa | 0

Risposto

plot the output c(t) using mathlab and show setting time on you graphf

Use "fplot" as done in your code above.

plot the output c(t) using mathlab and show setting time on you graphf

Use "fplot" as done in your code above.

26 giorni fa | 0

Risposto

Describe the solution trajectory of differential systems with mixed time-varying delay

To write your equations as a system of delay differential equations that can be handled by a MATLAB solver, set x3 = integral_{...

Describe the solution trajectory of differential systems with mixed time-varying delay

To write your equations as a system of delay differential equations that can be handled by a MATLAB solver, set x3 = integral_{...

26 giorni fa | 1

| accettato

Risposto

Plot between Nu vs Ra

Replace Ra_values = 10:100:100000; by Ra_values = (10:100:100000).';

Plot between Nu vs Ra

Replace Ra_values = 10:100:100000; by Ra_values = (10:100:100000).';

26 giorni fa | 0

Risposto

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...

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...

27 giorni fa | 1

Risposto

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...

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...

27 giorni fa | 0

Risposto

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...

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...

27 giorni fa | 0

| accettato

Risposto

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 ...

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 ...

27 giorni fa | 0

Risposto

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...

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...

28 giorni fa | 0

Risposto

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...

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...

29 giorni fa | 0

| accettato

Risposto

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 = ...

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 = ...

29 giorni fa | 1

| accettato

Risposto

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...

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...

29 giorni fa | 0

Risposto

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 ...

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 ...

30 giorni fa | 1

Risposto

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 + ...

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 + ...

30 giorni fa | 1

| accettato

Risposto

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) -...

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) -...

circa un mese fa | 1

| accettato

Risposto

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...

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...

circa un mese fa | 1