Risposto
How to graph coupled differential equations
% Parameters m1 = 4.0; m2 = 6.0; L = 1.5; k = 100.0; % Added missing parameter g = 9.81; F0 = 100; tF = 1; % Initial st...

circa un mese fa | 0

Risposto
Nonlinear Curve fitting with integrals
It takes too long to run the code online - try if it produces reasonable results. Maybe you have to set lower and upper bounds f...

circa un mese fa | 0

Risposto
Matlab function that tells the full path of from where the function has been called
What about passing the full path where the function is located from the respective functions 2,3,4,5... to function1 ? See htt...

circa un mese fa | 0

Risposto
how to deal with discontinuities in solutions to differential equations using ode45
Usually, the differential equations to be solved are much more complicated than yours which makes it impossible to foresee possi...

circa un mese fa | 2

| accettato

Risposto
How to rename a new matrix obtained by adding a row to the previous matrix
A = [1:5; 6:10] B = [A;11*ones(1,size(A,2))]

circa un mese fa | 0

| accettato

Risposto
eventfunction for stopping solver when output becomes complex doesn't work
Try options = odeset('Events', @(t,y)stop(t, y, Vb0, lambda, mb, rhog0, Cw)); [t, y] = ode45(@(t, y) DE(t,y, Vb0, lambda, mb,...

circa un mese fa | 0

Risposto
Drawing Bessel function is not possible in that it give me an error
Maybe this is what you want: x=(0:0.01:5).'; i=0:20; hold on for m=0:5 j=sum((-1).^i./(factorial(i).*gamma(i+m+1)).*(x/...

circa 2 mesi fa | 0

| accettato

Risposto
Facing problems in nonlinear system
"pdepe" gives a different solution than your method. Maybe the nonlinear systems have multiple solutions. %%%% Problem_01 %%%% ...

circa 2 mesi fa | 0

| accettato

Risposto
How to code a 1D ODE in MATLAB
Usually, the equation reads dW/dt + U*dW/dx = C (1) If U in this equation is positive, your flow goes from left to right,...

circa 2 mesi fa | 0

| accettato

Risposto
lsqcurvefit Surface Fitting Troubleshooting
There is a problem with your "poly_fun". You have to use poly_fun = @(c,XY) c(1)*XY(:,:,1).^1 + c(2)*XY(:,:,2).^1 + c(3)*XY(:,:...

circa 2 mesi fa | 2

| accettato

Risposto
I have this system of linear equations and i need to solve it for unknowns. how i can solve it in MATLAB? we have 8 unknowns and 8 equations.
@javeria If you want to avoid writing your system as A*X = B, you can use a nonlinear solver like "fsolve". As the system is l...

circa 2 mesi fa | 0

Risposto
Finding one real solution with vpasolve
lambda = vpasolve(LCG/B*1/lambda_t==0.75-1/(5.236*CV^2/lambda_t^2+2.4) == 0, lambda_t, [0,inf]);

circa 2 mesi fa | 0

Risposto
make function from matrix
t = linspace(0,10,3000).'; a = rand(3000,1); f = @(x)interp1(t,a,x); f(8.336)

circa 2 mesi fa | 0

Risposto
I am struggling to format and run a pde as part of a liquid-liquid packed column
K_M = 0.1252; xb_meth = 0.0054; Rho_FAME0 = 0.915; Phid_cho = 0.2352; Vd = 6.1006e-4; kod_meth = 1.0157e-5; ap = 204; Vc ...

circa 2 mesi fa | 0

| accettato

Risposto
trouver un angle avec l'equation de trajectoire d'un projectile
x0 = 0.08; y0 = 0.05; v0 = 2; g = 9.81; syms t x theta eqn1 = x == v0*t*cos(theta) + x0; eqn2 = 0 == v0*t*sin(theta) - 0.5...

circa 2 mesi fa | 0

Risposto
using bvp5c evaluate velocity and temperature at specific points
yeval = deval(S,[-1,-0.4,0,0.4,1]); veval = yeval(1,:) thetaeval = yeval(2,:)

circa 2 mesi fa | 0

| accettato

Risposto
Verifying my code for this problem, are these the right solution and graphs?
Look at the contour plot. The temperature at x = 0 should be at T_b = 50°C. This is not the case in your graphics. Maybe you int...

circa 2 mesi fa | 0

Risposto
Optimization method based on the Nonlinear least squares as well as fmincon for the defined objective function (Non-linear function).
You can add a multiple of the vector you obtain by the below command "double(null(A))" to "estimated_T", and you will still get ...

circa 2 mesi fa | 0

Risposto
Error in ode45 and 'Events' of a Satellite Tracking Radar Model
The "tspan" vector you provide for "ode45" has no influence on the internal time steps taken by the integrator. In zones with st...

circa 2 mesi fa | 0

Risposto
I have been trying to derivate equation under to find the formula for 3 parameters Weibull distribution as same as picture. Can someone please help me fix it?
This code works, but it will not give you the result you want. syms m sigma_i sigma_0 sigma_th i N lnL = symsum(log((m/sigma_0...

2 mesi fa | 0

Risposto
How to vectorize analytical derivative functions for bvp4c?
Only the computation of the vector-valued derivatives function f can be vectorized, not the computation of its analytical deriva...

2 mesi fa | 0

| accettato

Risposto
can't run applyBoundaryCondition
"applyBoundaryConditions" is used for general PDE models. You want to use the "thermal" model. The available properties and ob...

2 mesi fa | 0

Risposto
pdepe boundary condition help!!
Note that the f you define in "pdefun" doesn't fit the formula from the paper. You forgot to multiply both terms with m_d/lambda...

2 mesi fa | 0

| accettato

Risposto
I must solve ode equation system but dsolve gives me error
The third boundary condition is wrong as stated. You must use dv/ds at some position s (e.g. 200). %Geometric and Inertial Par...

2 mesi fa | 0

| accettato

Risposto
How can I change a variable and collect an vector of outputs?
V_array = 0.5:0.1:20; % Material Parameters c0 = 0.2; % Alloy composition k = 0.28; % Partition Coefficient (conc solute in s...

2 mesi fa | 1

Risposto
MLE of combined analytically and empirically defined distributions
MATLAB's "mle" accepts custom pdf's. If you write b as sin^2(p) and (1-b) as cos^2(p), it should work to estimate mu, sigma an...

2 mesi fa | 0

Risposto
Plot function for different values of two coupled variables
k5 = [0.1 0.2 0.3]; k6 = [0.3 0.5 0.7]; colors = ['b','r','g']; N = 3; hold on for i = 1:N [time,p,m] = fun(k5(i),k6(i...

2 mesi fa | 0

Risposto
I need to create a cicle in which, for a specific value, it solves one set of equation or another. It doesn't compute.
Works for me. Maybe you use an older MATLAB release ? Lcz = 10; %Contact lenght [mm] delta0 = 25; %Softening Lenght [mm] delt...

2 mesi fa | 0

| accettato

Risposto
Variable might be set by a nonscalar operator
%% Define varaibles theta_deg=10:1:80; %% theta in 1 degrees increments theta_incline=18.05; theta_decline= 40.65; M=1; % ma...

2 mesi fa | 1

| accettato

Risposto
Can I specify different thermal conductivities for different regions in a space using Matlap PDE Toolbox? For example, the space is divided into a thermal conductivity space o
There are two ways: A good method is to generate the geometry such that it consists of a certain number of zones to which you c...

2 mesi fa | 0

| accettato

Carica altro