# Alan Stevens

**Statistics**

RANK**56**

of 239 438

REPUTATION**2 456**

CONTRIBUTIONS

**0** Questions

**923** Answers

ANSWER ACCEPTANCE **
0.00%
**

VOTES RECEIVED**231**

RANK

of 16 816

REPUTATION**N/A**

AVERAGE RATING**0.00**

CONTRIBUTIONS**0** Files

DOWNLOADS **0**

ALL TIME DOWNLOADS**0**

CONTRIBUTIONS**0 Posts**

CONTRIBUTIONS**0** Public Channels

AVERAGE RATING

CONTRIBUTIONS**0** Highlights

AVERAGE NO. OF LIKES

**Content Feed**

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...

4 dagar 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...

12 dagar 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...

12 dagar ago | 1

Impose a numerical limit on the value of a sum of two vectors

One possibility is min(max(a1),a1+a2)

14 dagar 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...

15 dagar ago | 0

rounding the imaginary part of the number

Z = 0.0000e+00 + 7.8765e-04*1i; format shortG round(Z,4)

15 dagar 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]...

15 dagar 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...

18 dagar 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 ...

20 dagar 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...

21 dagar 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));

24 dagar 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...

25 dagar 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....

25 dagar ago | 0

| accepted

Need help on using trapz on a definite integral

Replace Final = Nach + Vor; by Final(i) = Nach + Vor;

26 dagar 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...

26 dagar 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.

26 dagar ago | 0

Plotting to test probability convergence

See my reply to your question at https://uk.mathworks.com/matlabcentral/answers/871048-how-to-do-integration-of-this-equation-an...

27 dagar 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...

27 dagar 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...

27 dagar 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...

28 dagar 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 ...

29 dagar 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...

ungefär en månad 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...

ungefär en månad 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...

ungefär en månad 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...

ungefär en månad ago | 0

| accepted

Solving a system of 12 ode’s with discrete data using Euler method

Look up interp1 doc interp1

ungefär en månad ago | 0

| accepted

How to save all variables in for loop

Try load(i,:)=com(i,1)*sin(com(i,2)*t);

ungefär en månad 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...

ungefär en månad ago | 0

| accepted

Issues in For loop

Have you intialised j (i.e. set j = 1;) before going into the for i = 1:38976 loop?

ungefär en månad ago | 0

How can I randomize cube positions in matlab?

Like this? E=1; N=3; Nsamples=10; d=randi(N,1,Nsamples); s=randi(2,1,Nsample...

ungefär en månad ago | 0

| accepted