Respondida
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]);

11 meses hace | 0

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

11 meses hace | 0

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

11 meses hace | 0

| aceptada

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

11 meses hace | 0

Respondida
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,:)

12 meses hace | 0

| aceptada

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

12 meses hace | 0

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

12 meses hace | 0

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

12 meses hace | 0

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

12 meses hace | 0

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

12 meses hace | 0

| aceptada

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

12 meses hace | 0

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

12 meses hace | 0

| aceptada

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

12 meses hace | 0

| aceptada

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

12 meses hace | 1

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

12 meses hace | 0

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

12 meses hace | 0

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

alrededor de 1 año hace | 0

| aceptada

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

alrededor de 1 año hace | 1

| aceptada

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

alrededor de 1 año hace | 0

| aceptada

Respondida
Finding x and y values of minimum z in 2-variable function.
z = @(x,y)(x.^2).*(y.^3)-1.4*(x.^2).*(y.^2)-3.91*(x.^2).*y+5.78*(x.^2)-... 11*x.*(y.^3)+15.4*x.*(y.^2)+43.01*x.*y-63....

alrededor de 1 año hace | 0

Respondida
Find best fit curve from multiple curves of different size
I parametrized the curves by normalized arc length. To compute the average curve, I took the mean of R- and Z-coordinates of the...

alrededor de 1 año hace | 0

| aceptada

Respondida
How can I get splines subject to constraints?
This code looks helpful for your purpose: https://uk.mathworks.com/matlabcentral/fileexchange/7976-random-number-from-empirical...

alrededor de 1 año hace | 0

| aceptada

Respondida
How do I edit a loop to be able to store the variables and outputs used in each loop in a way which can be graphed?
Be careful here. In correlations for friction factors, log is usually meant as log to the basis of 10, not the natural logarithm...

alrededor de 1 año hace | 0

| aceptada

Respondida
Question solve(eqn)
Replace solve (f1, f2, f3, f4, f5, f6) by [a0,a1,a2,a3,a4,a5] = solve (f1, f2, f3, f4, f5, f6) But note that extrapolating t...

alrededor de 1 año hace | 0

Respondida
find is not working in a loop, it skips numbers
Never test two double precision numbers for equality. Because of floating point arithmetic, they will almost always be classifie...

alrededor de 1 año hace | 0

| aceptada

Respondida
Can someone rectify this code?
My guess is that you work with an older MATLAB version where script and function parts cannot be mixed. Thus save function F =...

alrededor de 1 año hace | 0

Respondida
elseif that goes through previous statements but with different parameters
if e(t)==1 nn1 = nn11; nn2 = nn21; elseif e(t)==2 nn1 = nn12; nn2 = nn22; elseif e(t)==3 nn1 = nn13; nn2 = nn23...

alrededor de 1 año hace | 0

Respondida
syms to another .m file
Call the other .m-file with the syms-equation (or function) as input argument. syms x f = x.^2; result = fun(f) function res...

alrededor de 1 año hace | 0

Respondida
Clarification regarding Sensitivity Analysis of a coupled ODE system using the odesensitivity function
Start with the simple example y'(t) = p(1)*y(t), y(0) = 1. The solution is y(t) = exp(p(1)*t). The (local) sensitivity is ...

alrededor de 1 año hace | 1

| aceptada

Respondida
How to make ode45 store/return only results of every n-th time step to reduce memory use?
You can specify the times when the solution should be saved in the input variable "tspan" to "ode45". And you should test wheth...

alrededor de 1 año hace | 1

| aceptada