Réponse apportée
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 ...

plus de 2 ans il y a | 0

| A accepté

Réponse apportée
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...

plus de 2 ans il y a | 1

| A accepté

Réponse apportée
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....

plus de 2 ans il y a | 0

| A accepté

Réponse apportée
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...

plus de 2 ans il y a | 1

Réponse apportée
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...

plus de 2 ans il y a | 0

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

plus de 2 ans il y a | 0

| A accepté

Réponse apportée
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...

plus de 2 ans il y a | 0

Réponse apportée
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...

plus de 2 ans il y a | 0

Réponse apportée
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...

plus de 2 ans il y a | 0

| A accepté

Réponse apportée
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...

plus de 2 ans il y a | 1

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

plus de 2 ans il y a | 0

Réponse apportée
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...

plus de 2 ans il y a | 1

| A accepté

Réponse apportée
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...

plus de 2 ans il y a | 0

| A accepté

Réponse apportée
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...

plus de 2 ans il y a | 0

Réponse apportée
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))

plus de 2 ans il y a | 0

| A accepté

Réponse apportée
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...

plus de 2 ans il y a | 3

| A accepté

Réponse apportée
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; ...

plus de 2 ans il y a | 1

Réponse apportée
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...

plus de 2 ans il y a | 0

Réponse apportée
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

plus de 2 ans il y a | 1

| A accepté

Réponse apportée
Why Dsolve considers my equations' variables constants!?
syms t x1(t) x2(t) p1(t) p2(t) eq = [diff(x1, t) == x2, diff(x2, t) == - p2 - x2, diff(p1, t) == 0, diff(p2, t) == p2 - p1] ds...

plus de 2 ans il y a | 0

| A accepté

Réponse apportée
Difference Plot in matlab
You must either interpolate the output values of the first signal to the times of the second signal or the output values of the ...

plus de 2 ans il y a | 0

| A accepté

Réponse apportée
I'm suppose to plot the conventional integral , but I'm getting wrong values
syms t x(t) zeta = 0.1; omegan = 3; f1 = 1; f2 = 0; T1 = 7; T2 = 30; T3 = 40; D2x = diff(x,t,2); Dx = diff(x,t); eqn =...

plus de 2 ans il y a | 0

| A accepté

Réponse apportée
Unable to perform assignment because the left and right sides have a different number of elements.
r0=0.05; k1=0.5; k2=0.5; mu=0.5; rho=0.5; epsilon=0.25; K=1; alpha=0.1; q=0.1; eta=0.05; sigma=5; zeta=0.05; omega0=...

plus de 2 ans il y a | 0

| A accepté

Réponse apportée
How do I solve multiple bvps with coupled boundary values?
Use the facility of bvp4c to solve multipoint boundary value problems as in your case: https://uk.mathworks.com/help/matlab/mat...

plus de 2 ans il y a | 1

| A accepté

Réponse apportée
Solving a nonlinear system of equations (360 unknowns) with constraints
There are specialized algorithms for network simulation. Just blindly writing down the equations and applying a general-purpose ...

plus de 2 ans il y a | 1

Réponse apportée
Stacked Cuboid for PDE Domain
Usually in such cases the domain is divided into subcubes that share common faces. E.g. the upper cube could be divided into 9 ...

plus de 2 ans il y a | 0

Réponse apportée
Solving system of nonlinear equations using Matlab
You equations only have a negative solution for x3 as you can see from the function plot. syms x1 x2 x3 f1 = x1*24.57^2*1.6^2 ...

plus de 2 ans il y a | 0

Réponse apportée
Index exceeds the number of array elements. Index must not exceed 1.
num_samples = size(ITAE.ans , 1); ITAE_value = ITAE.ans(1, num_samples); num_samples is the number of rows of ITAE.ans, but yo...

plus de 2 ans il y a | 0

Réponse apportée
I have a set of non-linear equations and I have to find the values of r1, r2 and r3 in terms of q, where q=n. Please tell me how to solve it.
You have four equations in four unknowns. Usually, such a system has a unique solution and there is no way to express one variab...

plus de 2 ans il y a | 0

| A accepté

Réponse apportée
What happens when we call a function containing a parfor loop within another parfor loop?
This page has the answer to your question: https://uk.mathworks.com/help/parallel-computing/nested-parfor-loops-and-for-loops.h...

plus de 2 ans il y a | 1

| A accepté

Charger plus