Respondida
Changing ode solver tolerances on the fly
Is it possible to change (relax) the values of the "Absolute tolerance", Relative Tolerance" on the fly just for the "problem" t...

casi 3 años hace | 1

Respondida
Saving of state just before nonconvergence of ode solvers
The stepsize before the solver exists usually becomes small. So if you choose "tspan" in the call to the ode integrator as a 2-e...

casi 3 años hace | 0

| aceptada

Respondida
ODE45 not producing outputs test does not finish
Use a stiff integrator, namely ode15s instead of ode45.

casi 3 años hace | 0

| aceptada

Respondida
Unrecognized function error in bvp4c (in boundary conditions)
A = 2; B = 1.5; g = load('gr_in.mat'); gr = g.gr; rc = gr(:,2); figure(1) plot(gr(:,2),gr(:,1)) N = 100; r = linspace(rc...

casi 3 años hace | 1

| aceptada

Respondida
Write function for complicated recurrence relation.
I modified your formula slightly, and after this modification, it seems to give correct results: N = 14; L = zeros(N); L(1,1:...

casi 3 años hace | 0

Respondida
Fit ellipsoid to (x,y,z) data
https://uk.mathworks.com/matlabcentral/answers/2038016-how-to-plot-the-best-fitted-ellipse-or-circle?s_tid=srchtitle https://uk...

casi 3 años hace | 2

Respondida
script no longer works after update to 2023b (related to ode45)
Why do you use ode45 for a boundary value problem ? Use bvp4c instead. The below code doesn't give errors or warnings - at leas...

casi 3 años hace | 0

| aceptada

Respondida
how to code this type of objective function (optimization)
Use "fminimax" and read the documentation on how to apply the software also for maximizing the minimum instead of minimizing the...

casi 3 años hace | 1

Respondida
How to plot the best fitted ellipse or circle?
Compute the center of gravity of the point cloud. Call it (x',y'). Compute the point of your point cloud with the greatest dist...

casi 3 años hace | 1

Respondida
Plot Lagrange parabolas for a Simpson's composite rule
clearvars; close all; clc; f =@ (x) 2 + sin((pi/20)*x) + sin((pi/3)*x); x0 = 0;%input("Inferior limit in x to aproximate: ...

casi 3 años hace | 1

| aceptada

Respondida
Solving system of ODE and PDE in 1D
Except for DTwDt(1) = 293.0; which must of couse read DTwDt(1) = 0.0; it seems the equations as written are correctly implem...

casi 3 años hace | 0

Respondida
Error using sym/subs Inconsistency between sizes of second and third arguments.
You cannot subs a vector of size 1x2 for a vector of size 1x3 (see above).

casi 3 años hace | 0

Respondida
Find Coefficients of a 5th order ODE without ode45
syms t x(t) Dx = diff(x,t); D2x = diff(x,t,2); D3x = diff(x,t,3); D4x = diff(x,t,4); D5x = diff(x,t,5); eqn = D5x-5*D3x+4*...

casi 3 años hace | 0

| aceptada

Respondida
How to find mean of sin(t)/t with t [0, pi]?
syms x I = vpa(int(sinc(x),x,0,1))

casi 3 años hace | 0

Respondida
Error in fzero usage
bs_price = blackScholes(sigma,S,K,r,T,r,q); function resu = blackScholes(sigma,S,K,r,T,q) Do you see the different number of i...

casi 3 años hace | 0

| aceptada

Respondida
problem with symbolic factorization with two symbols
Here are the factors: syms s K eqn = s^3 + 10*s^2 + (21+K)*s + 4*K == 0; factors = solve(eqn,s,'MaxDegree',3)

casi 3 años hace | 1

Respondida
My while loop only runs once then hits an error
Maybe like this ? % Section 1 clear; clc; n=1; h(n) = 0; v(n) = 0; t(n) = 0; % initial values m = 0.5; % mass of rocke...

casi 3 años hace | 0

Respondida
Finding sum of for loop with if else statement
You mean this ? trangeF = 0:0.001:1.6; % Time range xt = zeros(size(trangeF)); for n = -2:2 if n == 0 s0 = 312...

casi 3 años hace | 0

Respondida
Solving a system of integro-differential equations
Define dI_n / dt = P1*exp(-lambda_n^2*D_a*t) , I_n(0) = 0 as additional ODEs to be solved. Then I_n = integral_{tau = 0}^{tau...

casi 3 años hace | 0

Respondida
"I'm trying to add an 'if' statement to analyze each result, but Matlab is 'ignoring' the condition and not performing the addition."
c is a vector. Do you want d to be rad2deg(theta) + 180 if all elements of the vector c are < 0 ? Then your code is correct. Oth...

casi 3 años hace | 0

Respondida
Replace a variable in a fuction handle by a function
Define "modelfun" for "nlinfit" not as a function handle, but as a normal function. Within this function, you can do everything ...

casi 3 años hace | 0

Respondida
Using fittype with a symbolic function
The problem was solved by converting the symbolic function "qme" to a numerical function handle by using qme_numeric = matlabFu...

casi 3 años hace | 0

| aceptada

Respondida
fmincon: optimize till nonlinear condition is false
Is there an option for fmincon or other optimizers to define this or maybe to ask the optimizer to go on till the nonlinear cond...

casi 3 años hace | 1

Respondida
i tried to find k which depends on w , but i cannot preallocate p3,p4.....,please help?
Use A=[e15*exp(k),-e15*exp(-k),(c44*p3(j)+e15)*s3(j)*exp(s3(j)*k),(c44*p4(j)+e15)*s4(j)*exp(s4(j)*k),(c44*p5(j)+e15)*s5(j)*exp...

casi 3 años hace | 0

Respondida
Error in bvp4c code
@(r,y)baseBC(ya,yRc,yRp,yRb,P0,We,m,Da,beta,alpha,phi) You don't want to use ya and yb (in your notation: r and y), namely the ...

casi 3 años hace | 0

| aceptada

Respondida
How to plot a current based on active voltage?
i = [100,10,100]; t = [0,0.2,0.4,0.6]; R = 8.314; F = 9.648e4; Temp = 80 + 273.15; %Conductivity reference temperature ipp_...

casi 3 años hace | 0

| aceptada

Respondida
graphs are not displaying
"initial_state" had to be defined. globals had to be included where the variables were needed. "y" has been removed from the g...

casi 3 años hace | 1

Respondida
Differential equation to find activation voltage
Matlab cannot find an analytical expression for u. You have to use a numerical approach. I'm not sure about the term exp(b_an*G...

casi 3 años hace | 1

| aceptada

Respondida
How do I need to configure Matlab's solver such that it returns the correct solutions when complex magnitudes are involved?
Better use solve(expression*expression'==1,w) instead of solve(abs(expression)==1,w) From your hand calculation it follows t...

casi 3 años hace | 0

| aceptada

Respondida
linsolve for coefficient matrix
A.'*X.' = B.'

casi 3 años hace | 1

Cargar más