## 25 Reputation

0 years, 318 days

 >
 >
 >

## @tomleslie thankyou! This Post...

@tomleslie thankyou! This Post is so helpful to me, really really helpful

## @tomleslie  i want the command parr...

@tomleslie  i want the command parrallel to solve equations...

system picks all equations automaiically after loop.

Maple Worksheet - Error

Failed to load the worksheet /maplenet/convert/model.mw .

## @tomleslie No. fiifteen equatiions ...

@tomleslie No. fiifteen equatiions fifteen unknowns..

## change function......

@acer if i change functtion
yfunc := proc (t) options operator, arrow; sqrt(202650*t*K/(phi*mu)) end proc

then how to write in form of matrix ??

Mfunc := Matrix(10, 10, proc (i, j) options operator, arrow; yfunc(a+(1/9)*(i-1)*(b-a)) end proc)

same like this or  somthing change ??? i can't understand how u write this form

explain please... i shall be very thankful to you

## u there???...

@acer explain please.  i need ur help.. thank u :)

@acer  thank you.. ur code is very helpfull for me..
can you expalin how you write Mfunc???
Mfunc := Matrix(10, 10, proc (i, j) options operator, arrow; yfunc(a+(1/9)*(i-1)*(b-a)) end proc)

thank you for your help :)

## how to get the values and write in form ...

@acer without equ[i,j] or seq cammond how to get the u[i,j] values and write in form of matrix?? Crank_scheme.mwCrank_scheme.mw
ans anminate with change in time k from 0.001..0.33??
and difference between function y plot...

## i explain what i want to do......

@acer
i want to find the difference between numerical values and experimental values...  for numerical values, discitization scheme use. my discritization scheme is this

```-lemda*u[i-1,j+1]+(2+2*lemda)*u[i,j+1]-lemda*u[i+1,j+1]=lemda*u[i-1,j]+(2-2*lemda)
*u[i,j]+lemda*u[i+1,j]```

estimate the value of u[i,j]  then write in form of matrix to get matrix plot.. i dont know is there an other way to get plot of u[i,j] and next i want to find chnage in time (k) koi k=0.0001 to 0.001... chnage in numerical values with change in time k...
and i change the data values into function y (time and position ) data..... and plot in 3d to find the difference.. here  explain what i want...
hope u understard what i want to do.. is there another way to plot discritized scheme then plz do favour..
your work error plot is not helpful for me....  plz suggest some other way to understand...... without equ[i,j] or [seq] can i find thevalue of u[i,j] and plot graph simply and find difernce .. minor error ??
hope you understand what i want... so plz see n check more what we do to simplify code ?? and get minnor difference

## help to solve problem...

@acer please answer..!! how to fin difference between matrix plot and function plot???

## how to plot difference/error graph of m...

@acer plz do one more help .. how to find the differnce between 2 graphs . i want to estimate error. i m a little bit confuse how to minus matrix plot into functioon plot ???

 >
 >
 >
 (1)
 >
 >
 >
 >
 >
 >
 >
 >
 > sys := ([seq])(seq(eq[i, j], j = 0 .. N), i = 1 .. N):
 > nops(sys);
 > vars:=indets(sys) minus {k}:
 > nn := Matrix(N+1, N+1,(i, j)-> u[i-1, j-1]):
 > ##
 > p:=proc(kk) local u_res,A;
 > u_res:=solve(eval(sys,k=kk),vars);
 > A:=eval(nn,u_res);
 > #plots:-matrixplot(A)
 > end proc:
 (2)
 > ## Testing p for k=0.001:
 > plots:-matrixplot(p(0.0001)-p(0.0009));
 >
 >
 >
 (3)
 >
 >
 >
 >

## chnge the No of points N:= 49 system ta...

@acer for N= 49 my system take 24 hours  and no results show... can u help me for geting  results plz

 >
 >
 >
 (1)
 >
 >
 >
 >
 >
 >
 > sys := ([seq])(seq(eq[i, j], j = 0 .. N), i = 1 .. N):
 > nops(sys);
 > vars:=indets(sys) minus {k}:
 > nn := Matrix(N+1, N+1,(i, j)-> u[i-1, j-1]):
 > ##
 > p:=proc(kk) local u_res,A;
 > u_res:=solve(eval(sys,k=kk),vars);
 > A:=eval(nn,u_res);
 > #plots:-matrixplot(A)
 > end proc:
 (2)
 > ## Testing p for k=0.001:
 > plots:-matrixplot(p(0.001)-p(0.0005));
 > ## Animating the plot for k=0.0001..0.001:
 > plots:-animate(K->plots:-matrixplot(p(k)-p(K)),[k],k=0.0001..0.001);
 >

## how to plot in 2d ???...

@Carl Love can we plot this graph in 2d..???

 >
 (1)
 >
 (2)
 >
 >
 >
 >
 >
 (3)
 >
 >
 (4)
 >
 >