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

Draw a common graph f ( x ) = |x|, Graph f ( x ) = |x| Solution The...

Graph f ( x ) = |x| Solution There actually isn't much to in this problem outside of reminding ourselves of what absolute value is. Remember again that the absolute value f

Surface area, Find the amount of sheet metal need to form a conical funnel ...

Find the amount of sheet metal need to form a conical funnel of base radius 30cm with a vertical height of 50cm, allowing for 0.5cm overlap. Find the total surface area?

gauss elimination method , Question: Use  Gauss elimination method to ...

Question: Use  Gauss elimination method to solve the following system of equations.  -y +3z=4  2x-y-2z= 2  2x-2y+z =6  4x-y-7z= 0

Sets, creative assignment about sets

creative assignment about sets

Maximax method-decision making under uncertainty, MAXIMAX method Maxima...

MAXIMAX method Maximax method is based upon 'extreme optimism' the decision maker chooses that particular strategy which corresponds to the maximum of the maximum pay off for e

Utilizes second derivative test to classify critical point, Utilizes the se...

Utilizes the second derivative test to classify the critical points of the function,                                               h ( x ) = 3x 5 - 5x 3 + 3 Solution T

Find out a series solution for differential equation, Find out a series sol...

Find out a series solution for the following differential equation about x 0 = 0 y′′ + y = 0.   Solution Note that in this case p(x)=1 and therefore every point is an or

Fractions, #how do I add fractions?

#how do I add fractions?

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