Advanced Search

Question: plot vectors

Question:plot vectors

dkrbabajee 49 Maple
vector plotting
December 28 2007
0
Hi I am getting a problem with plot. I am trying to plot ka (number of iterations) versus va (error norm) in this worksheet: > with(linalg):with(LinearAlgebra):with(plots):SOS:=proc(x) local n,k; n:=vecdim(x); sum(x[k]^2,k=1..n); end;NJ2Next:=proc(f,x,x0) local val,jac, xn, xn1, val1, val2; val:=eval(subs(x=x0,f)); jac:=eval(jacobian(f,x),x=x0):jac:=inverse(jac); xn:=evalm(x0-jac&*val); val1:=eval(subs(x=xn,f));xn1:=evalm(xn-jac&*val1); val2:=eval(subs(x=xn1,f)); evalm(xn1-jac&*val2); end: > NJ2:=proc(f,x,x0,N) local k, sol, valf, va, ka;sol:=array(0..N):sol[0]:=x0:va:=array(1..N):ka:=array(1..N):valf:=eval(subs(x=x0,f)): for k from 1 to N do sol[k]:=NJ2Next(f,x,sol[k-1]): valf:=eval(subs(x=sol[k-1],f)): va[k]:=evalm(log(norm(sol[k]-sol[k-1],2))):ka[k]:=eval(k):err:=evalm(norm(sol[k]-sol[k-1],2)): od: print(ka):print(va) :plot([ka,va]):end: > x:=array(1..2):Digits:=100;F1:=[x[1]^2-2*x[1]-x[2]+0.5,x[1]^2+4*x[2]^2-4];X1:=[-.3,1.]; Digits := 100 F1 := [x[1] - 2 x[1] - x[2] + 0.5, x[1] + 4 x[2] - 4] X1 := [-0.3, 1.] > st:= time(): NJ2(F1,x,X1,4): time() - st; [1, 2, 3, 4] [-2.55078506399648359592905174620552185528060918024541977391565768568635516\ 8364202167048764106251232207, -11.39740819108762660117294666038928030717478\ 929221870097017112980886449952632270964551725126295769979, -46.584115598915\ 422739808075169077071195792568837243840280671447846182066678716385354809725\ 54796113013, -187.345668309705534371600275407273119332095933488131206414042\ 3327883190073745036498254923506929068150] Warning, unable to evaluate the functions to numeric values in the region; see the plotting command's help page to ensure the calling sequence is correct The above message is displayed. Can anyone help me to obtain the plot of ka versus va? Thanks. I need it for my thesis.
Please Wait...