Torsten
Statistics
0 Questions
6,145 Answers
RANK
15
of 288,464
REPUTATION
15,164
CONTRIBUTIONS
0 Questions
6,145 Answers
ANSWER ACCEPTANCE
0.00%
VOTES RECEIVED
1,756
RANK
of 19,482
REPUTATION
N/A
AVERAGE RATING
0.00
CONTRIBUTIONS
0 Files
DOWNLOADS
0
ALL TIME DOWNLOADS
0
RANK
of 142,569
CONTRIBUTIONS
0 Problems
0 Solutions
SCORE
0
NUMBER OF BADGES
0
CONTRIBUTIONS
0 Posts
CONTRIBUTIONS
0 Public Channels
AVERAGE RATING
CONTRIBUTIONS
0 Highlights
AVERAGE NO. OF LIKES
Content Feed
Solving a system of differential equations with a variable stored in an array
This is what you get with your equations and your data for beta: T = 0:30:360; b = [0.50, 0.49, 0.45, 0.33, 0.27, 0.32, 0.37, ...
1 hour ago | 0
Aligning MATLAB and Mathematica Solutions for the Non-linear Klein-Gordon Equation
All you see in the graphs are floating point errors. Your system of equations has solution x(i) = 0 for 1 <= i <= numBases and...
15 hours ago | 0
| accepted
Defining a time-dependent boundary feature
https://uk.mathworks.com/help/pdf_doc/pde/pde.pdf Chapters: Specify Boundary Conditions > Nonconstant Boundary Conditions Sp...
22 hours ago | 0
How to write a Monte Carlo Simulation?
ps = 120; pb = 30; Profit = @(q,D) ps*min(q,D)-pb*q; number_of_trials = 1000000; expected_profit = zeros(301,1); for i = 0:...
22 hours ago | 0
| accepted
This is an iterative program, but the result is wrong. Why?
g = @(x) cos(x); xc = fpi(g,1,1e-6) g(xc) function xc = fpi(g,x0,tol) xold = x0; error = 2*tol; while error > tol ...
23 hours ago | 0
Solving simple integro-differential equation using dsolve in Matlab
If you insert the line simplify(diff(x,t)+2*x+2*x*t) you will see that this expression turns out to be equal to 0. So MATLAB i...
1 day ago | 0
Integrating the density function of Brownian motion
Maybe it helps to resolve the integration with respect to t: syms t T real positive syms x real g = 1/sqrt(t)*exp(-x^2/t); i...
2 days ago | 0
Helmholtz problem in circular disk
According to your code, data(i,j) is a zero of besselj(j-1,x) in the interval [(i-1)*pi i*pi]. So let's plot besselj(1,x) and b...
3 days ago | 1
How to build a graph using stem for a given equation?
Since your delta function does not support array inputs, you must use something like x1=a1*arrayfun(@(n1)delta(n1-b1),n1) inst...
3 days ago | 1
How do I code a Hamiltonian Cycle in an undirected weighted graph?
I didn't test it yet, but maybe this one from the File Exchange: https://uk.mathworks.com/matlabcentral/fileexchange/51610-hami...
3 days ago | 0
Find critical points of parametric function
The critical points are solutions of a polynomial of degree 6 in z. There is no solution formula for roots of polynomials of deg...
3 days ago | 1
| accepted
Central and Forward difference schemes heat conduction for two-layer materials
heattransfer() function heattransfer xstart = 0; xend = 1; xinterface = 0.25; xmesh1 = linspace(xstart,xinterface,25); xme...
3 days ago | 0
| accepted
Average Optimization using GA or intlinprog algorithms
https://uk.mathworks.com/matlabcentral/answers/2063957-optimization-problem-with-array-input?s_tid=srchtitle
3 days ago | 0
how to ues the exponentials and derivatives in bvp4c
yp1 = deval(S,1) ym1 = deval(S,-1); T1= exp(-b1*(1+m))*ym1(4); T2= exp(-b1)*yp1(4); N1= -(1+b2*(1+m))*ym1(5); N2= -(1+b2)*y...
3 days ago | 0
| accepted
Nonlinear time dependent 2D system of PDE's
I doubt there is a software capable to handle this kind of PDE system automatically. You will have to discretize the spatial d...
4 days ago | 0
How could I have these two erros when I did everything right (otherwise)?
c_p defined in your script does not mean that it is also visible in your functions. Either define the functions as nested funct...
5 days ago | 0
Problem with solving two differential equations depending on each other
Since W0 = 0, you get Re = 0 and thus lambda = something undefined (you divide by Re !). What does the c stand for in dm^3/c fo...
5 days ago | 0
How to plot function with different first x values?
I = @(R)0.5*(R>=0 & R<10)+5/R*(R>=10 & R<=100); U = @(R)0.5*R*(R>=0 & R<10)+5*(R>=10 & R<=100); P = @(R)I(R)*U(R); R = 0:0.1:...
6 days ago | 0
| accepted
Why eigen vector looks different than the calculated one
Eigenvectors are only unique up to a multiplicative constant, and MATLAB's eigenvectors are normalized to be of norm 1: [1;-3]/...
6 days ago | 0
| accepted
Paramter optimization using fmincon
It seems that the model is not able to reproduce your measurement data. The curves look qualitatively different. clc; clear; ...
6 days ago | 0
| accepted
i make a code of shooting method for 3rd order differential equation but i dont get the accurate result
y20 = 2; sol = fsolve(@fun_shooting,y20); [X,Y] = fun_ode(sol); plot(X,Y(:,1)) function res = fun_shooting(y2) [X,Y] = ...
6 days ago | 1
| accepted
Non linear system by Newton-Raphson
x1 = 1; y1 = 1; P = [x1 ; y1]; F= @(x,y) [atan(y)-y.*x.^2 ; x.^2+y.^2-2*x]; JF = @(x,y) [-2.*x.*y (1/(1+y.^2))-(x.^(-2)) ;...
6 days ago | 2
need help getting this to run
Maybe initial_strain_index = find(colum1 == min(colum1),'first'); peak_stress_index = find(colum3 == max(colum3),'first'); in...
7 days ago | 0
How to visualise the 'output' of a Neumann boundary condition?
I think the example "Heat Transfer Problem on a Square" under https://uk.mathworks.com/help/pde/ug/pde.stationaryresults.eval...
8 days ago | 1
| accepted
one of the terms "Rho_net" on the RHS of an set is varying with time over the tspan. How do I implement the set of the ODEs such that the term updates at every time step?
If you insert the line size(dYdt(tspans(1),initial_conditions)) before the call to ode45, you will see that your function dYdt...
9 days ago | 0
| accepted
Newton Raphson gives answers divided by 0
Your linear systems in the course of the Newton iterations cannot be solved since the Jacobian is rank-deficient (it has only ra...
9 days ago | 1
I am wanting to delete these entries of vector x when x = 0. How do I do this. Matlab gives me an error right now saying index cannot exceed 38 but I don't know what's wrong.
Use x(x==0) = []; instead of for i = 1:56 if x(i) == 0 x(i) = [] end end The length of the x-vector wil...
9 days ago | 0
I am trying to create a while loop where a marker can't leave a box of x = -10 and y = 10, y = -10. Once the marker reaches x = 11 I want the loop to stop.
Put the k = k+1 at the end of the while-loop, not at the beginning. And if x(k) == -10, you only set x(k+1), but not y(k+1). Th...
10 days ago | 0
| accepted
solving for y in equation
a = 85.147; h = 55; v = 50; u = 25; f = @(x,y)x-sqrt(a^2+h^2)-(h-y).*cot(asin((h-y)/v))+(h-u)*cot(pi/4-asin((h-y)/v))+v*sin(...
10 days ago | 0
Iteration loop for a function to reach a certain value.
g=9.81; NU=2.52; k=.0293;%from book L=.025;%plate to cover spacing m sigma=5.6697e-8; epsilonp=0.95;%plate emittance epsil...
10 days ago | 0
| accepted