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.

Write Your Message!

Name

Email id

Message

Verfication Code

User Account

All Pages

Initial value problems, math, Mathematics

Related Discussions:- Initial value problems, math

Ronding off numbers, how to round off numbers to the nearest tens and to th...

Calculate area of a square, The area of a square is given by the formula wi...

Derivation, d^2y/dx^2 if x=ct,y=c/t

Multiplying/dividing fractions and mixed numbers, -9/5 / 2

Definition of vertical asymptote, Vertical asymptote Definition : The funct...

#title., fixed cost of $1400 ,printing cost of .40 cents -each item to sell...

INVESTING MONEY, HOW MANY SHARES CAN I BUY WITH 1000 DOLLARS

Marketing, What are the Input and Output of Marketing

Ann, What was last years salary if after a 3% increase the salary is 35,020...

ALgebra, Please quote me a price

Write Your Message!

Featured Services

Online Tutoring

Project Development

Exam Preparation

Course Help

Assignment Help

Popular Subjects

Assured A++ Grade

Academics

Major Subjects

Majors

Get In Touch

TERMS & POLICIES

HELP & SUPPORT