How can I skip the rest of an if statement?
Also in for-loops, "break" only exits the inner structure, but the outer loop runs until the end (see below). Thus you cannot ...

9 days ago | 0

Surface plot doesn't correspond to optimization solution
The problem is now, when I select a point randomly in the surface plot, and then calibrate the model for these values, I do not ...

9 days ago | 0

| accepted

Solving an equation with symbolic matrix

10 days ago | 1

A group of 12 people all roll a pair of dice 12 times, how often does each person roll a total number that is greater than or equal to 7?
In my opinion, it should be roll1=randi(6,12); roll2=randi(6,12); total=roll1+roll2; GT7 = total>=7; totalGT7 = sum(GT7,1);...

10 days ago | 1

| accepted

error in for loop
u(j-3) is only defined for j>=4, u(j+3) is only defined for j<=length(R)-3. Similar restrictions hold for the other indices. Thu...

10 days ago | 0

I'm suppose to plot the response of the system but i'm getting wrong values
alpha is a constant , and this is the value of it , I doubt that the rate of change at t=0 for x(5) should be 2.55e17 as you g...

11 days ago | 0

How can I multiply a row vector with each of two column vectors?
t = [1 2]; s = [(1:3)',(4:6)']; kron(t,s)

11 days ago | 3

final value of x in PDE
i want to solve this equation for r at given time and concentration I assume that the time for which you want to get r is in th...

12 days ago | 0

How to prevent substrate concentration going below minimum limit in fed batch fermentation model.
You can use the "Events" option of the ode integrators to interrupt and restart integration with different settings if a certain...

12 days ago | 0

| accepted

How can I ensure that the initial solution (x0) for fsolve does not result in Inf or NaN values?
You know your "external_fun" best. Try to deduce in advance which inputs will lead to failure. Maybe you can use "lsqnonlin" ins...

13 days ago | 0

Huge minimum value from expected results
I cannot run this code because it takes too long or maybe the matrix exponential cannot be computed analytically. Does it give ...

13 days ago | 0

| accepted

ODE45 solver returning constant value as the input variable.
As you can see when you run the below code, p is in the order of 1e12 which makes dL/dt effectivly 0. Check your units. %Def...

13 days ago | 1

| accepted

The code does not converge
I seems your function does not vary much with rho. The evaluation cannot find a value for rho for which the function has a zero....

13 days ago | 0

| accepted

How to plot the error of two numerical methods on the same graph?
%% Given data f=@(x) 8-4.5*(x-sin(x)); df=@(x) -4.5*(1-cos(x)); x0=1; tol=0.0001; n=50; [x_newton,i_newton,Error_newton]=n...

13 days ago | 1

Modifying the objective function in order to meet the optimization variable boundaries
if (A <= 1e-04) & (A >= 4e-04) How can A be <= 1e-4 and >= 4e-4 at the same time ? And if-statements have to be done elementwi...

14 days ago | 0

Solve a matrix equation
D = eig(A) will give you the values for D for which det(A-DI)=0

14 days ago | 0

| accepted

Solving goal programming problem with MATLAB
x0 must be a vector of size 1x5, not 1x4. Further, your function must return a vector of length 4, not 3. Further, the first o...

15 days ago | 0

plot() and fplot() give different results for the same function
x = linspace(0, 100, 10000); y = 2*sqrt(1+x)+sin(sqrt(pi)*x).*sin(sqrt(1+x)*pi).*(2*x+1)./sqrt(x) - 2*sqrt(1+x).*cos(sqrt(x)*pi...

15 days ago | 0

Vpasolve can not find a solution that I know it exist
Maybe it must read Eq2 = ((a21*W2)/(c2*(a21+a22))) - x12 == 0; instead of Eq2 = ((a12*W2)/(c2*(a21+a22))) - x12 == 0; ? c...

15 days ago | 0

| accepted

How to build a graph of the solution of a second-order differential equation in MATLAB
Since sin(pi) = 0, your differential equation reduces to y'' + y = 0 with general solution a*sin(t) + b*cos(t) Now incorpor...

15 days ago | 1

"Busy" when running this code with ODE45
Use ode15s instead of ode45 - you have a stiff system of differential equations.

16 days ago | 0

How to solve complex equations?
Your equation has only complex solutions. Note that you used sind instead of sin in the function definition. Thus assuming the...

17 days ago | 1

| accepted

how to output quantities involving time derivatives in pdepe
You don't have access to the spatial discretization of pdepe, thus no access to the exact time derivatives. But if you choose th...

17 days ago | 0

| accepted

lsqcurvefit claiming complex-valued function
Although the fiting function fit_function is the log of a real valued part of the function, why is it considering that the funct...

17 days ago | 0

How to test whether a function has roots within a given range?
F = @(x) sin(x); xmin = 0; xmax = 6*pi; x = linspace(xmin,xmax); Fx = F(x); I = find(diff(sign(Fx))) F(x(I)) F(x(I+1))

17 days ago | 0

| accepted

how to output quantities within pdefun, bcfun, icfun in pdepe
After pdepe has finished, call "heatcyl" at the output times. x = linspace(0,1,25); t = linspace(0,1,25); m = 1; sol = pdepe...

17 days ago | 2

| accepted

user defined function including a for loop taking vector input giving vector output
p = 10; result = myfunc2(p) function [y2]=myfunc2(p) y2 = zeros(1,p); y2(1) = 1; for t=1:p-1 y2(t+1)=y2(t)+t; ...

18 days ago | 1

how to input variable with set of value to the PDE
The x-interval of integration was missing, and your boundary condition at x = xl would not work. Further, I changed the T-vector...

18 days ago | 0

what's the relation between A , B and M,G for this Nonlinear system of equation ?
A^2 + B^2 - 2*A*B*cos(2*Phi) = G^2 + M^2

19 days ago | 1

| accepted