## 10 Reputation

13 years, 209 days

## "ERROR " Crout matrix decomposition...

hi , can anyone tell me what's wrong with this!! please!

> crout := proc (A::array, n::integer)

local L, U, j;

L := matrix(n, n, 0);

U := matrix(n, n, 0);

L[1, 1] := A[1, 1];

L[2, 1] := A[2, 1];

U[1, 2] := A[1, 2]/L[1.1];

for j from 2 to n-1 do

## Crout matrix decomposition...

Maple

did someone know the program of crout matrix decomposition?!

## cholesky decomposition...

hi , can anyone tell me what's wrong witht his! please!!

> choleski:=proc(S::array,n::integer)
> local L,i,j,k,t:
> L:=matrix(n,n,0); # matrice nulle d’ordre n
> for j to n do # calcul de L, colonne par colonne

## Newton's interpolation polynomial...

hi , can anyone tell me what's wrong with this!

>  interpol:=proc(f)

local i,j,k,n,epsilon,x,y,b,p;

## "ERROR " Fixed point method...

hi , can anyone tell me what's wrong withthis!! plzzz!thank you in advance!

Fixed point method:

> fixe := proc (p, g)

local n, nmax, e;

e := 1/100000000; nmax := 500;

while (evalb(n < nmax) and e < abs(g(p)-p) do

n := 1;

p := g(p); print(n, "             ", evalf(p, 10...

 Page 1 of 1
﻿