## 115 Reputation

4 years, 211 days

## @Mariusz Iwaniuk  thank you soooo ...

thank you soooo mach

## @Mariusz Iwaniuk  than you very ma...

than you very mach .

## 0<x^2-y^2-4x<0.001i dont want to s...

0<x^2-y^2+4x<0.001

i dont want to solution exact 0

## @Rouben Rostamian thank you very very ...

@Rouben Rostamian

thank you very very mach

## @tomleslie  thanks for answer. i c...

thanks for answer. i correction my mistake. excuse me for my mistakes

i get error

restart;with(LinearAlgebra);Digits := 50;with(plots);

K := Matrix(5, 5, [[1.82780000000*10^7, 0., 0., 0., 0.], [0., 8.4785000000*10^6, 0., 0., 0.], [0., 0., 3.8748000000*10^6, 0., 0.], [0., 0., 0., 4.947100000000*10^8, 0.], [0., 0., 0., 0., 7.7075000000*10^6]]);
M := Matrix(5, 5, [[4.1197000000*10^6, 0., 0., 0., 0.], [0., 175900.0000, 0., 0., 0.], [0., 0., 52796., 0., 0.], [0., 0., 0., 2.0029000000*10^6, 0.], [0., 0., 0., 0., 21711.]]);
L := Vector[column](5, [400689.48, 36882.21, 24223.175, 7197.4865, 5007.6466]);
V := Vector[column](5, [4.4368, 48.201, 73.39160, 246.999, 355.012]);
ic := [seq({x[i](0) = 0, (D(x[i]))(0) = 0}, i = 1 .. 6, 1)];
eq := seq(diff(x[i](t), t\$2)+(2*0.5e-1)*sqrt(V[i])*(diff(x[i](t), t))+V[i]*x[i](t) = L[i]/M[i, i], i = 1 .. 5, 1);
res := evalf(seq(dsolve({eq[i], ic[i][]}, numeric, method = rkf45), i = 1 .. 5, 1));
cols := [black, blue, red, green, cyan];display([seq(odeplot(res[k], [t, x[k](t)], t = 0 .. 4, color = cols[k]), k = 1 .. 5)]);
Z := `<,>`(res[1], res[2], res[3], res[4], res[5]);
Phi := Matrix(5, 5, [[1., 1., 1., 1., 1.], [1.96506022575420, 1.62041320563413, 1.42204109823483, 0.548765310483432e-1, -.795724394004826], [9.21473910149806, 2.36597630710861, -.837387679777313, -8.23805381555997, 0.930594200639915e-1], [10.5672146719479, 1.26937656539014, -.710851068914949, 12.5801876194058, -0.564138989400088e-1], [12.4282433098880, -2.02516824673162, .481291188486771, -1.71479964168293, 0.513709374343217e-2]]);
Y := Phi . Z;
#odeplot(rhs(Y[1]), t = 0 .. 4)

## @Rouben Rostamian   this is a exam...

this is a example of finite element method . i want to writte a code in maple to solve faster this examples.

yes i mistake for a negative sign

ok. thanks a lot

## i draw it. i want to get a num...

i draw it.
i want to get a numeric of this plot for draw in excel . can i ?

## @vv  thank you can i transfer to ...

thank you

can i transfer to this picture

## @kambiz1199  ...

 1 2 3 4 Page 3 of 4
﻿