folsen

function bvp( N )

r0 = 0.02; %m
h = 20; %W/m^2*K
k = 14; %W/m*K
TA = 20; %degrees C
T0 = 100; %degrees C
L = 0.1; %m

a = @(x) -1./(L-x/2);

folsen

/*
* AppController.j
* NewApplication
*
* Created by You on July 5, 2009.
* Copyright 2009, Your Company All rights reserved.
*/

@import <Foundation/CPObject.j>
@import "FileUpload.j"

folsen

function errVSstep(A, y0, t0, tf)
errors = [];
range = 1:20;
h = (tf-t0)./2.^range;

for i=range
    [approx, err] = eulerint(A,y0,t0,tf,2^i);
    errors = [errors, norm(err)];
end