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

Polynomials, zeroes of polynomial 2x2-3x-2

zeroes of polynomial 2x2-3x-2

Solve 2 ln (x) - ln (1 - x ) = 2 single logarithm, Solve 2 ln (√x) - ln (1 ...

Solve 2 ln (√x) - ln (1 - x ) = 2 . Solution: Firstly get the two logarithms combined in a single logarithm. 2 ln (√x) - ln (x  - l) = 2 ln ((√x) 2 ) ln (1 - x ) = 2

Normal to y=f(x) , If the normal to y=f(x) makes an angle of pie/4 with y-a...

If the normal to y=f(x) makes an angle of pie/4 with y-axis at (1,1) , then f''(x) is eqivalent to? Ans) The normal makes an angle 135 degree with the x axis. also f ''(1)

Demonstrate that dijkstra algorithm - digraph, Demonstrate that Dijkstra's ...

Demonstrate that Dijkstra's algorithm does not necessarily work if some of the costs are negative by finding a digraph with negative costs (but no negative cost dicircuits) for whi

Variation and proportion, i am not getting what miss has taught us please w...

i am not getting what miss has taught us please will you will help me in my studies

Logarithmic differentiation, Logarithmic Differentiation : There is one...

Logarithmic Differentiation : There is one final topic to discuss in this section. Taking derivatives of some complicated functions can be simplified by using logarithms.  It i

The ratio of boys to girls at the dance was 3:4, The ratio of boys to girls...

The ratio of boys to girls at the dance was 3:4. There were 60 girls at the dance. How many boys were at the dance? Use a proportion comparing boys to girls at the dance. Boys/

find out the dimensions which will minimize, We desire to construct a box ...

We desire to construct a box whose base length is three times the base width. The material utilized to build the top & bottom cost $10/ft 2 and the material utilized to build the

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