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

Rolle''s theorem, The curve (y+1) 2 =x 2 passes by the points (1, 0) and ...

The curve (y+1) 2 =x 2 passes by the points (1, 0) and (- 1, 0). Does Rolle's Theorem clarify the conclusion that  dy dx  vanishes for some value of x in the interval -1≤x≤1?

Trignometery., using the formula sin A =under root 1+ cos2A /2 . find value...

using the formula sin A =under root 1+ cos2A /2 . find value of 30 degree, it is being given that cos 60 degree =1/2.

Derive a boolean first-order query, Consider a database whose universe is a...

Consider a database whose universe is a finite set of vertices V and whose unique relation .E is binary and encodes the edges of an undirected (resp., directed) graph G: (V, E). Ea

Draw tangent graph y = sec ( x ), G raph y = sec ( x ) Solution: As wi...

G raph y = sec ( x ) Solution: As with tangent we will have to avoid x's for which cosine is zero (recall that sec x =1/ cos x) Secant will not present at

Calculus, find or evaluate the integral integrate((e^2x + e^x + 1)/(e^x))dx...

find or evaluate the integral integrate((e^2x + e^x + 1)/(e^x))dx

Describe simplifying fractions with example, Describe Simplifying Fractions...

Describe Simplifying Fractions with example? When a fraction cannot be reduced any further, the fraction is in its simplest form. To reduce a fraction to its simplest form, div

Partial Differentiation, If the sides angles of a triangle ABC vary in such...

If the sides angles of a triangle ABC vary in such a way that it''s circum - radius remain constant. Prove that, da/cos A +db/cos B+dc/cos C=0

Determine the area of the shaded region, The diagram below shows the cross ...

The diagram below shows the cross section of a pipe  1/2  inch thick that has an inside diameter of 3 inches. Determine the area of the shaded region in terms of π. a. 8.75π i

Rationalize the denominator, Rationalize the denominator for following.  Su...

Rationalize the denominator for following.  Suppose that x is positive. Solution We'll have to start this one off along with first using the third property of radica

Continuous Probability Distributions, Ask questioOn average, Josh makes thr...

Ask questioOn average, Josh makes three word-processing errors per page on the first draft of his reports for work. What is the probability that on the next page he will make a) 5

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