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

Trig identities, What is the exact vale of sin(theta/2) when sintheta=3/5, ...

What is the exact vale of sin(theta/2) when sintheta=3/5, pi/2

Pair of straight line, The equation ax2 + 2hxy + by2 =0 represents a pair o...

The equation ax2 + 2hxy + by2 =0 represents a pair of straight lines passing through the origin and its angle is tan q = ±2root under h2-ab/(a+b) and even the eqn ax2+2hxy+by2+2gx+

Determine the equation of the line, Example :  Determine the equation of th...

Example :  Determine the equation of the line which passes through the point (8, 2) and is, parallel to the line given by 10 y+ 3x = -2 Solution In both of parts we are goi

Relative measures of dispersion-illustration, Illustration 2 In a ...

Illustration 2 In a described farm located in the UK the average salary of the employees is £ 3500 along with a standard deviation of £150 The similar firm has a local

Calcilate the height of the cone of which the bucket , A bucket of height 8...

A bucket of height 8 cm and made up of copper sheet is in the form of frustum of right circular cone with radii of its lower and upper ends as 3 cm and 9 cm respectively. Calculate

Find the area of shaded region of circle of radius, Find the area of shaded...

Find the area of shaded region of circle of radius =7cm, if ∠AOB=70 o , ∠COD=50 o and ∠EOF=60 o . (Ans:77cm 2 ) Ans:    Ar( Sector AOB + Sector COD + Sector OEF) =  7

The fisher’s index, The Fisher's index The index of Fisher acts as a c...

The Fisher's index The index of Fisher acts as a compromise between Paasche' index and Laspeyre's index. This is calculated as a geometric mean of the two indexes.

Evaluate infinity limit into the polynomial , Example   Evaluate following...

Example   Evaluate following limits. Solution Here our first thought is probably to just "plug" infinity into the polynomial & "evaluate" every term to finds out the

Parabola, please i need the answers to x^_7x+10 i want the vertex,axis of s...

please i need the answers to x^_7x+10 i want the vertex,axis of semetery,y intersect and the x intercept

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