# Question:solve numeric stiff the ODEs and plot

## Question:solve numeric stiff the ODEs and plot

Maple

Hi every body:

I have six ODEs in two case with their initial conditions:

case 1:

eqs1:= diff(r[1](t), t, t)+.68041478*(diff(r[1](t), t))+732.61274*r[1](t)-0.11751824e-2*V(t) = -0.52084940e-2*q(t), 2.03*10^(-8)*(diff(V(t), t))+4.0650406*10^(-7)*V(t)+0.11751824e-2*(diff(r[1](t), t)) = 0, diff(q(t), t, t)+(9.5199780*(q(t)^2-1))*(diff(q(t), t))+1006.9998*q(t) = -2192.5091*(diff(r[1](t), t, t));
ics1 := V(0) = 0, q(0) = 0.1e-2, (D(q))(0) = 0, r[1](0) = 0, (D(r[1]))(0) = 0;

case 2:

eqs2:= diff(r[2](t), t, t)+3.3737665*(diff(r[2](t), t))+52672.540*r[2](t)-0.93177795e-2*V(t) = 0.18358493e-2*q(t), 2.03*10^(-8)*(diff(V(t), t))+4.0650406*10^(-7)*V(t)+0.11751824e-2*(diff(r[1](t), t))+0.93177795e-2*(diff(r[2](t), t)) = 0, diff(q(t), t, t)+(9.5199780*(q(t)^2-1))*(diff(q(t), t))+1006.9998*q(t) = -2192.5091*(diff(r[1](t), t, t))+772.79848*(diff(r[2](t), t, t));

ics2 := V(0) = 0, q(0) = 0.1e-2, (D(q))(0) = 0, r[2](0) = 0, (D(r[2]))(0) = 0;

How can i solve numeric stiff the ODEs and plot V(t) in case 2 for t=0..2?

tnx...

﻿