已回答
paretosearch does not satisfy nonlinear constraints
You may get feasible initial points if you run this code before calling "paretosearch" (maybe for different x0 to get several fe...

1 month 前 | 1

| 已接受

已回答
When does an ODE integrator reach the best numerical accuracy?
Solve the problem for different error tolerances and see when the solution stabilizes. If available, also test different solve...

1 month 前 | 1

已回答
Plot grraph , srink is problem
I don't know what you mean by "I want to change the Y axis value so that the box and plot on x -axis does concide". Maybe like...

1 month 前 | 0

已回答
double and triple integral for a complex function. I coded but did not get any output
According to the literature, the Monte Carlo Integration method seems to be the usual approach: https://iopscience.iop.org/arti...

1 month 前 | 0

| 已接受

已回答
How to update a function with output from another function?
% Solve the coupled ODEs and PDE [t, y] = ode45(@(t, y) odes(t, y, R, C_p_gas, k_liner, k_CFRP, rho_liner, rho_CFRP, C_p_liner,...

1 month 前 | 0

已回答
Unable to recognise function from called script
The script should be t = [1, 2, 3, 4]; y = lab1_function(t,0.1) the function should be function y = lab1_function(t, c) y...

1 month 前 | 1

| 已接受

已回答
I have matlab code, when i run the program i could 't get the values and also graph.Can u help me in this regards?
Don't solve the equations separately (first for T and C, then for f and h) , but all four together in one call to bvp4c. They a...

1 month 前 | 0

| 已接受

已回答
How to Extract Delayed State Terms in a Model with Distributed Delay?
So in the distributed case you want to extract y2_delayed(t) = integral_{tau = t-1.5]^{tau = t-0.5} y2(tau) dtau thus the te...

1 month 前 | 0

| 已接受

已回答
Solution of a 2nd order non linear implicit differential equation using ode15i implicit solver
% Initial conditions y0 = [0 2]; % initial displacement and velocity %yp0_guess = [0; 0]; % Initial guess for derivatives %...

1 month 前 | 0

已回答
Different behaviour when using symbolic simplification on norm() for scalars and vectors?
If you want to define the norm of the difference of two real-valued n-dimensional vectors, use norm_x_minus_y = sqrt(sum((x-y)....

1 month 前 | 1

| 已接受

已回答
Unable to perform assignment because value of type 'optim.problemdef.OptimizationExpression' is not convertible to 'double'.
If you use optimvars on the right-hand side of the assignment hapsCoreAvailability(coreIdx) = T_start + execTime; (T_start is ...

1 month 前 | 0

已回答
Estimation the parameters of a non-linear equation
If only reproducing the x-P curve matters, you can use the following code. If the density of the x-values in certain regions sh...

1 month 前 | 0

已回答
Find the combination that maximize the objective function values from the data set.
filename = 'Objective Optimization MATLAB.xlsx'; data = xlsread(filename); data(:,7:end) = []; n = 1; while(size(data,1) > 0...

1 month 前 | 0

| 已接受

已回答
Source code for solving second order ODE/PDE using PINN
Why do you want to use code for partial differential equations if you have an ordinary differential equation to solve ? syms ps...

2 months 前 | 0

已回答
To RESHAPE number of elements must not change
It means that the number of elements in the original and the reshaped array must be the same: a = rand(15,1); a_ok = reshape(a...

2 months 前 | 0

已回答
how to modify code for distributed delay
r = 1.5; fun = @(t,x)[r*x(1)*(1-x(2));x(1)]; x0 = [0.5;0]; tspan = [0 20]; [t,x] = ode45(fun,tspan,x0); plot(t,x(:,1)) gri...

2 months 前 | 1

| 已接受

已回答
Finding the closest coordinate from a surface plot based on a X, Y location
lat_GMM and lon_GMM are not within the rectangle in which data for z are supplied. In this case, interp2 returns NaN because i...

2 months 前 | 0

已回答
Global stability of an epidemiology
clear; clc; % Define the parameters of the epidemic model alpha1 = 0.0444; % transmission rate alpha3 = 0.2; % fixed rec...

2 months 前 | 0

| 已接受

已回答
How to fit kinetic model to estimate kinetic parameter from experimental data
I replaced C = b(11)-(0.9*y(2))-(0.947*y(3))-(0.9*y(4)/0.511)-(1.137*(y(1)-1.04)); by C = b(1)-(0.9*y(2))-(0.947*y(3))-(0.9*y...

2 months 前 | 1

已回答
I am trying to solve the system of coupled partial differential equations described in the attachment using the function pdepe. My code runs into this error:
"pdepe" is a solver for parabolic-elliptic partial differential equations (with a second-order spatial derivative term modelling...

2 months 前 | 1

| 已接受

已回答
How to input a novel boundary condition for a coupled PDE system
xOde = [0 1].'; % Define coupling points for ODE definition function [pl, ql, pr, qr] = PDE_PSw_EK_BC_2(xl, yl, xr, yr, t, v...

2 months 前 | 0

| 已接受

已回答
quadprog output: this problem is non-convex
Maybe there are Inf or NaN coefficients in your A ? The following small example works. A = rand(53,66); [m,n] = size(A); opti...

2 months 前 | 0

已回答
Problem with finding the global minimum with fmincon
What are x(1) - x(4) in your mathematical description ? In other words: What are your optimization parameters ? You don't defin...

2 months 前 | 0

已回答
Why do my laptop and work computer produce different results when adding a small term to a matrix?
In the representation with 4 decimal places, the results don't seem to differ. But it's not unusual that computations on diffe...

2 months 前 | 0

| 已接受

已回答
When i call and run this code it just save the Phase1 results in CSvV file not other results .
Use PhaseTable = struct2table(PhaseResults); writetable(PhaseTable,'Phase.csv'); instead of % Write Phase re...

2 months 前 | 0

已回答
Calculate the normal vector between two nodes in the space
n = [0 0 1] points from A to B. So I don't understand what you mean by "normal" to the line connecting two points in 3d. Usual...

2 months 前 | 0

| 已接受

已回答
Parse error at observed data
% Observed data (replace with actual data) % Format: [time, infected, isolated, treated, recovered] observed_data = [ 0, ...

2 months 前 | 0

已回答
Multivariate regression for constrained dependent variables
Define a Nx12 matrix A as A = [x11,x12,x13,x14,x15,x16,x21,x22,x23,x24,x25,x26] and a Nx1 vector b as b = ones(N,1) and sol...

2 months 前 | 0

已回答
Newton's method gives NaN. Can someone improve my code?
Your undamped Newton's method throws you from x = 0.5 to x = 46 appr. in the next step. Here, your function cannot be evaluated ...

2 months 前 | 0

| 已接受

已回答
3D integration over a region bounded by some planes.
f = @(x,y,z)x.^2+y+2*z; N = [100,1000,10000,100000,1000000,10000000,100000000]; Fi = zeros(numel(N),1); Fo = Fi; for i = 1:n...

2 months 前 | 1

| 已接受

加载更多