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

I am bad at math, i dont know how to do probobility iam so bad at it

i dont know how to do probobility iam so bad at it

#title.simpal harmonic motion., #questionShow that the system oscillates in...

#questionShow that the system oscillates in simple harmonic motion demonstrated by; , for which the general solution where X = (x – x0)..

Characteristics of exponential smoothing, Characteristics of Exponential Sm...

Characteristics of Exponential Smoothing 1. More weight is described to the most recent data. 2. All past data are incorporated not like in moving averages. 3. Les

Theory of quadratic equations.., solve the following simultaneous equations...

solve the following simultaneous equations x+y=a+b ; a/x_b/y

Fractions, What is two-thirds plus two-thirds?

What is two-thirds plus two-thirds?

Find out the dimensions of the field-optimization, We have to enclose a fie...

We have to enclose a field along with a fence. We contain 500 feet of fencing material & a building is on one side of the field & thus won't require any fencing.  Find out the dime

Disjointed sets or mutually exclusive, Disjointed Sets or Mutually Exclusiv...

Disjointed Sets or Mutually Exclusive Two sets are said to be mutually or disjointed exclusive whether they have no elements in common. Sets P and R underneath are disjointed

+, what is 2+2=

what is 2+2=

Explain the vertex formula, Explain the Vertex Formula ? The vertex for...

Explain the Vertex Formula ? The vertex formula is a convenient way of finding the vertex of the graph for any quadratic function. The graph of the quadratic equation f(x) = ax

find the present age, 5 years however, a man's age will be 3times his son'...

5 years however, a man's age will be 3times his son's age and 5 years ago, he was 7 times as old as his son.    Find their present ages.

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