Initial value problems, math, Mathematics

Assignment Help:
Write a Matlab function MyIVP that solves an initial-value problem (IVP) for a system of ordinary
differential equations (ODEs) of the form
x ?(t) = f (t, x(t)),
where f : R × Rn ? Rn is an arbitrary function with one one-dimensional input (for time t) and one n-dimensional input x, and n-dimensional output. The function should implement a Runge-Kutta formula (for example, the rk34 formula or the Dormand & Prince formula).
The initial value x0 is provided by the user of MyIVP. The first line of MyIVP (saved in a file MyIVP.m) should look like this
function [xend,t,xt]=MyIVP(f,x0,tspan,N) Inputs
• f: function defining the right-hand side of the ODE. f should accept two arguments: t (a number) and x (an n-dimensional vector). The function f should return an n-dimensional vector y (the time derivative). Typical calling sequence: y=f(t,x), returning the value of f at time t in position x.
• x0: initial value where integration starts from (n-dimensional vector).
• tspan: Starting time and end time for integration. Integration has to run from time t =tspan(1)
to time t =tspan(2).
• N: number of steps for integration. The integration stepsize h=(tspan(2)-tspan(1))/N should
be small.
Outputs
• xend: result of integration at t =tspan(2).
• t: vector of times at which intermediate values have been computed (this should have N + 1
entries).
• xt: intermediate values (n × (N + 1)-array). xt(:,k) should be the solution at t(k).
You can check the built-in variable nargout inside your function to see if the user wants to get three outputs or only the end value xend. If nargout==1 you don’t need to store the intermediate values.

Related Discussions:- Initial value problems, math

Show that positive integers is divisible by 6, Show that the product of 3 c...

Show that the product of 3 consecutive positive integers is divisible by 6. Ans: n,n+1,n+2 be three consecutive positive integers We know that n is of the form 3q, 3q +1

3-d geometry, Q) In 3D-geometry give + and - signs for x,y,z, in all eight ...

Q) In 3D-geometry give + and - signs for x,y,z, in all eight octants Ans) There is no specific hard rule for numbering the octants. So, it makes no real sense to ask which octan

Measures of dispersion- measures of central tendency, Measures of Dispersio...

Measures of Dispersion - The measures of dispersion are extremely useful in statistical work since they indicate whether the rest of the data are scattered away from the mean

Arithmetic progression, the radii of circular base of right circular cylind...

the radii of circular base of right circular cylinder and cone are in the ratio of 3:4 and their height are in the ratio of the 2:3 what is the ratio of their volume?

Geometric progression (g.p.), Learning geometric progression ...

Learning geometric progression vis-á-vis arithmetic progression should make it easier. In geometric progression also we denote the first t

Regression, Regression line drawn as Y=C+1075x, when x was 2, and y was 239...

Regression line drawn as Y=C+1075x, when x was 2, and y was 239, given that y intercept was 11. calculate the residual

Quadratic equation, If roots of (x-p)(x-q) = c are a and b what will be th...

If roots of (x-p)(x-q) = c are a and b what will be the roots of (x-a)(x-b) = -c    please explain? Ans) (x-p)(x-q)=c x2-(p+q)x-c=0 hence,   a+b=p+q  and      a.b=pq-c

AREA, How do you find the distributive property any faster?

How do you find the distributive property any faster?

What is a lattice, What is a lattice? Which of the following graphs are lat...

What is a lattice? Which of the following graphs are lattice and why? Ans:  Let (L, ≤) be a poset. If each subset {x, y} consisting of any two elements of L, comprises a glb (I

Write Your Message!

Captcha
Free Assignment Quote

Assured A++ Grade

Get guaranteed satisfaction & time on delivery in every assignment order you paid with us! We ensure premium quality solution document along with free turntin report!

All rights reserved! Copyrights ©2019-2020 ExpertsMind IT Educational Pvt Ltd