numerical integration and solving for limit
Here is the maximum value you can insert for t: syms x f = 1/sqrt((1-x^2)+0.01/2*(1-x^4)); vpaintegral(f,x,-1,1) And here is...

13 days ago | 0

| accepted

errors on function vs time plot
t = linspace(-1,1,100); v = 12*exp(-8*t).*(t>=0); plot(t,v)

13 days ago | 0

Why doesn't my plot "go to infitnity"?
What finite y value do you want -inf to be in the plot ? Since nobody knows, MATLAB just skips these data points (same for NaN...

13 days ago | 1

k1 is a vector. Thus with the commands dAdt(i) = dAdt(i) + k1.*yita_mn(i,j).*(At(j))*cos(Ot(j)-Ot(i)); dOdt(i) = dOdt(i) + k1....

13 days ago | 0

Double-check for any mistakes
- rho + gamma*(abar(t-1,1) - a(t-1,1)) - gamma*(a(t-1,1) - c(t-1,1)) - delta - (delta*a(t-1,1) + abar(t-1,1)*g_abar - gamma*(aba...

13 days ago | 0

| accepted

2nd order lineair non-homogenous differential equation with EULER not working properly
The number of steps to take to reach tend is (tend-tstart)/dT, not tend/dT. Further, you apply the initial conditions at t = -3...

14 days ago | 0

Solution of the equation
This code doesn't work since there does not seem to be a solution for some or all of the d values you supplied. I suggest you f...

14 days ago | 0

ans = Empty sym: 0-by-1 error in solving equation
syms x y Y = 0:1:100; avd=5; s = 0.25; paco2=40; r= 0.8; hb=14; pb=740; eqn = ((hb*(((((pb-47)*y-paco2/r).^3+150*((pb-47...

14 days ago | 0

Numerical integration of second moment of area for given data
x and y data describe the curves of the upper and lower profile ? https://math.stackexchange.com/questions/3126310/calculating-...

15 days ago | 0

Delete 0's in a matrix
% First case (see above) A = [0,0;0,0;800,1096;0,1461;1462,1701]; A( any(~A,2), : ) = []; A % Second case (see above) A = [...

15 days ago | 1

| accepted

Methods used inside "roots" and "solve" functions to solve polynomial equations.
For numerial computation of polynomial roots, see https://math.mit.edu/~edelman/publications/polynomial_roots.pdf https://de.m...

15 days ago | 0

x = 8e6 + 2e6*rand(1,1000) value = 3e5; x(x>9e6) = x(x>9e6) + value; x

16 days ago | 0

| accepted

How do I remove the Ans from the result
How do I remove the Ans from the result By correctly coding the if-clause: Alpha = [1,0]; if Alpha(1,1) == 1 disp('Ok') e...

17 days ago | 0

How to solve Phase Change Stefan problem with the PDEtool?
There is no chance to solve the problem with a standard tool to solve partial differential equations (like the PDE Toolbox). I ...

17 days ago | 0

| accepted

Plotting a 4D graph with the fourth dimension using color (a 3D graph with the fourth dimension represented using colors)
Use "slice". This will display R0 on predefined planes in 3d space.

17 days ago | 1

Solve a multiobjective optimization problem by problem-based approach in Matlab2021a
% teste_multiobjective file x1 = optimvar('x1','LowerBound',0); x2 = optimvar('x2','LowerBound',0); s1 = 5*x1 + 2.5*x2 >=...

17 days ago | 0

Plotting the function by the points that need to be determined
%% initial conditions global d k0 h_bar ksi m E; Ef = 2.77*10^3; Kb = physconst('boltzmann'); % 1.38*10^(-23) T = 0.12:0.24...

17 days ago | 0

How to find the maximum value of the power and its position on the time axis (x,y)?
s=1.25; t1=0:0.01:pi/2; t2=pi/2:0.01:pi; t=[t1 t2]; p1= 4*(sin(2*t1)).^2.*exp(-t1/s).^2 + 2; p2= 4*(sin(2*t2)).^2.*exp(-t2/...

17 days ago | 1

How to plot these three fucntion X,S and P in one graph with lines
t=0:0.5:14; X0=0.036; Xm=2.24; u0=1.14; S0=41.44; Ms=1.56; Yxs=0.37; Ypx=0.958; dt=0.90; X=(X0*Xm*exp(u0*t))./(Xm-X0+...

17 days ago | 0

| accepted

How to make 'for loop' for three variables (i,j,n)?
for i=1:m-1 instead of for i=2:m-1

18 days ago | 0

goodness of fit data for a fit
Statistical information about the goodness of fit can be obtained by using "fit" instead of "lsqcurvefit" and the second output ...

18 days ago | 0

why plot not generate for the following codes?
creation = 1; syms n alpha = sym(0.03); ketalpha1 = symsum(exp(-abs(alpha)^2 * alpha^n / sqrt(factorial(n)) ), n, 0, 10); ...

18 days ago | 0

| accepted

How to generate N independent, identically distributed random numbers from the Reciprocal Normal Distribution?
% Generate 10 random numbers from the inverse normal distribution % (mu=0, sigma=1) % % G(y) = 1 - F(1/y) % is the cumulati...

19 days ago | 0

| accepted

'integral' with piecewise expressions
L = 1; Beta = @(y,L) pi/2.*((0<=y) & (y<L/2)) + (pi/2+0.3491).*((L/2<=y) & (y<3/4*L)) + (pi/2-0.3491).*((3/4*L<=y) & (y<=L)); ...

19 days ago | 1

| accepted

All the combinations for 3 variables
https://de.mathworks.com/matlabcentral/fileexchange/5475-cartprod-cartesian-product-of-multiple-sets https://de.mathworks.com/m...

19 days ago | 0

| accepted

Plotting the function by the points that need to be determined
Maybe something like this ? %% initial conditions global d k0 h_bar ksi m E; Ef = 2.77*10^3; Kb = physconst('boltzmann'); %...

20 days ago | 0

I don't know whether you and Sona Ghulami are the same person, but i suggest you start with the last code under https://de.math...

20 days ago | 0

plotting a signal which has two equations
None of the two if conditions is satisfied for your t-array. Thus the array S_if does not exist after the for-loop - even if you...

20 days ago | 0

defining a fitting type
My guess is fitType = @(x,a,b,c,d,f)Even_fit(x,a,b,c,d,f,N); fit(x,y,fitType)

20 days ago | 0

| accepted