# Question:plot vectors

## Question:plot vectors

Maple
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.
﻿