Respondida
Warning: Non-finite result. The integration was unsuccessful. Singularity likely. > In integral2Calc>integral2t (line 131)
fun(Phi0(1),Phi0(1)) gives NaN + 1i*NaN

alrededor de 1 año hace | 0

Respondida
How to solve rachford - race function using newton method
Multiply the equation by product_{i=1}^{n_c} (1+f_Rg*(K_i-1)). This gives you a polynomial of degree n_c-1 in f_Rg which can be ...

alrededor de 1 año hace | 0

| aceptada

Respondida
How can I evaluate the transition matrix using the fourth order Runge-Kutta ODEs
If it's still of interest: this code seems to work correctly. I don't know why using n = size(A,1); E = A(t+dt/2)*(eye(n)+dt/2...

alrededor de 1 año hace | 0

Respondida
Why do I receive error using BVP5C "The boundary condition function BCFUN should return a column vector of length 1"?
You have one first-order ODE, but two boundary conditions. That's mathematically impossible - you can impose only one condition....

alrededor de 1 año hace | 0

| aceptada

Respondida
Event function and ODE.
Don't use both events simultaneously in your event function. Start the solver with the first event function and integrate until...

alrededor de 1 año hace | 0

| aceptada

Respondida
Trouble solving system of pdes with pdepe
"pdepe" is a solver for parabolic-elliptic partial differential equations. The positions with zero elements in the f-vector indi...

alrededor de 1 año hace | 0

Respondida
Solving the differential equation gives unexpected results. How to fix that?
Better ? v(t) = v(0) - g*t - v0*log(m0/(m0-k*t))

alrededor de 1 año hace | 0

Respondida
Define function for transient transient boundary condition
F1= 3.843e-8; F2= -3.482e-5; F3= 1.179e-2; F4= -1.95; F5= 23.69; Trefrig = @(location,state)F1*state.time^4 + F2*state.time...

alrededor de 1 año hace | 0

Respondida
How to write an event function to stop at second or third x axis crossing.
I'd make the solver stop at each x-crossing. After it returns control to the calling program, add 1 to the number of crossings a...

alrededor de 1 año hace | 0

Respondida
Why don't I need a comment character after an (...) ellipsis?
Most probably there is the convention that ... is the symbol to connect two subsequent lines and that everything behind the ... ...

más de 1 año hace | 1

Respondida
AYUDA EN CODIGO DE MATLAB PARA UN ESTRUCTURA DE PUNETE EN 3D
If you insert the lines % Calcular desplazamientos en los grados de libertad libres rank(K_aa) size(K_aa) before the command...

más de 1 año hace | 0

Respondida
Why Is Colon Operator Not Equivalent To linspace Command?
Both arrays seem to be generated differently such that floating point issues come into play: z1 = 0:.08:50000; z2 = linspace(0...

más de 1 año hace | 0

Respondida
How to draw the function relationship graph between two or three variables, where these variables are solved by a system of equations.
Here are the three possible solutions for VL given tc. % 定义符号变量 VL=20; syms tc t2 Q_out Vco Vcn VL %参数关系 L=100e-9; C=50e-6...

más de 1 año hace | 0

Respondida
Plotting a variable versus space and time
Use "slice": https://uk.mathworks.com/help/matlab/ref/slice.html

más de 1 año hace | 0

| aceptada

Respondida
checkGradients, but the objective function has two inputs: x and xdata?
valid = checkGradients(@(p)semiCircle(p, Q),p0)

más de 1 año hace | 0

| aceptada

Respondida
How can I evaluate the transition matrix using the fourth order Runge-Kutta ODEs
Seems to work. What's your problem ? Don't use the transition matrix to integrate in one pass. If you have to, generate it with...

más de 1 año hace | 0

Respondida
problem with vpa and symsum
The numerical bessel functions are not able to compute your sum correctly (for higher values of k, they return NaN) (see below)....

más de 1 año hace | 0

| aceptada

Respondida
getting singular jacobian error while using bvp4c solver
I can't interprete your solution curves, but it seems to be a problem for the solver that the Scaled Mass goes to 0 with increas...

más de 1 año hace | 0

Respondida
Empty sym: 0-by-1 : Don't understand why
I multiplied by the denominators and squared to get rid of the square roots. Your system seems to have 4 (complex-valued) solut...

más de 1 año hace | 0

| aceptada

Respondida
Non-linear least squares estimation with linear constraints in Simulink?
I vaguely remember that using "lsqcurvefit" within Simulink requires to use the Levenberg-Marquardt algorithm. I can't answer wh...

más de 1 año hace | 1

| aceptada

Respondida
Do I have a convergence issue?
B0 = 0.5; B = fsolve(@fun,B0) function residual = fun(B) g = 10; %m/s^2 Q = 70*1000; %N/m tau_L1 = 24*1000; %Pa t_w = 30...

más de 1 año hace | 0

Respondida
HOw to find displacement integrating velocity?
The equation from above describes the stationary velocity profile over the length L. Thus replacing v by dx/dt does not make sen...

más de 1 año hace | 1

| aceptada

Respondida
Integration symbiolic form in matlab gives error
I think you will have big problems if you try to integrate a function with values in the order of +/- 1e70. syms r theta Part_...

más de 1 año hace | 0

Respondida
Why the symbolic method fail when I tried to expand a ellipse
I get the correct result (at least for coe = [1;1;1;1;1;1]) : aex = expandA(ones(6,1)).' function aex = expandA(coe) A = co...

más de 1 año hace | 0

Respondida
Why did i receive error with my for loop ode45 function?
I used numerical values for L1, L2, mblade, kt and Jblade right from the beginning and got expressions for dL_dtheta2_numeric an...

más de 1 año hace | 0

Respondida
Why do i get the error (Array indices must be positive integers or logical values) when trying to calculate the determinant
I don't get any error. But you shouldn't overwrite the function "det" by using the variable name "det" for the result. Name it ...

más de 1 año hace | 0

Respondida
having if/else statement go to previous input
r = rand(); if opt == 1 || (opt == 3 && r < 0.5) abc else xyz end

más de 1 año hace | 0

Respondida
6x6 system of multivariate quadratic equations ... non-negative real solutions
I think there is no specialized solver for a system of quadratic equations. Thus a general nonlinear solver ("fsolve","lsqnonlin...

más de 1 año hace | 0

Respondida
Color Burr distribution function line
pd = makedist('Burr','alpha',1,'c',2,'k',5); R = random(pd,100,1); pd = fitdist(R,'Burr'); figure(1) plot(pd) x = 0:0.01:2;...

más de 1 año hace | 0

| aceptada

Respondida
Writing ODE equations for ode45 to solve;
The form of the equations that can be solved by ode45 is M*y'(t) = f(t,y). So you can decide whether you choose M to be the iden...

más de 1 año hace | 0

| aceptada