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

1 year ago | 0

| accepted

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

1 year ago | 0

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

1 year ago | 0

| accepted

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

1 year ago | 0

| accepted

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

1 year ago | 0

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

1 year ago | 0

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

1 year ago | 0

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

1 year ago | 0

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

1 year ago | 1

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

1 year ago | 0

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

1 year ago | 0

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

1 year ago | 0

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

1 year ago | 0

| accepted

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

1 year ago | 0

| accepted

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

1 year ago | 0

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

1 year ago | 0

| accepted

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

1 year ago | 0

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

1 year ago | 0

| accepted

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

1 year ago | 1

| accepted

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

1 year ago | 0

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

1 year ago | 1

| accepted

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

1 year ago | 0

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

1 year ago | 0

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

1 year ago | 0

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

1 year ago | 0

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

1 year ago | 0

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

1 year ago | 0

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

1 year ago | 0

| accepted

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

1 year ago | 0

| accepted

Answered
How to use a solver to get the terminal velocity of a sphere.
v0 = fsolve(@fun,1) function res = fun(v0) d_p =...; rho_air = ...; nu_air = ...; rho_w = ...; g = ...; Re ...

1 year ago | 1

| accepted