Respondida
Heat equation using Boundary condition
x=linspace(0,1,20); t=linspace(0,2,10); u=pdepe(0,@eql,@initial,@ba,x,t); surf(x,t,u) title('Surface plot of solution')...

casi 3 años hace | 0

Respondida
unable to get code to run please help
This should be the rough structure of the optimization code. But you don't use your optimization variable P in the objective fun...

casi 3 años hace | 0

Respondida
¿Cómo puedo hacer el campo de direcciones y resolver el PVI y(0)=1 de la siguiente EDO dy/dx=e^(-0.01xy^2)?
The message means that MATLAB cannot find an analytical solution to your differential equation. You will have to use ode45 to ge...

casi 3 años hace | 1

| aceptada

Respondida
how get it?
a = magic(5) b = [1 3 2 4 5] v = arrayfun(@(i)a(b(i),i),1:size(a,2))

casi 3 años hace | 1

Respondida
How do you set up an ode solver with a more than one function in the ode?
If you have two unknown functions (f and g), you need two equations.

casi 3 años hace | 0

Respondida
dimensions of array being concatenated are not consistent
b = [779.3-780*R(k); 76]; instead of b = [779.3-780*R; 76]; But why do you use such a complicated code for 2x2 linear systems...

casi 3 años hace | 0

Respondida
Error using odearguments: returns a vector of length 1, but the length of initial conditions vector is 5.
tspan = [0 100]; SOL0 = [10;0;0;0;0]; [W,SOL] = ode15s(@(w,sol)fun(sol(1),sol(2),sol(3),sol(4),sol(5)),tspan,SOL0); plot(W,SO...

casi 3 años hace | 1

| aceptada

Respondida
Automatic substitution of a symbolic variable with a numerical value
The a in the symbolic equation exp = str2sym('c = a + b'); is not associated with the a = 1 in the first line. This would...

casi 3 años hace | 0

Respondida
Estimating the parameter of a complex equation using maximum likelihood estimation (MLE) method
And this is a probability distribution for arbitrary values of q1, q2, A, lambda, beta1 and beta2 ? I can't believe. If you wer...

casi 3 años hace | 0

Respondida
Output issues with deformation calculation
Use double(DEFAB) etc. to see the result in the usual notation.

casi 3 años hace | 0

| aceptada

Respondida
contour plot with a circular boundary
If you want to plot it within a circular contour, you have a 3d-plot. For a 3d-plot, you can only plot cuts in length direction ...

casi 3 años hace | 0

Respondida
How do I get largest distance between two rows in an nx3 matrix?
max(pdist(X),[],"all") where X is your matrix.

casi 3 años hace | 0

Respondida
interpolation of values gives Inf
InterpCorr.(TitleStr)(1) is equal to 0, so you divide by 0 every time. And my guess is that your line of code should be some ki...

casi 3 años hace | 0

| aceptada

Respondida
How to find the value of x if y=0 from two arrays of data?
Use interp1 if y_temperature_soil is monotonic with depth: y0 = interp1(y_temperature_soil,x_depth,0)

casi 3 años hace | 0

| aceptada

Respondida
How do you solve a coupled ODE when one of the ODE results in a vector of length 3 and the other results in a scalar of length 1?
All solution components have to be aggregated in one big vector z, and also the derivatives have to be supplied in this vector f...

casi 3 años hace | 0

| aceptada

Respondida
How to get best fitting Parameters for a function for all datapoints
You cannot estimate lh and lo separately in your approximating model, only their sum (lh+lo). But your model seems to return c...

casi 3 años hace | 1

Respondida
Linear Interpolation code without using interp1
How can I write a code that works faster? Allocate ynew = zeros(size(n2)) at the beginning of your code. Check whether n is mo...

casi 3 años hace | 1

Respondida
creating smaller matrix from a large matrix
You mean rng("default") A = rand(16); for i = 1:4 for j = 1:4 B = A((i-1)*4+1:i*4,(j-1)*4+1:j*4); A_c...

casi 3 años hace | 0

Respondida
indexing in the nested loop
sx = 0; sy = 0; for s = T-1:-1:0 sx = sx + x^(T-1-s); sy = sy + y^(T-1-s); ratio(s+1) = sy/(sx+K^(1/gamma)*x^(T...

casi 3 años hace | 0

| aceptada

Respondida
nested loop shows only the result of last iteration
... W_t(n+1)=(1+r*(1-tau)).^(-(T-n)).*C_0*epsilon.^T.*K^(1/gamma); for m=(n+1:1:T-1) W_t(n+1)=W_t(n+1)+((C_0.*epsilon.^m...

casi 3 años hace | 0

| aceptada

Respondida
I'm not sure as to why MATLAB is not outputting the figure, as well as the vector [T_Euler]. When I run the code it is fine, but when publishing I don't see either. Any fixes?
%% problem description % Determine the temperature vs time for hot iron with convection and % radiation. HW3_prob4_driver()...

casi 3 años hace | 0

Respondida
I don't know why this code doesn't stops :(
res(3) == 0 and res(4) == 0 are two nonlinear equations in the unknowns yp(3) and yp(4). ode15i seems to have problems following...

casi 3 años hace | 0

Respondida
Delay differential equations where every variable can take on multiple delays
I don't understand your mathematical formulation. Should the summation be over j instead of i ? And shouldn't there be a bracket...

casi 3 años hace | 0

| aceptada

Respondida
Solving Multi-variable Partial Differential Equation
The only thing you can get with the "boundary condition" you gave is the solution on the characteristic curve of the PDE through...

casi 3 años hace | 1

| aceptada

Respondida
Curve fitting with more than 3 variables
Do you really think that there is a causal correlation between the housenumber and the energy from the grid ? So if the house nu...

casi 3 años hace | 1

Respondida
i have some error in my MATLAB code .
B*delta_e(t) is undefined since B and delta_e(t) are both 4x1. Maybe you mean B.*delta_e(t) instead.

casi 3 años hace | 0

Respondida
How to draw ellipse around my data
https://uk.mathworks.com/help/stats/fitgmdist.html After fitting, you can plot a 90% ellipse around your data. What does it me...

casi 3 años hace | 0

Respondida
Event function with DAE system
options = odeset('Mass', M, 'MassSingular', 'yes', 'RelTol', 1e-4, 'AbsTol', 1e-6,'Events',@event); function [position,isterm...

casi 3 años hace | 1

| aceptada

Respondida
I get error in unrecognized function or variable 'T'.
I changed the bold parts below: T = t; for j = 1:length(t) %with drugs T(j+1)=T(j)+dt; M13(j+1) = M13(j)+...

casi 3 años hace | 0

| aceptada

Respondida
Error in ode15s: incompatible sizes arrays, while no errors in function
if Re(i,j)/Re_check(i,j)<0.95 || Re(i,j)/Re_check(i,j)>1.05 fprintf('Re difference too big \n') return end if Re(i...

casi 3 años hace | 0

| aceptada

Cargar más