Answered
How to optimize the run time in my optimization problem.
I don't see anything obviously wrong in your coding - except for the last term in MTE_i which should be (TotalProduct(6)-Exper...

2 years ago | 0

| accepted

Answered
I'm getting the 'Not enough input arguments' error event though I am pretty sure I am giving all the arguments needed at all times
format long [x,it] = NewtonSystemSolver(5, 10^-12, 100, [5e4,0,0,0,5e4].', @Slammer1, @Slammer1Jacobian) Slammer1(x) function...

2 years ago | 1

| accepted

Answered
ode15s giving integration tolerance error and loop stopping at zeroth iteration
The evaluation of the variable "ros" for your vector of initial conditions y0est gives NaN. The reason is that you start with mL...

2 years ago | 0

Answered
Multiple iterations for plotting
Maybe you mean h= 0.5; b = 0.5; d1= 0.474; d2= 0.035; As= 2512*10^-6; Asc= 2512*10^-6; fy= 400*10^6; Es= 200*10^9; fc...

2 years ago | 0

| accepted

Answered
There is a matrix operation, A*B=c. would there be a way to find a specific elemetn in C and call out variables that were involved in the calculation of the specific element
Since you set C = A.*B and not C = A*B, the elements that influence the value at position (i,j) in C are the elements (i,j) in A...

2 years ago | 0

Answered
Set up a gradient descent algorithm for a multivariable economic dispatch problem
There must be something wrong in the computation of t in your code. For t=1, the usual gradient descent works: clc; clear; sym...

2 years ago | 0

| accepted

Answered
Setting Boundary Conditions & Using bvp4c
xcastle = 20; % x-coordinate of target v0 = 20; % Prescribed (...

2 years ago | 0

| accepted

Answered
Solving ODE system with a constraint/boundary condition
p_A + p_B should remain constant since dp_A/dz + dp_B/dz = 0 according to your equations. But this is not the case for p_A + p_...

2 years ago | 0

Answered
How do I write time dependent variables in my equations of motion?
Integrate from 0 to 5 with F = 0, Vin = 280. Restart with the solution obtained and integrate from 5 to 10 with F = -10, Vin = ...

2 years ago | 0

| accepted

Answered
What is the mathematical form called for PDEs used with the MATLAB function pdepe?
Would this be considered an elliptical/parabolic PDE in conservative form? No. "Conservative" usually refers to the first spati...

2 years ago | 1

| accepted

Answered
Energy optimization robot paths
toolPositionHome = T_home(1:3,4); % Define a Set of Waypoints Based on the Tool Position waypoints = toolPositionHome + ... ...

2 years ago | 0

Answered
I'm having an issue defining the initial conditions for solving a system of differential equations
You have two second-order ODEs for th and w. Thus you need 4 boundary conditions (not 6) that can be prescribed for w, dw/dx, ...

2 years ago | 1

| accepted

Answered
Received this error, not sure how to fix
Maybe you mean for ni=1:j instead of for ni=1+j If not, you try to access N(7) which does not exist.

2 years ago | 0

| accepted

Answered
Vector Length Error when Calling Function
You overwrite MATLAB's intrinsic "length" in the lines length = data(2:end,2) length = cell2mat(length)

2 years ago | 0

Answered
How to match matrix elements for a condition?
If A becomes larger, this brute-force way of solving will become intractable. A = [15 25 36 17 48 59 31 64 18 21 97 84 31 64 15...

2 years ago | 0

| accepted

Answered
How do I use my data set in an equation?
Say lm is 2x6, e.g. lb = [1 2 3 4 5 6] lm = [7 8 9 10 11 12;13 14 15 16 17 18] What would you like to compute ? [log(1/7) lo...

2 years ago | 0

Answered
error >> Matrix dimensions must agree. How to fix
z doesn't need to be prescribed in array form because you integrate with respect to z. Here I assume that the integration variab...

2 years ago | 0

Answered
How to get the formal expression of the expansion of power of a polynomial with symbolic calculation
syms x n = 10; k = 2; a = str2sym(compose("a" + (0:n-1))) f = x^n + sum(a.*x.^(0:n-1)) g = f^k; coefficients = coeffs(coll...

2 years ago | 0

Answered
Using MATLAB to do Numeric Calculus
Your code is correct. Here is a simpler way to define s and v. Just substitute tcrash by tmax. tcrash = 1596.2/(2*4.905)+sqrt...

2 years ago | 0

| accepted

Answered
probability that values are belong to an interval
It depends on how x is distributed. If f is the (continuous) probability density function for x, you get the probability p as ...

2 years ago | 0

| accepted

Answered
How to use ode15s to solve a stiff ode with mass?
tf = 200; x0 = [0 0 0 0 2500 0 0.8 0 50 50 50 50 5]; M = zeros(13,13); M(1,1) = 1; M(2,2) = 1; M(3,3) = 1; M(4,4) = 1; ...

2 years ago | 1

Answered
I am getting error:Warning: Failure at t=0.000000e+00. Unable to meet integration tolerances without reducing the step size below the smallest value allowed (7.905050e-323)
If you remove the semicolon behind %% mass balances dy_dz=[1/u*(-r3) 1/u*(...

2 years ago | 0

| accepted

Answered
vectorized operations on symbolic functions
rng("default") syms x1 x2 x3; % symbolic variables y = x1^3/3 + x2^2/2 - x3; % symbolic function y X = rand(500,3) % each row...

2 years ago | 0

Answered
Dot product does not work for symbolic vectors
Since you cannot assume functions to be real-valued, you must use dot(r,r) - (r_x*r_x' + r_y*r_y' + r_z*r_z' ) instead of dot...

2 years ago | 1

Answered
NaN in converting a sym to double
exp(lambda1*c0*eta0) and exp(lambda3*c0*eta0) evaluate to Infinity for s big enough (s slightly bigger than 0 already suffice...

2 years ago | 0

| accepted

Answered
How to plot a three variable complex function
You cannot plot the complete result of a function that depends on three independent variables. But you can plot it on slices (i...

2 years ago | 0

| accepted

Answered
PDEPE solver failure (spatial discretization) ; partial differential/derivative temperature modelling
Please include the constants defined as globals so that we can test your code. One obvious error is the boundary condition for...

2 years ago | 0

| accepted

Answered
Packed bed storage modeling probelm
rhoS = 2750; % Solid density [kg/m^3] cpS = 850; % Solid specific heat capacity [J/(kg*K)] ks = 1.28; ...

2 years ago | 0

| accepted

Answered
Extracting values from optimization variable and setting constraints
Defining vectors QTYup and QTYdown as integer vectors with 0<=QTYup(i)<=5 and -5<=QTYdown(i)<=0 (i=1,...,3 for tomatoes, onions ...

2 years ago | 0

Answered
I am trying to form a conditional loop that shows the growth of money, I put in $1000 every year and it grows by 85 for ten years but keep returning one value
money = zeros(10,1); supersaver = 1000; Return = 0.08; money(1) = supersaver*(1+Return); % Money at end of year 1 for i = 2...

2 years ago | 0

| accepted