Answered
I am trying to get the answer for the sum to infinity for a geometric series. My output should be 3/2 and a loop should be used. I've done that but it won't work. Please help
Sum2Infinity = 0.0; for i = 0:1:10000 Sum2Infinity = Sum2Infinity + FirstTerm * CommonRatio^i; end

3 years ago | 0

| accepted

Answered
I cant get the correct plotting
% Given p = 7500; % kg/m^3 Cd33 = 16*10^10; %N/m^2 vp= 4620; % Ns/m Ceb = 1.43E-9; % F phi = 2.62; % As/m % Piezoelectri...

3 years ago | 0

| accepted

Answered
Matrix size mismatch Matlab Function
S(j) is a single value (1x1). exp(-norm(r-mu(:,j))^2/(2*sigma^2)) is a vector (50x1). You try to save a vector in a scalar val...

3 years ago | 0

| accepted

Answered
Help with solving a system of coupled equations of motions in time domain?
How did you come up with the sdot vector for the equations so far ? If you use a mass matrix, you simply have to add the additi...

3 years ago | 1

Answered
Error using odearguments ; Error in ode15s (line 153); Error in anand_shrivastan (line 625) ;below is the code
Add dZdt = dZdt.'; as last line in function "SpeciesODE".

3 years ago | 0

| accepted

Answered
How to solve large system of nonlinear equations using fsolve?
n = 20; A = rand(n); B = @(X)X*A-X.*X; X0 = rand(1,n); X = fsolve(B,X0) % or n = 20; A = rand(n); B = @(X)quadratic(X,A)...

3 years ago | 0

| accepted

Answered
Cannot use solver function with poisscdf function
syms lambda n eqn = exp(-lambda)*symsum(lambda^n/factorial(n),n,0,35) == 0.95; solve(eqn,lambda)

3 years ago | 0

Answered
Very slow function, excessive overhead?
The inner while loop is not necessary: if flag(i,2)-flag(i-1,2) == 1, then of course flag(i,2)-flag(i-1,2) ~== -1. Thus the incr...

3 years ago | 0

Answered
Problem during the use of 'fzero'
I'm surprised the initial value x0 is part of your function definition. Are you sure this is correct ? %% Global variables decl...

3 years ago | 1

Answered
How to solve 2nd order BVP using Chebyshev Spectral Method?
alpha = 0.1; xmesh = linspace(0,1,50); solinit = bvpinit(xmesh, [0 0]); bvpfcn = @(x,y)[y(2);0.75/(1+6*alpha*y(2)^2)]; bcfcn...

3 years ago | 0

Answered
Magnetized Hybrid Nanofluid Flow , MATLAB Code has some problem. Please help to Rectify. Highly Appreciated
MHN function MHN % Initialization of paramters beta=1.5; lambda=1.5; wt=2.5; wb=1.5; ks1=0.5; ks2=0.1; a=0.5; epsilo...

3 years ago | 0

| accepted

Answered
optimization in smart energy management system.I have used MILP algorithm but i got infeasible solution no point satisfy the constraints.how could i sort out these issues.
I have used MILP algorithm but i got infeasible solution no point satisfy the constraints.how could i sort out these issues. It...

3 years ago | 0

Answered
how to supply a cost function gradient based on some of optimisation paramters
Would you please tell me if there is a way to supply a gradient based on some optimization parameters and use the finite differe...

3 years ago | 1

Answered
What is the steady state criteria in crank nicolson scheme, especially impilicit finite difference method? and How to stop time level when steady state is reached?
You can define several stopping criteria. Try max(abs(1-CN(:)./COLD(:))) < tol with tol = 1e-4 ... 1e-6

3 years ago | 0

Answered
Nested Function is not working correctly with ode45
You can also get an analytical solution for v since your piecewise equation is easily integrated. EGA215_1234567_A1_Q2 functio...

3 years ago | 0

Answered
Constrained Optimization Problem With Obj Function Which Is Numerical Solution To A System Of ODE's
You never use the 53 solution variables of u and v in your code. And you solve the following ODEs for u and v du/dt = u, u(0)...

3 years ago | 0

Answered
Each element of INTCON must be an integer in the range [1 numVars], ([1 1]). What is the meaning??
f = c' * P_max + c0; % objective function is a scalar value. Thus intlinprog assumes you want to solve a problem with only one ...

3 years ago | 1

Answered
Error using symfun/subsref (line 203) Invalid argument at position 2.
N(i) = symfun(((1/(2*A))*(a(i)+b(i)*x+c(i)*y)),[x y]); instead of N = symfun(((1/(2*A))*(a(i)+b(i)*x+c(i)*y)),[x y]); And N...

3 years ago | 0

Answered
I have these 3 columns and I want to divide p/z and put the p/z values in a new columns
p=[2080,1885,1620,1205,888,645]; z=[0.759,0.767,0.787,0.828,0.866,0.900]; gp=[0,6.873,14.002,23.687,31.009,36.207]; new_row =...

3 years ago | 0

Answered
How to obtained matris vector solution using fminunc function?
You will have to use fminunc three times: c_1=[3;5;8]; c_2=[2;6;8]; x1_sol = zeros(size(c_1)); x2_sol = zeros(size(c_1)); f...

3 years ago | 0

| accepted

Answered
Applying updated initial conditions at every points of tspan and solve system of ODE using ode45
Then why do you solve for v1 if you want to prescribe it ? Simply set v1 = u0x + vp*p1 in your code and determine displacement...

3 years ago | 0

Answered
Using lsqnonlin for deconvolution
lsqnonlin tries to adjust x, and you reset it to 0: function F=insulin(x) ... x=zeros(89,1); This destroys the complete fitt...

3 years ago | 1

Answered
The following error was reported evaluating the function in FunctionLine update: Unable to convert expression into double array.
Your calculations partially produce complex and NaN values. And you cannot prevent this by trying to force them to be real. pv1...

3 years ago | 0

Answered
Solve of implicity equations
You should first check whether a zero of your function really exists. It doesn't seem to be the case: alpha = 0.14; A = 3300; n...

3 years ago | 0

Answered
why do I get this error in the simulation
Before calling "interp1", use % Sort array oxygen_energy_tot in ascending order [oxygen_energy_tot_sorted,idx_sort] = sort(oxy...

3 years ago | 0

| accepted

Answered
ode45 with a array of vector
a1 = ...; a2 = ...; a3 = ...; b1 = ...; b2 = ...; b3 = ...; c1 = ...; c2 = ...; c3 = ...; vp = ...; tValues = linspace...

3 years ago | 0

| accepted

Answered
function handle array Problem
result = integral(@fun,0,1000,'ArrayValued',1) function value = fun(a) m = [1,2]; U(3) = a + 2; V(3) = a - 2; for...

3 years ago | 0

| accepted

Answered
i want to solve a set of homogeneous linear equation
Or use the following code to produce an optimal solution in the least-squares sense: A =1.0e+04*[6.6064,-3.5642,0,0;-3.5642,6.6...

3 years ago | 0

Answered
is there any way to get a 2nd order polynominal through 2 points(each end point)
rng("default") syms a b c x y n = 30; % number of data points xdata = sym('xdata',[1 n]); % x values ydata = sym('ydat...

3 years ago | 0

Answered
Error using plot Vectors must be the same length
As you can see above, S and f have a different number of elements. So they cannot be plotted against each other.

3 years ago | 1