Community Profile

# Alan Stevens

Last seen: ein Tag ago Active since 2020

#### Content Feed

View by

Runge-Kutta integration of the Taylor-Maccoll eq
Shouldn't k1 = h*F1(t, y1(i), y2(i)); m1 = h*F1(t, y1(i), y2(i)); k2 = h*F2(t+0.5*h, y1(i)+0.5*h, y2(i...

6 Tage ago | 0

| accepted

How can I shade the area between a curve and a vertical line?
Try help fill or help patch Although they refer to polygons, you could easily represent your curve as a multi-segment polyg...

14 Tage ago | 0

Finding best parametric function estimation for ODE of first order
First plot the points, ui vs xi to see what sort of curve it might be. If it looks like it could be a polynomial (as in the cas...

14 Tage ago | 1

Impose a numerical limit on the value of a sum of two vectors
One possibility is min(max(a1),a1+a2)

16 Tage ago | 1

| accepted

How read Matrix Market data into MATLAB?
Do you mean like this: %%MatrixMarket matrix coordinate real symmetric % System: M x_dotdot + E x_dot + K x = B u % y...

16 Tage ago | 0

rounding the imaginary part of the number
Z = 0.0000e+00 + 7.8765e-04*1i; format shortG round(Z,4)

16 Tage ago | 1

| accepted

How can I model second order ODE with matrices and external forcing?
Like this perhaps (I've made up arbitrary data; you will obviously have to replace it with your own) x0 = [1; -1]; v0 = [0; 0]...

16 Tage ago | 0

Randomisation of ramberg osgood equation
You don't need a loop. Yu can just do deltaepsilon = r/E + 2*(r/(2*Kp)).^(1/np); where Kp = K' and np = n'. Note that you nee...

20 Tage ago | 0

How to solve a ode within a for loop?
Like this A = 0:0.05:0.35; tspan = 0: 0.1: 1; p32 = zeros(length(tspan),length(A)); r = 0.7; r0 = 2; p0 = [r; r0]; for k ...

22 Tage ago | 0

| accepted

How can I solve a system of equations with exponential terms in Matlab?
There is almost certainly no analytical solution here! For a numerical solution (given the values of the other parameters), try...

23 Tage ago | 0

To find integral of a function in loop
Try replacing q = integral(f,6e-7, x); with q(i) = integral(f,6e-7, x(i));

26 Tage ago | 0

| accepted

Numerical integration of the missile dynamic model
I don't know the actual problem, but if you look at your values of Y against time you will see that some of them are clearly div...

26 Tage ago | 0

| accepted

Monte Carlo simulation for probability of a multi part machine working
Here's one possibility: % each A part works 98% of the time. Each B works 91% of the time. % The machine has 12A and 7B parts....

26 Tage ago | 0

| accepted

Need help on using trapz on a definite integral
Replace Final = Nach + Vor; by Final(i) = Nach + Vor;

27 Tage ago | 0

Multiple parameter Non Linear Curve Fit- Error using ^
You probably need to use element by element multiplication and division etc in your definition of fun. I.e. use a.*b and a...

28 Tage ago | 0

| accepted

How to separate plot titles on single line? [solved]
Define the titles with an extra space at the end of each word.

28 Tage ago | 0

Plotting to test probability convergence

29 Tage ago | 0

| accepted

How to do integration of this equation and plot graph between T and y
Like this perhaps: u=0.020; x=0.025; z=0; t=5; eta= 0.6; P=200; C=500; ro= 1000; l=0.008; Alpha=17e-6; r = 0.01; % Gu...

29 Tage ago | 0

| accepted

Making relative histogram for data, and a convergence diagram to see if enough trials have been run for monte carlo simulation.
Something like the following? N = 1000; Throwsum = zeros(N,1); Relativefreq = zeros(N,1); count = 0; for i=1:1000 throwD...

29 Tage ago | 0

| accepted

Runge Kutta and Binary Search to analyze system of ODEs
At the moment you have dj/dt = -K*d0^2; which means rate of reduction of radius is constant. Since d0^2 is a constant this m...

30 Tage ago | 0

| accepted

Gaussian Curve Plotting Issues
If you have 10000 individuals with a probability of 0.013 of being in group A, you would expect your curve to peak close to 130 ...

etwa ein Monat ago | 0

Solve equation with 2 unknowns using solve
There are an infinite number of solutions! Restrict the range. One possibility is as follows bfn = @(a) acos(1.6-cos(a)); fn...

etwa ein Monat ago | 1

Is this formula correct for simulate 3D Water rocket
Probably vy = v0*cos(alpha); should be more like vy = v0*cos(alpha)*cos(beta); where beta is an angle that will point the r...

etwa ein Monat ago | 0

Need to Produce a different size matrix on each iteration
Try the following wqe=[1;2;3;4;5;6;7;8;9;10;11;12;13;14;15]; kli=[51;52;53;54;55;56;57;58;59;60;61;62;63;64;65]; hsd=size(kli...

etwa ein Monat ago | 0

| accepted

Getting error using syms to implement power formula
Like this? % Arbitrary values - replace with your own fs = 1; x_h = @(n) 1/(n+10)^2; %Power formula syms N...

etwa ein Monat ago | 0

| accepted

Solving a system of 12 ode’s with discrete data using Euler method
Look up interp1 doc interp1

etwa ein Monat ago | 0

| accepted

How to save all variables in for loop

etwa ein Monat ago | 0

How to handle a variable inside matrix without using syms toolbox?
Are you looking for somethig like this; f = @(w) [1,-0.1409*w^2;0,1]*[1,0;1/286.48,1]*[1,-0.05493*w^2;0,1]... *[1,0;1/1793...

etwa ein Monat ago | 0

| accepted