## 150 Reputation

5 years, 229 days

## i want to solve a equation with newton m...

Maple

hi . i want to solve a equation with newton method but i cant,

restart;
f(x):=x^2-3;
df(x):=diff(f(x),x);
x:=1;
for i from 2 to 5 by 1 do
x:=x-(f(x)/df(x));
print(x);
end do

## i have problem in a solve example,int fo...

Maple

hi

i want to solve a example in this picture but i get error.

thanks

1.mw

example

## how can i simplify with definition two n...

Maple

how can i simplify with definition two new parameter

i show in this picture

## how can i solve with ShowSolution...

Maple

restart;
with(Student[Calculus1]);
a := Int(exp(cos(x)), y = 0 .. sin(x), x = 0 .. (1/2)*Pi);
ShowSolution(a);
i get error

## how can i draw a plot rkf45 ...

Maple

i want to draw plot Y[1] . but i cant. thank you

restart;
with(LinearAlgebra); Digits := 15; with(plots); with(Optimization);;
NewM := Matrix(5, 5, [[4119700.0000, 0., 0., 0., 0.], [0., 175900.0000, 0., 0., 0.], [0., 0., 52796., 0., 0.], [0., 0., 0., 2002900.0000, 0.], [0., 0., 0., 0., 21711.]]); NewK := Matrix(5, 5, [[18278000.0000, 0., 0., 0., 0.], [0., 8478500.0000, 0., 0., 0.], [0., 0., 3874800.0000, 0., 0.], [0., 0., 0., 494710000.0000, 0.], [0., 0., 0., 0., 7707500.0000]]); L := Vector[column](5, [400689.480747934, 36882.2103608425, 24223.1756570268, 7197.48654698287, 5007.64668721342]); V := Vector[column](5, [4.43679672962542, 48.2014976972537, 73.3916065733549, 246.999805581163, 355.012621460930]);
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]]);
NS := ImportMatrix("D:/tabas.txt", datatype = float[8]);
t_NS := NS[() .. (), 1];
acc_NS := NS[() .. (), 2];
plot(NS, t = 0 .. 4.5, size = [800, 400]);
acc := unapply(CurveFitting:-Spline(t_NS, acc_NS/(9.81), t, degree = 1), t);
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]*acc(t)/NewM[i, i], i = 1 .. 5, 1);
ic := [seq({x[i](0) = 0, (D(x[i]))(0) = 0}, i = 1 .. 5, 1)];
sol := seq(dsolve({eq[i], ic[i][]}, numeric, output = listprocedure), i = 1 .. 5, 1);
Z := `<,>`(seq(sol[i][2], i = 1 .. 5));
#Y := Phi . Z;
#odeplot(Y[1], t = 0 .. 10);
tabas.txt

1.mw

 5 6 7 8 9 10 11 Page 7 of 11
﻿