Answered
Trying to find area under my plot
https://de.mathworks.com/matlabcentral/answers/514107-how-to-find-area-enclosed-by-a-closed-curve To get the positive and negat...

3 years ago | 0

Answered
I want to use a matrix and a function to then create another function
theta1 = -50; theta2 = 30; theta3 = 15; theta4 = -75; theta5 = -45; theta6 = -60; DH = [theta1, 486.5, 150, -90;... ...

3 years ago | 0

Answered
Solving system of 9 nonlinear equaitons in 16 variables
x0 = -10*ones(16,1); AB = [13 15 -1;9 24 0;7 2 35]; options = optimset('TolFun',1e-16,'TolX',1e-16); x = fmincon(@(x)fun(x,AB...

3 years ago | 2

| accepted

Answered
the same value resulted from the linear optimization problem
You define w(k)=beta(k)*(-k_(k) *w_k(k)-k_(k)*log2(lamda(k)*(1-tau^2))*Pt + 2*k_(k)*log2(sqrt(segma_squared))-v_k(k)); and w...

3 years ago | 0

Answered
Hello guys, I need to solve this system of nonlinear symbolic equation; the variables known are X,Y,Z and I already know the solution for q1; so I just need to find q2 and q3.
syms d1 d5 a2 a3 a4 q1 q2 q3 X Y Z % Square eq1 and eq2, add both sides of the equations and take sqrt % Check for other solu...

3 years ago | 0

| accepted

Answered
What's the simplest way to determine the mass of a sphere from a plot of its density varying in r?
If r and rho are column vectors where r is ordered as r(1) < r(2) < ... < r(n) and rho is the density at position r, an approxim...

3 years ago | 2

| accepted

Answered
Index exceeds the number of array elements. Index must not exceed 1.
Pb is a scalar. So Pb(ii) does not exist. You can use Fbump(ii)=Pb/VA(ii); %this line of code is givig me the error but I don...

3 years ago | 0

Answered
Erro in calling function using horner()
I don't understand what you intend by using the lines [pz,b] = horner(a,x); [dpz,b]= horner(b,x); [pz ,a] = horner(a,x) The...

3 years ago | 0

| accepted

Answered
Solving an integral when the parameters is a matrix - Controllability Gramian
You mean syms t t0 tf A = [0 1 0 0;0 0 -4.91 0;0 0 0 1;0 0 73.55 0]; B = [0;2;0;-10]; [M,J] = jordan(A); phi_s=M * expm(J...

3 years ago | 0

| accepted

Answered
sum of Bessel function
syms x syms m integer A = 2/((2*m+1)*pi)*besselj(0,pi/2*m)*sin((2*m+1)*1000*x-(2*m+1)*pi/2) g = symsum(A,m,0,Inf) fplot(g...

3 years ago | 0

Answered
Solving Laplace equation using Finite difference
If L = 3 and H = 1, you cannot choose the same number of points in x and y direction if you update U as if you use a grid with e...

3 years ago | 0

| accepted

Answered
solve matrix differential equations with ode45
syms t y1(t) y2(t) A = [0 1;8 -2]; Y = dsolve(diff([y1;y2],t)==A*[y1;y2]) y1 = Y.y1; y2 = Y.y2; diff(y1,t)-y2 diff(y2,t)-8...

3 years ago | 0

Answered
How do I use meshgrid and quiver3 in a system of differential equations?
https://web.sha1.bfh.science/Labs/F2/ODE/VectorFields.pdf

3 years ago | 0

Answered
How can I solve this matrix question using for loop?
Not the most efficient solution, but maybe best to understand what's happening. A=[0 0 0 ;6 5 7;12 14 30;35 40 42;75 55 50;85 6...

3 years ago | 0

| accepted

Answered
myfunction not working when applied on a meshgrid
x= 0:0.01:2*pi; %629 element. y= 0:0.001:1; %1001 element. [X, Y]= meshgrid(x,y); z=zeros(1001,629); %% for X and Y are bot...

3 years ago | 1

| accepted

Answered
How to find the actual value of d,b,t and D in the maximum value of sac?
close all clear all clc % Constants c = 340; % Speed of sound rho = 1.204; % ...

3 years ago | 0

Answered
Need to solve this with Runge-Kutta 4th order method in MatLab.
fy=@(x,y,z) z; fz=@(x,y,z) (1-y^2)*z-y; x(1)=0; z(1)=0; y(1)=1; h=0.1; xfinal=10; N=ceil((xfinal-x(1))/h); for j=1:N ...

3 years ago | 0

| accepted

Answered
If statement should continue for some values even if "else" is true
You must have a flag which indicates whether the engine is on or off. function [w,T] = EnginePower(B_SoC) persistent on ...

3 years ago | 0

Answered
Make a loop with a function to vary one value
tspan = [0 :0.01:20]; y0 = 0; opt = odeset('RelTol',1e-12,'AbsTol',1e-16,'MaxStep',0.001); Eta = [10, 50, 100, 500, 1000]; f...

3 years ago | 0

| accepted

Answered
Verify the Divergence theorem in MATLAB
Nice exercise. syms X Y Z x y z syms a phi positive F = [X^2*Y^2,Y^2*Z^2,Z^2*X^2]; divF = diff(F(1),X)+diff(F(2),Y)+diff(F(...

3 years ago | 0

Answered
How to get the derivate using bvp4c
I assumed Alpha = A in your code. clc; p=0.01; Betaf= 207; Betas = 17; Beta = 0.5; kof = 0.613; kos = 400; m = 1; b2 = ...

3 years ago | 0

| accepted

Answered
How can I more fully vectorize to eliminate for loops?
Wans2 = (sum(G.*(W.'-V),2)).'

3 years ago | 0

| accepted

Answered
Trying to use fsolve works over lower part of range (values < 10) but fails in upper part of range (values > 10)
It's very important that you have a good initial guess for pHv given Vbv. So looping over the complete range of Vbv and taking ...

3 years ago | 1

| accepted

Answered
List of all Binomial permutations with a certain number of events
n = 8; A = dec2bin(0:2^n-1)-'0' A4 = A(sum(A,2) == 4,:) 8*7*6*5/(4*3*2*1) % more general Ai = arrayfun(@(i)A(sum(A,2) == i,...

3 years ago | 0

Answered
Code for solving first-order ode with symbolic initial condition
syms y(h) k ode = diff(y) == k*y; sol_ode_general = dsolve(ode) vars = symvar(sol_ode_general) eqn1 = subs(sol_ode_general,v...

3 years ago | 0

| accepted

Answered
How to convert PDE into matrix form
Can you write this in matrix form ? m1*diff(x1(t), t, t) + 0*m2*diff(x2(t), t, t) + (c1 + c2) * x1d - c2 * x2d + (k1 + k2) * x...

3 years ago | 0

Answered
Why do I receive the error in horzcat?
blkdiag(ddf(x), zeros(n)) is 2x2, A' is 1x1. Thus it's not possible to concatenate the two horizontally in the command H...

3 years ago | 0

Answered
Determine the coefficients in symbolic
I don't know why the coefficients appear in the wrong order. syms In In_1 On_1 On = 91/100*In-7/8*In_1+1/2*On_1; coeffs(On,[I...

3 years ago | 0

| accepted

Answered
how to solve non-linear equations in a nozzle
Maybe you can give better initial guesses in x0 for the solution than I can ... x0 = 10*rand(7,1); options = optimset('Max...

3 years ago | 1

| accepted

Answered
Solving Nonlinear System of Equations with 10 variables
W=100; s=10; S=23; %F=13; %m=9; fun = @(X1,X2,X3,X4,X5,X6,X7,X8,X9,X10)... [5-(W+S*9)*X1+S*9*X2-S*20*X1^2+S*20*X...

3 years ago | 0