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

oltre un anno 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...

oltre un anno 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...

oltre un anno 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...

oltre un anno 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...

oltre un anno 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...

oltre un anno 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...

oltre un anno 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...

oltre un anno 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...

oltre un anno 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...

oltre un anno 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...

oltre un anno fa | 0

| accettato

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

oltre un anno fa | 0

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

oltre un anno fa | 0

| accettato

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

oltre un anno fa | 0

| accettato

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

oltre un anno fa | 0

| accettato

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

oltre un anno fa | 0

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

oltre un anno fa | 0

| accettato

Risposto
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 =...

oltre un anno fa | 0

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

oltre un anno fa | 0

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

oltre un anno fa | 0

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

oltre un anno fa | 1

| accettato

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

oltre un anno fa | 1

| accettato

Risposto
Energy Harvesting of Cantilever Beam with MSMA alloy1
I inserted "warning('off')" at the start of your code because your matrix M in your function "system_eqs" is singular and MATLAB...

oltre un anno fa | 0

| accettato

Risposto
Sampling from inverse gamma distribution by specifing mean and standard deviation
It works for mu = 0.5, sig = 0.5, e.g. . It seems that your distribution parameter "a" is too close to the critical value of a =...

oltre un anno fa | 0

| accettato

Risposto
How to code an equation with that integrates a vector.
If there is no closed-form expression for u (e.g. if u depends on x), you can solve the vector-differential equation dx/dt = u ...

oltre un anno fa | 0

Risposto
The surface area won't round to one decimal place and I don't know how to get rid of the 'e's.
https://uk.mathworks.com/help/matlab/ref/fprintf.html https://uk.mathworks.com/help/matlab/ref/format.html

oltre un anno fa | 1

Risposto
Need to correct the following script in part of Infinitely Long Horizontal Cylinder (to be of perfect circular base)
%% 2. Infinitely Long Horizontal Cylinder subplot(1,3,2); theta = linspace(0, 2*pi, 30); x_v = linspace(-10, 0, 20); [Theta,...

oltre un anno fa | 0

| accettato

Risposto
Grey box model of 1d heat diffusion in a rod (official mathworks example)
I think it will be simpler to couple "pdepe" for the solution of your PDE and "lsqnonlin" for optimizing your model parameters. ...

oltre un anno fa | 0

| accettato

Risposto
trouble using pdepe to solve system of pdes
Your initial condition for u(2) at x = 0 is not consistent with your boundary condition. Use function u0 = icfun(x) global ...

oltre un anno fa | 0

| accettato

Risposto
How plot 2d of more than 3 function and substitute the number in each function separate parameter?
y=@(x,t,m)real(((-0.385e1 .* tanh(-x - 0.310e1) - 0.385e1 * coth(-x - 0.310e1)) .^ (0.1e1 ./ m ./ 0.2e1)) .* exp(i * (0.431e1 * ...

oltre un anno fa | 1

| accettato